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1 . INTRODUCTION 


1.1  SCOPE  OF  THE  REXOR  PROGRAM 

REXOR  is  a rotorc’-aft  analysis  tool  which  has  resulted  from  applying  an 
interdisciplinary  math  modeling  philosophy.  Reference  1.  The  REXOR  math 
model  has  been  written  for  a single  four-bladed,  gyro-controlled, 
hingeless-rntor  helicopter  with  additional  capability  for  analysis  of 
teetering  or  hinge-offset  rotor  systems  with  conventional  controls  and  two 
or  four  blades.  This  helicopter  may  be  conventional  in  design,  winged,  or 
compounded.  The  specific  analysis  is  limited  to  a maximum  of  four  blades, 
but  it  can  be  expanded  to  include  more  blades  by  following  the  detailed 
mechanical  derivation  procedure  established  for  the  analysis . The  model  is 
broken  down  into  uhe  three  major  categories  shown  in  Figure  1-1.  These 
categories  are  the  control  system,  the  rotor,  and  the  body. 


Figure  1-1  indicates  the  manner  in  which  these  components  are  related  to 
one  another  as  utilized  in  the  analysis.  The  analysis  is  the  simulation  of 
an  entire  aircraft,  which  includes  a detailed  dynamic  description  of  the 
rotor  and  control  system  as  well  as  a conventional  six-degree-of- freedom 
body  dynamic  description  which  operates  in  two  modes  identified  as  TRIM 
and  FLY.  In  the  TRIM  mode,  the  aircraft  is  constrained  to  a prescribed 
static  flight  condition  while  the  controls  are  activated  and  the  rotor  is 
allowed  to  respond  to  obtain  a force  and  moment  equilibrium  of  the  aircraft 
at  that  static  condition.  In  the  FLY  mode  the  entire  aircraft  is  free  to 
respond  dynamically  to  a maximum  of  30  degrees  of  freedom  to  control  inputs 
or  to  any  other  arbitrary  inputs  such  as  gusts.  Pilot  inputs  can  be  any 
single  or  multiple  control  manipulation  in  the  form  of  simple  steps  or 
pulses,  doublets,  stick  stirs,  or  other  transient  input  within  the  capa- 
bilities of  the  control  system  simulated.  As  a result,  transient  loads  and 
resulting  aircraft  and  rotor  dynamic  response  can  be  obtained.  For  corre- 
lation purposes,  actual  flight  test  control  motions  can  be  used  as  input 
to  provide  comparative  response  data..  For  specialized  applications,  an 
analytic  autopilot  may  be  used  to  control  the  flight  path  of  the  aircraft. 
Additionally,  gust  inputs  and  other  types  of  external  excitations  could 
be  applied  directly  to  the  rotor  and/or  airframe. 


1.2  REXOR  CAPABILITIES 

REXOR  is  a detailed  rotorcraft  math  model  simulation  with  particular 
emphasis  on  the  main  rotor  mechanics.  The  program  is  particularly  valuable 
in  a detailed  exploration  of  rotor  characteristics  of  proposed  uesigns,  in 
identifying  problem  areas  and  verifying  fixes  in  flight  test  development 
programs.  A case  history  is  given  in  Reference  2. 
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Typical  REXOR  applications  are  listed  below. 

Dynamics : 

• Rotor  stability  in  regions  below  Ip  as  function  of  flight 
speed,  maneuvers,  rotor  rpm,  nonlinear  blade  aerodynamics 

• of^0^  f nSr°iVity  and  diss 'Pation  capacity  as  a function 
of  gusts  and  pilot  control  inputs 

• Effects  of  design  parameters  (mechanical  and  elastic 
couplings,  controls,  etc.)  on  rotor  stability  and  load 
sensitivity 

• Correlation  and  check  of  specialized  dynamic  models. 

Handling  Qualities : 

• Vehicle  response  to  pilot  control  inputs  for  vehicle  flight 
conditions,  speed,  altitude,  rotor  rpm,  design  parameter 
variations 

• Vehicle  stability  as  function  of  speed,  rotor  rpm,  flight 
conditions,  design  parameters 

• Effect  of  design  parameter  variations  on  handling  qualities 

• Development  and  checking  of  handling  qualities  models. 

Failure  Analysis: 

• Effect  of  loss  of  one  inplane  damper  on  subsequent  flight 
time  history 


• Blade  projectile  hit  and  ensuing  events 

* Blade  strike  and  resulting  rotor  track. 


Performance 


• Correlation  and  independent  check  of  performance  models, 
particularly  in  regions  of  highly  nonlinear  blade  aerodynamic 
operation  (retreating  blade  stall  and  compressibility  effects) 

• Develop  data  for  performance  models  for  use  in  nonlinear 
areas 


"rx  * 


E 


I 


Loads : 


Steady-state  rotor  loads  as  a function  of  rotor  rpm,  flight 
velocity,  control  trim  settings 


Dynamic  rotor  loads  as  a function  of  rotor  rpm,  flight  velo 
city,  vehicle  maneuvers,  pilot  control  inputs 


Rotor/fuselage  clearances  as  a function  of  speed,  vehicle 
maneuvers,  rotor  rpm,  pilot  control  inputs,  flight 
configuration 


Rotor/ fuselage /wing  design  characteristics  requirements  as 
functions  of  maneuver  load  factor,  control  commands  (see 
Reference  3). 


1.3  IMPROPER  APPLICATION  OF  REXOR 


While  REXOR  is  capable  of  performing  a number  of  analysis  tasks,  the 
program  range  of  use  is  certainly  not  all  inclusive.  Examples  of  types 
of  use  where  REXOR  either  wouldn't  work  well  or  would  be  impractical  are 
given  below. 


REXOR  is  an  extensive  math  model  and,  as  such,  may  consume  a considerable 
amound  of  computer  time  to  execute  a case.  Therefore,  the  program  is  not 
intended  as  a parametric  design  analysis  tool,  but  rather  as  a device  to 
verify  the  correctness  of  a parametric  selection  process. 


REXOR  does  not  treat  blade-to-blade  vortex  interaction,  ^is  condition 
limits  the  validity  of  the  vibration  solution  in  the  transition  flight 
regime . 


REXOR  typically  uses  twenty  or  less  blade  radial  stations.  The  computer 
blade  deflections  show  good  correlation  to  measured  data  with  this  model- 
ing. However,  since  shear  is  a first  derivative,  and  moment  is  a second 
derivative  of  deflection  data,  care  needs  to  be  exercised  in  their  use 
(Reference  4). 


1.4  THE  REXOR  REPORT  AND  ITS  USE 


This  report  is  presented  in  three  volumes. 
• Volume  I 


A development  of  rotorcraft  mechanics  and  aerodynamics  including 
a derivation  of  the  equations  of  motion  from  first  principles. 
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e Volume  II 


V 


The  development  and  explanation  of  the  computer  code  required  to 
implement  the  equations  of  motion. 

0 Volume  III 

A user's  manual  containing  a description  of  code  input/output  and 
instructions  to  operate  the  program. 

Volume  I is  intended  to  be  a self-sufficient  guide  to  the  math  development 
of  the  equations  of  motion  and  is  the  reference  background  as  duch.  Volume 
II  gives  the  location  of  computation  elements , and  serves  to  locate  elements 
for  inspection  or  modification.  Volume  III  presents  normal  program  operation 
plus  troubleshooting  guide  material  required  for  day-to-day  program  use. 


I 


i 

i 


2 . BACKGROUND 


THE  DERIVATION  OF  REXOR 


The  REXOR  analysis  is  the  result  of  a natural  evaluation  of  capability 
arising  from  a desire  to  develop  a rotorcraft  handling  qualities  e.nalysis. 
This  required  a model  of  a complete  aircraft  which,  as  it  was  refined, 
exhibited  unique  advantages  for  analysis  of  dynamic  stability,  loads  and 
performance.  The  eventual  outgrowth  is  an  interdisciplinary  analysis  with 
broad  areas  of  application. 


The  initial  development  effort  was  called  Rotor  Junior  (Reference  5)-  This 
program  had  a reasonably  detailed  nonrotating  airframe  representation,  but 
used  a steady-state  representation  for  the  main  and  tail  rotors.  (Bailey  T 
Coefficients,  Reference  6.)  A major  rework  of  this  program  especially  in 
the  main  rotor  led  to  Rotor  Senior.  Here  the  dynamic  response  of  rotating 
blades  was  calculated  (References  7 and  8).  Many  simplifying  assumptions 
were  made  in  the  construction  of  the  blade  dynamics.  Starting  afresh  in 
this  area  plus  a general  overhaul  of  the  model  by  an  interdisciplinary  team 
of  specialists  led  to  Revised  and  Extended  RotOR  Senior  or  REXOR.  A key 
feature  was  incorporating  the  solution  of  generalized  force  equations  for 
incremental  accelerations  using  a generalized  mass  matrix,  and  thus 
coupling  of  all  the  degrees  of  freedom. 


GROWTH  OF  CAPABILITY 


The  original  model  (Rotor  Junior)  was  intended  strictly  for  stability 
derivatives.  With  the  addition  of  blade  dynamics,  some  evidence  of  and 
useful  ir'.ight  into  subharmonic  rotor  instability  problems  was  obtained 
(Rotor  Periior).  With  the  advent  of  very  detailed  blade  and  blade  inter- 
face modeling,  many  new  design  problem  areas  have  been  avoided,  and 
deficiencies  of  existing  designs  corrected  (REXOR). 
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BASIC  COMPUTATIONAL  IDEA 


3.1  MODAL  SOLUTION  - OVERVIEW 

The  aircraft  is  described  dynamically  in  approximately  thirty  fully-coupled 
degrees  of  freedom.  In  addition  to  the  six  degrees  of  freedom  of  the  hub 
principal  reference  axes,  the  fuselage  is  related  to  the  rotor  hub  through 
rotor  rotational  speed,  shaft  pitch,  and  shaft  roll  degrees  of  freedom. 

The  control  gyro/swashplate  combination  has  three  degrees  of  freedom. 

Motion  of  each  of  the  four  main  rotor  blades  is  described  by  three  coupled 
flapwise  and  inplane  modes  and  a pitch  horn  bending  degree  of  freedom 
which  couples  blade  feathering  to  the  control  gyro  or  swashplate.  These 
four  degrees  of  freedom  per  blade  make  sixteen  for  the  four  blades. 


The  blade  modes  are  primitive  moues  in  that  they  are  determined  from  a 
lumped  parameter  analysis  of  a rotating  cantilever  blade  at  a selected 
rotor  speed  and  collective  blade  angle,  hereafter  referred  to  as  the 
reference  feather  angle.  The  generalized  stiffness  matrix  is  computed 
using  these  rotating  modes  and  contains  only  the  structural  stiffness  of 
the  blades  and  hub.  This  formulation  ensures  proper  internal  and  external 
force  and  moment  balance.  The  modal  deflections  outboard  of  the  feather 
hinge  are  rotated  through  tne  actual  feather  angle  less  the  reference 
feather  angle.  Thus,  blade  element  deflections  outboard  of  the  feathering 
hinge  due  to  modal  displacements  are  defined  to  remain  aligned  with  a 
coordinate  axis  system  which  is  orthogonal  to  a plane  containing  the 
instantaneous  deformed  feather  axis  and  rotated  tnrough  the  instantaneous 
feather  angle  less  the  reference  feather  angle.  As  a result,  the  internal 
strain  energy  in  the  blade  due  to  unit  modal  displacements  is  invariant 
with  variation  in  blade  angle.  This  technique  permits  the  highest  resolu- 
tion of  motion  and  forces  for  the  blade  with  ar.  assumed  mode  solution  for 
a given  number  of  modes. 


3.2  ENERGY  METHODS  DEVELOPMENT 

The  equations  of  motion  for  REXOR  are  developed  from  Lagrange's  equations, 
which  is  an  energy  approach.  If  one  can  express  the  kinetic,  potential, 
and  dissipative  energies  of  a system  in  addition  to  the  work  done  by  exter- 
nal forces,  then  Lagrange's  equations  provide  a powerful  method  foi  devel- 
oping the  equations  of  motion. 

The  dynamic  equations  of  motion  are  written  in  matrix  form  as 

- [A]  {q  ( + |g(=  0 (3-1) 
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where  [a]  is  approximately  a 30  by  30  matrix  of  generalized  mass  ele- 
ments, { o.  {■  is  a column  matrix  of  accelerations  of  the  generalized  coor- 
dinates and  | G | is  a column  matrix  derived  from  the  Lagrangian  energy 

functions,  dissipation  function  and  generalized  forces,  which  take  the 
form: 


|g|  = - [b]U  |-  - [c]  hi*  [Q]  {f  (t)( 


(3-2) 


The  equations  of  motion  are  solved  as  a time  history  at  rotor  azimuth  angle 
increments  required  to  provide  a stable  solution  for  the  highest  frequency 
mode  present  in  the  solution. 


ft: 


i 
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3-3  CALCULATION  OF  ROTOR  MODE  DISPLACEMENTS , VELOCITIES  AND  ACCELERATIONS 

In  a rotor  simulation  of  this  type,  it  is  difficult  to  compute  the  proper 
displacement  velocities  and  accelerations,  and  associated  inertia  and 
aerodynamic  forces  and  moments  which  are  required  for  high  resolution  of 
the  blade  feathering  moments.  This  requires  exacting  aerodynamic  data  as 
well  as  a precise  statement  of  the  inertial  loadings.  To  establish  the 
feathering  moments  due  to  these  loads,  the  relationship  between  the  feather 
axis  and  the  point  of  application  of  the  loads  must  be  precisely  determined. 
This  is  accomplished  by  a very  accurate  analytic  constructioi  of  the  unde- 
formed blade  and  a superposition  of  the  blade  elastic  bending  on  this 
shape.  In  order  to  achieve  the  highest  resolution  of  the  predicted  blade 
shape  and  feather  axis  position,  the  blade  modes  are  defined  it  approxi- 
mately the  trim  collective  blade  angle.  The  blade  static  position  is  also 
constructed  at  this  blade  angle.  Blade  element  displacements,  velocities, 
and  accelerations  are  then  computed  from  the  combined  static  shape,  the 
elastic  blade  motion,  and  blade  feathering  with  respect  to  the  reference 
feather  angle. 

The  aerodynamic  description  used  in  the  analysis  is  composed  of  a rotor 
inflow  model,  nonlinear  steady  and  unsteady  blade  element  aerodynamics, 
nonlinear  body  aerodynamic  characteristics,  rotor /body  aerodynamic 
interference,  and  auxiliary  airloads  from  the  tail  rotor  and  propeller. 

The  auxiliary  airloads  are  contained  in  modular  subroutines  and  are  func- 
tions of  advance  ratio  and  propeller  and  tail  rotor  collective  pitch.  The 
main  rotor  dovnwash  effect  on  the  wing  and  horizontal  tail  angles  of  attack 
is  an  empirical  function  of  rotor  thrus;  and  advance  ratio.  The  nonlinear 
body  aerodynamics  may  be  inputted  as  tables  of  actual  wind  tunnel  test 
data. 

The  aircraft  primary  control  systems  are  simulated  from  the  pilot  control 
levers  operating  through  a boost  system  in  all  control  axes . Gearing  and 
gains  in  the  control  path  are  inputs  to  the  analysis  and  may  be  easily 
changed  for  studying  the  effects  of  design  changes  in  the  control  system. 
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Control  servos  are  simulated  by  first-order  lags  with  rate  limits  and  with 
soft  and  hard  physical  stops.  Control  stiffnesses  in  collective  and  cyclic 
pitch  axes  of  the  main  rotor  are  included  in  the  dynamic  equations  of 
motion. 


3.1+  OUTPUT 

The  analysis  is  a time-history  solution  of  the  equations  of  motion.  The 
standard  output  format  provides  plots  of  up  to  40  output  parameters  in  TRIM 
and  6 0 in  FLY.  In  addition,  a tabulation  is  provided  at  the  end  of  each 
mode  of  operation.  The  analysis  also  provides  plot  capability  at  the  end 
of  TRIM  to  show  loads  at  various  points  in  the  system  over  a single  rotor 
revolution  on  an  expanded  linear  scale.  These  loads  are  harmonically 
analyzed  and  the  harmonic  components  printed  out. 

In  addition  to  the  tabbed  or  plotted  time-history  output  in  the  FLY  mode, 
a tape  or  disc  pack  record  can  be  made  at  selected  time  points  for  use  in 
a fast  Fourier  transform  analysis.  The  record  pickup  interface  and  FFT 
are  not  supplied  as  part  of  REXOR,  but  are  included  by  the  user  to  match 
the  computation  facility  used.  The  FFT  analysis  gives  a quantitative 
evaluation  of  the  mode /member  oscillation  in  question.  Accurate  damping 
ratio  and  resonant  frequency  data  are  available  from  this  procedure. 
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4 . SYMBOLS 


The  notation  used  in  REXOR  generally  follows  what  could  be  termed  NASA 
notation.  In  general: 


Axis  systems  use  a right-hand  triad  X,  Y,  Z 


• Rotations  about  these  axes  are  also  a right -nand  triad  0,  <}>,  i|> 


Rotation  rates,  again  a right-hand  triad, are  p,  q,  r 


Velocity  components  of  X,  Y,  Z are  u,  v,  w. 


4.1  SUBSCRIPTING  NOTATION 


Subscripting  is  used  as  a rule  in  REXOR  to  further  identify  a variable. 
Superscripts  except  in  a few  column  vectors  are  reserved  to  denote  raising 
to  a power.  The  subscripting  can  mean: 


• Type  of  element;  F for  fuselage,  SP  for  swashplate,  TR  for 
tail  rotor,  R for  rotor,  etc. 


• Coordinate  system  reference;  BLn  for  blade  axis,  H for  hub  axis, 
R for  rotor  axis,  etc. 


Modal  identifiers, 


4.1.1  Blade  Number 


The  blade  modal  identifier  typically  is  of  the  form  A . Where  n is  the 

mn 


blade  num.be r. 


4.1.2  Mode  Number 


Also  from  A , m is  the  mode  number,  and  is  keyed  to  the  symbol  A.  A 
mn  - J 


represents  blade  bending  modes  (3).  Therefore  m can  be  1 to  3. 


4.1.3  Mode  Type 


Other  than  blade  bending  the  remaining  blade  mode  is  torsion,  and  is  sepa- 


rately identified  as  B„Tr  . Nonblade  modes  are  identified  by  the  direction 

PHn 


and  subscripted  axis  of  motion.  Examples  are  ^ r for  rotation  of  the 


rotor  and  9 for  longitudinal  shaft  bending. 

Q 
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U.l.it  Generalized  Mass,  Damper,  Spring,  Forces 

The  generalized  masses  are  denoted  as  M doubly  subscripted  by  the  two 

modes  active  for  that  mass.  Examples  are  M „ „ and  M,  . . This 

* bS  $S  Amn  6„ 

n 

scheme  is  also  used  for  other  elements  of  the  equations  of  motion,  dampers 
(C),  springs  (K),  forces  (F).  Note  the  forces  are  a column  vector  and 
singly  subscripted. 

U.1.5  Forces  and  Moments 

In  the  process  of  forming  the  equations  of  motion  many  subelements  of 
forces  and  moments  are  formed,  translated  and  combined.  Several  layers 
of  subscripting  may  exist  in  performing  this  process.  The  guidelines  to 
the  layering  are: 

• First  layer  denotes  the  direction  or  axis  system  that  the  quantity 
is  formed  in.  Examples  are  X and  BLE. 

• Second  is  the  axis  system  involved  or  axis  system  being  translated 
to,  depending  on  the  specification  of  the  first  level.  The  second 
level  may  also  be  specified  as  0 or  nought,  to  indicate  the 
value  is  at  the  coordinate  system  origin.  This  notation  is  used 
to  show  an  inertial  reference  and  blade  root  summation  quantities. 

• The  third  layer,  usually  outside  a series  of  bracketed  quantities, 
shows  the  blade  number  being  computed,  or  the  overall  coordinate 
system  in  use  for  the  computation  at  hand. 
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5.  COORDINATE  SYSTEMS  AND  TRANSFORMATIONS 


5.1  INTRODUCTION 

Prior  to  developing  the  equations  of  motion,  a system  of  coordinate  sets 
with  a description  of  the  elements  of  the  system  in  these  sets  and  the 
interrelationship  of  the  sets  is  required. 


5.2  COORDINATE  SETS 

5.2.1  Fuselage  Coordinates  (X  , Y , Z„) 

r r r 

The  fuselage  set  origin  is  located  along  the  undeflected  main  rotor  shaft 
line,  and  at  a convenient  waterline  (the  eg  for  instance).  See  Figure  5-1. 

The  coordinates  form  a right-hand  triad  X^,  Yp,  Z„.  Notations  for  veloc- 
ities with  respect  to  earth  of  these  coordinates  are  either  Xp,  Yp, 
or  Up,  Vp,  Wp.  A conventional  double  dot  notation  is  used  for  accelera- 
tion. Euler  rotations  of  the  set  follow  conventional  practice  of  roll  right 
<l> p,  pitch  up  0p,  and  yaw  right  Rates  of  rotation  are  either 


denoted  by  dot  notation  or  p , q„,  r . Angular  acceleration  is  double 

r I*  f ...... 

dot  notation  of  the  rotation  or  dot  notation  of  the  rates,  <j>  , 0 , i|> 
or  pF,  qp,  rp. 

Numerous  aerodynamic  terms  are  referenced  to  the  fuselage  set.  Figure  5-2 
shows  the  relationship  of  airflow  to  this  set.  The  components  of  airflow, 

also  noted  as  u , v , w , are  defined  with  respect  to  the  fuselage  set 

r r r 

by  an  angle  of  attack  a,  and  a sideslip  angle  g.  The  angle  of  attack 
is  the  arcsin  of  the  ratio  of  the  vertical  component  and  the  vector  sum  of 
the  X and  Z components.  The  sideslip  is  the  Y component  of  airflow 
in  relation  to  the  total  vector  airflow  sum.  The  angle  of  attack  is 
positive  (pitch  up)  of  the  fuselage  set  with  respect  to  the  airflow.  The 
sideslip  is  positive  (yaw  left)  for  the  airflow  relative  to  the  set.  The 
airflow  is  the  vector  sum  of  the  fuselage  set  inertial  motion  and  flow 
fields  from  other  parts  of  the  vehicle,  such  as  main  rotor  downwash. 

5.2.2  Hub  Coordinates  (Principal  reference  set)  (X^,  Y^,  Z^) 

The  hub  set  origin  is  at  the  top  of  the  main  rotor  mast,  but  does  not 
rotate  with  the  mast.  The  mast  top  location  represents  the  optimum  choice 
as  a summing  point  for  loads,  and  a reference  point  to  track  relative  posi- 
tions of  model  elements  with  a minimum  of  algebraic  operations. 
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Figure  5 1.  Coordinate  Systems  Fuselage  Set 


For  an  undeflected  rotor  mast,  the  hub  set  origin  is  in  line  with  the 

axis  a distance  Z up  from  the  fuselage  set.  See  Figure  5-3.  At  this 
Or 

point,  the  hub  set  X,  Y axes  are  in  the  same  direction  and  sense  as  the 
fuselage  set  X,  Y. 

Due  to  small  angle  assumption  shaft  (s)  set  bending,  the  hub  set  is  rotated 
and  6g  (Euler  angles)  from  the  fuselage  set.  As  the  hub  set  is  the  prin- 
cipal reference  axis,  the  fuselage  set  rotates  under  the  hub  set.  The  fuselage 
set  origin  moves  to  the  right  and  back  for  positive  and  e 0 . The  corr=?3pon- 
ding  translations  are  (3Xp/38g)  6'  and  (3Yp/34>g)  4> g . The  partials  are  con- 
stants, and  are  zero  for  a virtual  pivot  point  at  the  fuselage  set  origin. 

Airflow  information  is  also  referenced  to  the  hub  set  for  use  in  the  main 
rotor  aerodynamic  calculations.  The  reference  scheme  is  shown  on 

Figure  5-4.  For  components  of  airflow  u , v , w with  respect  to  the 

hub  set,  an  angle  of  attack  cx^,  and  sideslip  ip  are  defined.  The 

generation  conventions  are  different  from  the  fuselage  airflow  reference. 

5.2.3  Rotor  Coordinates  (X„,  Y„,  Z„) 

R R R 

The  undeflected  rotor  set  has  the  same  origin  as  the  hub  set  and  a common 
Z axis.  See  Figure  5-5-  However,  the  X axis  rotates  with  the  blade 

A 

number  1 reference  axis  system.  The  YR  axis  points  to  the  blade  num- 
ber 1 leading  edge.  The  rotation  of  the  rotor  set  is  measured  counter- 
clockwise (CCW)  from  the  -X  axis  by  the  angle  . The  rowor  set 

n R 

serves  as  a datum  basis  for  blade  number  one  and  has  a common  Z axis 
with  the  BLn  sets. 

5.2.4  Blade _ Coordinates  (X^,  Y^,  Z^) 

To  bookkeep  the  deflections  properly  of  all  the  main  rotor  blades,  sets 
equivalent  to  the  rotor  set  are  created  for  each  blade.  These  are  the 
BLn  sets,  where  n is  the  blade  number  (counted  clockwise  from  blade 
number  one).  All  BLn  sets  are  identical  except  for  an  azimuthal  rota- 
tion (n-l)  A\jj,  where  is  the  interblade  angular  spacing.  The 

rotation  is  about  the  Z axis.  Note  that  BLn  sets  are  rotating 

coordinates  and  have  a common  Z axis. 


VIRTUAL  PIVOT 
POINT  \ 


(FORWARD)  XF 


YP  (RIGHT) 


V i 


:H»  f 


(?) 


(?) 


Zp  (DOWN) 


THERE  IS  NO  YAW  (Z)  ROTATION  BETWEEN  H AND  F SETS 


Figure  5-3.  Coordinate  Systems  - Hub  (Nonrotating  Shaft  Top) 
to  Fuselage  Axis  (Flexible  Shaft) 
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5.2.5  Blade  Element  Coordinates 


^ XBLE ’ YBLE’  ZBLE^ 


The  blade  element  set  origin  is  located  at  the  center  of  gravity  of  an 
element  of  a particular  blade.  See  Figure  5-6.  Reference  to  a column 
vector  subscripted  by  BLE  is  used  to  denote  the  blade  element  located  by 
the  blade  element  set  origin.  The  right-hand  coordinate  triad  of  this  set 
has  the  X axis  parallel  to  the  local  quarter  chord  line,  the  Y axis 
along  the  chord  line  toward  the  leading  edge.  The  Z axis  is  mutually 
perpendicular  and  pointed  up.  The  BLE  set  is  used  to  track  the  local 
feather  angle,  to  develop  aerodynamic  and  dynamic  loading  terms. 

The  BLE  set  origin  for  each  blade  element  specifies  the  element  c.g.  with 
respect  to  the  quarter  chord,  and  in  terms  of  the  BLE  directions,  i.e., 
for  the  Kth  element  the  position  coordinates  are  SY(K)  and  SZ(K). 

SX(K)  is  the  blade  radial  station.  Transformations  to  the  neutral,  no- 
stretch axis  are  made  for  X deflections.  Note:  The  quarter  chord  is 
merely  a convenient  reference  datum,  and  does  not  convey  any  model 
limitations  or  assumptions. 


5.2.6  Free stream  (Earth)  Set  (Xg,  Yg,  Zg) 


The  freestream  set  is  essentially  the  earth  or  inertial  set  inasmuch  as 
the  axis  alignments  are  the  same.  However,  the  freestream  set  can  assume 
any  origin.  Thus  the  use  of  the  set  is  to  reference  the  local  gravity 
vector  and/or  an  absolute  angular  displacement  or  linear  velocity  accel- 


eration of  another  set.  As  shown  on  Figure  5-7,  the 


axis  points  down 


toward  local  gravity.  Other  sets  reference  to  the  E set,  as  the  H set 
shown  here,  may  assume  any  starting  value  of  roll  and  pitch  such  as  the 
trim  initial  conditions.  The  relative  orientation  changes  with  progressing 
time  of  flight. 

With  the  freestream  set  origin  located  coincident  with  the  fuselage  set, 
the  components  of  fuselage  set  velocity  in  E set  are  u_,  v , w . These 

Cj  h £j 

components  combine  into  a trajectory  velocity  u and  path  X^.  The 

trajectory  path  is  yawed  rignt  ^ and  pitched  up  y^  from  the  E set. 
See  Figure  5-8. 


5.2.7  Gyro  Coordinates  ( X , Y , Z ) 

u u li 

A gyro  set  is  used  for  modeling  an  internal,  isolated  control  gyro.  This 
set  is  shown  in  Figure  5-9.  The  G set  origin  is  coincident  with  the 
fuselage  (F)  set  origin  and  has  the  same  sense  of  direction  and  rotation. 
Rotations  are  measured  relative  to  the  fuselage  set. 


C/4  LIME 


Figure  5-6.  Coordinate  Systems  - Blade  Element  Set 
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-7-  Coordinate  Systems  - Freestream  (Earth)  to  Hub  Axis 
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5.2.8  Swashpla'-.e  Coordinates  (X„„,  Yori,  ZOTJ 
t)ir  bir  Sir 

The  swashplate  set  serves  one  of  two  possible  functions.  For  a conventional 
control  system,  or  a system  using  an  isolated  control  gyro  (Section  5.2.7), 
the  motion  of  the  set  describes  the  motion  of  the  swashplate.  For  systems 
using  a control  gyro  directly  in  the  feathering  angle  control  path,  this 
set  describes  the  control  gyro  motions. 

As  shown  on  Figure  5-9,  the  SP  set  origin  is  located  in  line  with  the 
Z axis  and  above  the  hub  set  a distance  Z . The  SP  set  does  not 

Ubr 

rotate  with  the  rotor  shaft . For  no  deflection  of  the  SP  set,  the 
X and  Y axes  have  the  same  alignment  as  the  X and  Y of  the  hub  set. 


5.3  DEGREES  OF  FREEDOM 


The  degrees  of  freedom  of  the  REXOR  equations  are  defined  as  the  general- 
ized coordinate  variables  of  the  set  of  equations  of  motion  to  be  devel- 
oped in  Section  6.  These  degrees  of  freedom  fully  describe  the  motion  of 
the  physical  elements  of  the  modeled  helicopter,  but  each  direction  of 
motion  of  the  helicopter  may  not  have  a degree  of  freedom  directly  asso- 
ciated with  it.  The  physical  motions  may  be  described  by  a series  of 
| modal  variables  (Section  6.4)  or  through  a set  of  transformations  and  com- 

l binations  of  the  degrees  of  freedom  as  developed  in  Sections  5.4  and  5-5. 
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The  REXOR  rotorcraft  simulation  analysis  can  be  applied  to  describe  the 
vehicle-rotor-control  system  dynamic  response  for  up  to  thirty-two  fully- 
coupled  degrees  of  freedom.  These  include  the  normal  six  rigid  body  or 
vehicle  degrees  of  freedom;  rotor  speed;  and  additional  provisions  for  up 
to  twenty-three  degrees  of  freedom  defining  rotor  blade  motion  (one  mode 
per  degree  of  freedom) , flexible  swashplate  and  rotor  shaft  or  pylon 
motion,  and  a flapping  moment  feedback  control  gyro.  The  equations  of 
motion  are  written  in  a general  form  so  that  additional  degrees  of  freedom 
can  be  added  if  desired.  The  current  thirty  degrees-of- freedom  are  listed 
in  Figure  5-10,  followed  by  a discussion  describing  them  in  detail. 

5.3.1  Vehicle  or  Rigid  Body 

The  six  rigid  body  degrees-of- freedom,  three  translations,  and  three  rota- 
tions are  defined  as  motions  of  the  hub  or  principal  reference  axis  sys- 
tem, Section  5-2.2,  relative  to  freestream  (inertial)  reference  datum. 

Translational  displacements  (X,  Y,  Z)  of  the  origin  of  the  hub  coor- 

OH 

dinate,  and  rotational  displacements  ($,  0,  i|»)  about  the  hub  axes  describe 

H 

these  degrees  of  freedom.  See  Figure  5-7-  As  mentioned  in  Section  5-2.6, 
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ITEM 


SYMBOL 


TYPE  OF  MOTION 


MAIN  ROTOR 
HUB  AXIS 


X0H>  Y0H-  Z0H 

0 A ]b 
H<  VWH 


TRANSLATION  AND  ROTATION 
WITH  RESPECT  TO 
INERTIAL  REFERENCE 


ROTOR 


SHAFT  OR  PYLON  0S,  0§ 
BENDING 


BLADES 
(n  = 2 or  4) 

SWASHPLATE 

GYRO 


A1n'  A2n'  A3n 
^PHn 


ROTATION  OF  ROTOR  SE1 
WITH  RESPECT  TO  HUB 
AXES 

DEFLECTION 

OF  FUSELAGE  REFERENCE 
WITH  RESPECT  TO  HUB  AXES 

BLADE  BENDING  MODES  AND 
FEATHER/PITCH  HORN 
BENDING  OR  TORSION  WITH 
RESPECT  TO  BLADE 
ROOT  AXES 

SWASHPLATE  AXES  MOTIONS 
WITH  RESPECT  TO  HUB 
AXES 

GYRO  (INTERNAL)  AXES 
ROTATION  WITH  RESPECT 
TO  FUSELAGE  AXES 


. 5ure  5-10.  Degrees  of  Freedom 
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the  freestream  set  may  instantaneously  assume  any  reference  point; 

therefore,  only  the  time  derivatives  of  (X,  Y,  Z)_„  and  (<f>,  6,  tp)  have 

UH  H 

significance . 

In  order  to  locate  the  direction  of  the  gravity  vector  relative  to  the 
huh,  a running  calculation  of  the  Euler  angles  8^,  must  be  made. 

hi  h hi 

Since  these  are  not  degrees  of  freedom  and  therefore  not  calculated  in  the 
equations  of  motion,  they  must  be  calculated  outside  the  dynamic  equations 
as  the  time  history  proceeds . When  the  initial  orientation  of  the  hub  is 

defined,  <J>  , 0 , and  are  known  and  their  changing  values  may  be  cal- 

hi  h .hi 

culated  by  integrating  the  hub  rotation  rates  in  the  earth  or  freestream 
axes . 

5.3.2  Rotor 

The  rotation  for  the  rotor  degree  of  freedom  ij;  is  defined  as  motion  of 

R 

the  rotor  coordinate  system  relative  to  the  hub  axis  system.  This  is 
shown  in  Figure  5-5-  This  figure  also  indicates  the  change  from  Z down 
to  Z up  axis,  which  is  equivalent  to  a 180-degree  positive  rotation 
about  the  Y axis. 


5.3.3  Shaft  or  Pylon  Bending 

Shaft  or  pylon  bending  degrees-of- freedom  are  defined  as  motions  of  the 
fuselage  coordinates  relative  to  the  hub  or  principal  axis  system.  As 

shown  in  Figure  5-3,  fuselage  translations,  Xgp  and  Ygp , are  dependent 

variables  which  are  functions  of  the  two  shaft  bending  Euler  rotations , 

<}>  and  0 . These  rotations  are  about  a virtual  hinge  of  the  shaft  or 

O Q 

pylon . 


Shaft  bending  or  pylon  bending  is  assumed  to  be  small  enough  such  that 
displacement  along  the  Z axis  is  negligible.  Thus,  when  computing  the 


translation  of  the  origin  of  the  fuselage  coordinate  system  in  the  X and 

r 


Y dire-,  tions  due  to  shaft  bending,  small  angle  approximations  for  the 

r 


sine  and  cosine  of 


<f>s  and  8g 


can  be  used. 


3Y, 


These  translations  are  thus 


defined  as  products, 


in  the  Y direction  and 


,30, 


0, 


in 


the  Xp  direction.  The  partial  derivatives  represent  the  distances  from 
the  virtual  hinge  in  each  axis  to  the  fuselage  reference. 
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5.3.1*  Blades 


Each  blade's  motion  relative  to  the  rotor  coordinate  system  is  defined  in 
terms  of  four  generalized  coordinates.  These  consist  of  three  blade  bend- 
ing modes  and  a combined  feathering,  pitch  arm  bending  mode,  or  a torsion 
mode . 


5.3.1*.l  Blade  Bending 


Blade  motion  due  to  blade  bending  is  defined  by  the  following 
generalized  modal  coordinates: 


coupled  first  inplane  bending  mode 


coupled  first  flapwise  bending  mode 


coupled  second  flapwise  bending  mode 


Ordinarily  in  a modal  analysis,  the  effects  of  centrifugal  and  struc- 
tural stiffness  are  lumped  together  into  a generalized  stiffness 
which  is  simply  the  modal  natural  frequency  squared  times  the  gen- 
eralized mass.  In  contrast  to  this,  the  REXOR  analysis  separately 
treats  the  strain  energy  or  structural  stiffness  in  each  mode  and  the 
stiffening  due  to  the  centrifugal  force  field.  This  provides  the 
capability  of  being  able  to  account  for  the  periodic  variation  of 
stiffness  in  the  modes  due  to  the  reorientation  of  the  centrifugal 
force  field  with  respect  to  the  blade  principal  axis  due  to  varia- 
tions in  blade  angle.  This  feature  can  be  important  in  the  study  of 
subharmonic  stability  where  the  periodic  variation  of  coefficients 
may  be  important,  but  it  also  permits  being  able  to  make  rather  large 
changes  in  rotor  speed  and  collective  blade  angle  without  having  to 
change  blade  modal  data. 


Mode  shapes  and  natural  frequencies  are  initially  determined  for  a 
twisted  blade  at  or  near  the  collective  blade  angle  and  rotor  speed 
to  be  analyzed.  Such  effects  as  precone,  blade  sweep,  blade  droop, 
and  blade  angle  variation  are  included  in  the  REXOR  analysis  and 
couple  the  initially  orthogonal  modes.  The  elastic  bending  contribu- 
tion due  to  the  modal  deflections  is  calculated  relative  to  the 
blade's  static  shape. 


As  previously  noted,  the  blade  modes  are  initially  defined  at  some 


reference  feathering  angle,  As  time  progresses  in  the 

niijr 

analysis,  the  blade  feather  angle  varies  about  this  reference  posi- 
tion. The  mode  shapes  are  correspondingly  transformed  to  account  for 
the  difference  between  the  instantaneous  feathering  angle  and  the 
reference  feathering  angle,  at  the  same  time  accounting  for  other 


m 
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effects  such  as  the  static  and  instantaneous  shape  of  the  blades. 
This  yields  the  modal  coefficients  (partial  derivatives)  that  relate 
blade  element  motion  to  the  blade  bending  generalized  coordinates  as 
a function  of  time. 


The  vertical  and  inplane  blade  element  variational  motions,  6Y^  and 


<$Z. 


can  be  written  as  follows 


9Y. 

6Y.  = 

in  3A 


In 


(q^,t)6A 


3Y. 


In  + 3A 


3Y . 
1 


2n 


(Vt)6A2n  + 3A 


3n 


( ,t)6A 


3n 


(5-D 


and 


3Z.  3Z.  3Z. 

5Zin  - ST  <Vt)5Al„  * if  <Vt)iA2n  + SaT  (Vt)SA3n  ^ 

In  2n  3n 


where  the  given  or  input  partial  derivatives  are  the  true  modal 
coefficients  of  the  orthogonal  modes  for  the  blade  in  an  undeformed 
shape,  with  no  static  geometry  accounted  for,  and  at  the  rotor  speed 
and  collective  angle  for  which  the  blo.de  modes  were  initially 
calculated. 


The  orthogonal  bending  modes  used  in  the  analysis  are  illustrated  in 
Figures  5-11,  5-12,  and  5-13.  Observe  that  the  root  boundary  condi- 
tions for  the  modes  may  be  cantilevered  or  articulated. 


Note  that  in  addition  to  the  normal  bending  responses, 


Y. 

l 


and  Z. 


the  spanwise  motion  of  each  blade  element  is  also  determined,  and 
blade  feathering  due  to  pitch-lag  and  pitch-flap  kinematic  coupling 
effects  are  also  accounted  for  in  each  blade  bending  mode.  This 
feathering  is  added  to  that  due  to  swashplate  motion  as  is  blade 


feathering  due  to  flexibility. 


This  modal  data  is  developed  to  the  form  used  in  the  blade  equations 
in  Section  5.5.5.  The  discussion  of  modes  is  carried  on  from  a math 
viewpoint  in  Section  6.b. 


5 . 3 . ^ . 2 Pitch  Horn  Sending  - Dynamic  Torsion 


The  remaining  mode  per  blade,  pitch  horn  bending,  is  comprised  of 
either  a blade  feathering  drive  flexibility  with  a torsionally 
rigid  blade  or  an  uncoupled  torsion  mode.  Examining  the  first 
alternative,  the  swashplate  position  determines  the  primary  blade 
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feathering  motion.  In  addition,  the  linkage  between  the  swashplate 
and  the  blade  (see  Figure  5-1*0  has  flexibility  in  the  pitch  link, 
pitch  horn,  and  cuff.  The  feathering  or  pitch  horn  bending  degree- 
of-freedom  therefore  can  be  rigid  blade  feathering  motion  outboard  of 
the  blade  cuff  coupled  with  a net  inboard  stiffness.  Inboard  of  the 
blade  cuff,  feathering  flexibility  results  from  the  pitch  link,  pitch 
link  bearings,  pitch  horn,  and  cuff.  The  relationship  between  blade 

feathering,  <f>  , and  motion  of  this  degree-of-freedom, 


defined  as  the  partial  derivative, 


PH  /n 


Alternatively,  this  degree  of  freedom, 


can  be  a distributed 


torsional  response  of  the  blade  based  upon  defining  an  uncoupled 
dynamic  torsion  mode.  The  selection  of  the  degree-of-freedom  repre- 
sentation is  made  on  the  basis  of  the  type  of  analysis  being  per- 
formed. The  mode  defined  is  uncoupled  in  the  sense  that  it  is  not  a 
function  of  the  flapping  or  lead-lag  modes. 

An  optional  quasi-steady  torsional  response  of  the  blade  may  be  used 
in  conjunction  with  pitch  horn  bending.  This  is  superimposed  on  the 
rigid  blade  feathering  and  permits  a distributed  torsional  response 
alternative  of  the  blade  reacting  the  spanwise  variation  of  applied 
torsional  moments  from  aerodynamics,  coriolis,  and  centrifugal  force 
terms.  The  blade  torsional  response  at  the  ith  blade  station  is 
computed  from  the  following  equation: 


= Pi 

TT  S+1  J TOO 


rxi  dx  rtip 

, / Vx)dx 

root  Jx± 


(5-3) 


where  S is  the  Laplace  operator,  and  x^  is  the  time  constant 

associated  with  blade  torsional  response.  This  equation  is  imple- 
mented numerically  in  the  REXOR  program. 

To  aid  in  program  trouble  shooting  the  pitch  horn  bending  representa- 
tion (with  or  without  quasi-static  torsion)  may  also  be  operated  as 
a quasi-static  degree  of  freedom  without  second-order  response. 

5.3.5  Swashplate 

The  swashplate  has  three  degrees  of  freedom:  $ , 6Sp,  and  Zgp.  Rota- 

tions <J>Sp  and  8^p  are  Euler  angles  defining  the  orientation  of  swashplate 


• ' — ,***f!«**^*W£<W 


Figure  5-1^.  Blade,  Pitch  Horn  and  Feather  Hinge  Geometry 


coordinates  relative  to  the  hub.  Likewise,  the  translation  Zgp  defines 

vertical  displacement  of  the  swashplate  relative  to  the  hub  axis.  These 
are  shown  schematically  in  Figure  5-9. 


In  a gyro-controlled  feathering  moment  feedback  system,  such  as  the 
Lockheed  concept,  the  swashplate  becomes  the  control  gyro  and  the  swash- 
plate  degrees  of  freedom  correspond  to  the  degrees  of  freedom  associated 
with  the  control  gyro.  The  selection  of  this  or  other  control  system 
configurations  is  optional  in  the  program;  fundamentally  by  data  changes 
in  terms  of  inertias  and  spring  rates . 

5.3.6  Gyro 

A gyro  is  included  in  the  present  version  of  the  analysis  which  permits 
the  description  of  a direct  flapping  moment  feedback-gyro  controlled  rotor 
system.  As  already  noted,  the  inclusion  of  this  control  system  is  a 
selectable  option  in  the  program.  The  degrees  of  freedom  associated  with 

this  gyro  are  <Pn  and  0_,  and  are  Euler  angles  defining  the  orientation 
G G 

of  the  control  gyro  relative  to  the  fuselage  axis  system.  These  degrees 
of  freedom  are  shown  schematically  also  in  Figure  5-9- 


5.4  GENERAL  MOTION  AND  COORDINATE  TRANSFORMATIONS 

In  development  of  the  equations  of  motion,  it  is  convenient  to  write  the 
forces,  moments,  velocities,  and  accelerations  „ coordinate  systems 
related  to  separate  elements  of  the  system.  Consider  the  concept  of 
general  space  motion  of  a particle. 


5.4.1  General  Case  of  Space  Motion 

For  the  general  case  of  space  motion, a particle,  p,  moves  with  respect 
to  a reference  axis  system  which  is,  in  turn,  in  motion  with  respect  to  a 
fixed  coordinate  system.  This  is  illustrated  in  Figure  5-15  where  the 
fixed  or  inertial  coordinate  system  is  designated  by  capital  letters 
X,  Y,  Z,  and  the  moving  coordinate  system  is  designated  by  lower  case 
letters  x,  y,  z.  The  moving  coordinate  system  is  rotating  at  an  angular 

, 4 •>  # 

velocity,  a).  The  vector  w may,  in  general,  vary  m magnitude  and 

direction,  both  of  which  can  be  referenced  with  respect  to  the  fixed 
X,  Y,  Z axes. 


Thus,  the  absolute  motion  of  the  particle 


referred  to  the  inertial 


X,  Y,  Z axes,  is  equal  to  the  motion  of  the  particle  relative  to  the 
moving  coordinate  axes  x,  y,  z plus  the  motion  of  the  moving- axis  system 
with  respect  to  inertial  space. 





To  visualize  the  motion  of  the  particle  p,  let  its  motion  with  respect 
to  the  moving  axis  system  be  indicated  along  a curve  s fixed  in  the 
moving  axis  system,  x,  y,  z.  An  observer  sitting  on  the  moving  axis 
system  would  therefore  see  only  the  motion  of  p along  the  curve  s . 

From  Figure  5-15  5 the  position  of  p relative  to  the  x,  y,  z axes  is 
represented  by  the  vector 


r = xi  + yj  + zk 


(5-fc) 


where  i,  j,  and  k are  unit  vectors  along  x,  y,  z,  and  therefore  must 
be  treated  as  variables  due  to  their  changing  direction.  Differentiating 

r results  in 


r = xi  + yj 


• , , •,  . di  dj  dk 
yj  + zk  + x — + y-f-  + z — — 


dt  J dt  dt 


(5-5) 


c . di  dj 

Since  — = u)  x 1,  -r? 

dt  dt 

written  as 


-»■  -f  dk  ->■  -> 

m x j and  — = u)  x k,  this  expression  can  be 


• 1 


r = xi  + yj  + zk  + id  x (xi  + yj  + zk) 


(5-6) 


r = r + w x r 


(5-7) 


In  this  equation,  the  first  term,  r represents  the  velocity  p relative 

t m *>• 
to  the  rotating  axis,  x,  y,  z.  The  second  term,  u>  x r,  is  the  velocity 

of  the  point  in  the  moving  coordinate  system  due  to  the  rotation  to.  The 
absolute  or  inertial  velocity  R of  the  point  p is  obtained  by  adding 
the  velocity  of  the  origin  RQ  of  the  moving  axis  system  to  r , or: 


where 


R = t r + u x r 
w = pi  + qj  + rk 


0)  = pi  + qj  + rk 


(5-8) 


The  inertial  accelerations  of  the  point  p can  now  be  determined  by 
simply  differentiating  this  expression  with  respect  to  time.  Performing 
this  differentiation  yields 


R:Rg  + r + ijix(i!xr  + uxr  + 2axr 


(5-9) 


where  T-he  terms  w x to  x r and  w x r represent  accelerations  of  the 

-V 

coincident  point  in  the  moving  axis  system,  r is  the  acceleration  of  p 

relative  to  the  moving  axes,  x,  y,  z,  and  2ux  r is  the  coriolis 

acceleration  which  is  directed  normal  to  the  plane  containing  the  vectors 
-*  ->• 

u>  and  the  relative  velocity  r,  as  given  by  the  right-hand  rule. 

The  vectors  expressed  in  the  preceding  equations  are  in  the  most  general 
form  for  defining  the  motion  of  a particle  moving  in  a moving  coordinate 
system.  All  special  cases  can  be  deduced  from  these  equations. 

For  convenience,  the  time  derivative  equations  can  be  expanded  in  matrix 
form.  The  inertial  or  absolute  velocity  and  accelerations  of  the  particle 
p,  written  in  expanded  matrix  form,  are  given  by: 


0 -r 


+ 1 y f + r 


(5-10) 


0 -r 


+ y + r 


0 -r 


0 -q 


0 -r  q II  x 


0 z 


1 x 


+ I r 0 -p  y +2  r 0 -P  ly 


0 z 


(5-11) 


1 0 


cos0  0 sin0  cosiJj  -sinty  0 


T -1  = 0 

a-b 


cos<j>  -sin$  0 1 0 II  sintj>  cosiJj  0|  (5-18) 


_0  sin<()  cos  L-  sin6  0 cos0  JL  o 


By  inspection,  then,  it  can  be  seen  that 


inverse  of  [T]  = transpose  of  [t] 


h-J'1  - [VJT  ■ [Va] 


(5-19) 


Carrying  out  the  indicated  matrix  multiplication  yields  the  transformation 
matrix  [T]  : 


W ■ 


(cos^cosO)  ( sin<fisin0cosij'+cos(J)siniJj ) (-sin0cos<f>cos»|;+sin<)>siniiO 


(-sini|>cos0 ) (-sinifjsintj>sin0+cost|icos!f> ) ( sint^sin0cos<j>+cosiJjsin<J> ) 


(sin0 ) 


(-cos8sin$ ) 


( COS<j>COS0  ) 


(5-20) 


Using  this  transformation,  the  inertial  velocities  and  accelerations  of  a 
point  or  particle  be  written  in  one  coordinate  system  in  terms  of  those  in 
the  other  coordinate  system  as  follows : 


= M 


(5-21) 


- - — g £. — — 


c^^.^?'>?r,'^52^^3jj5 


= vj  5 


2 b 


(5-22) 


and  inversely, 


I 

X 

= [ Va]  ‘ 2 


(5-23) 


- W 5 


(5-24) 


5.4.3  Angular  Velocities  and  Accelerations  - General 


For  the  general  case,  consider  the  coordinates  in  the  previous  section, 


and  let  (p,  q,  r)&  and  (p,  q,  r)  be  the  respective  angular  velocities 


of  and  about  the  (x,  y,  z)  and  (x,  y,  a)  axis  systems.  Also,  assume 

8.  D 


that  the  Euler  angles  are  varying  with  time  (<f>,  0,  and  <p)  and  let 


(x,  y,  z)  be  the  reference  coordinate  set  with  (x,  y,  z),  coordinate 

0.  D 


set  moving  relative  to  it.  This  is  illustrated  in  Figure  5— IT . 


From  this  figure,  the  following  can  be  written. 


Differentiating  this  expression  with  respect  to  time  results  in  angular 
accelerations  (p,  q»  in  terms  of  the  reference  coordinate  system 
angular  velocities  and  accelerations.  This  results  in  the  following: 


These  equations  represent  a general  form  for  defining  angular  velocities 
and  accelerations  of  one  axis  system  rotating  relative  to  another  axis 
system, which  in  turn  is  in  motion. 


5 . 5 RELATIVE  MOTIONS  AND  TRANSFORMATIONS  USED  IN  THE  EQUATIONS  OF  MOTION 

In  this  section  are  presented  the  inertial  linear  and  angular  velocities 
and  accelerations  of  major  components  of  the  vehicle,  including  motion  of 
the  principal  reference  system,  the  fuselage,  the  swashplate,  the  control 
gyro,  the  rotor,  and  blade  elements.  Also  included  is  the  development  of 
coordinate  transformations  that  relate  motion  in  one  axis  system  to 
another.  Motion  of  the  principal  reference  axis  system  in  relation  to 
the  earth  is  described.  Motion  of  each  component  or  reference  axis  system 
is  then  defined  in  terms  of  the  degrees  of  freedom. 

5.5-1  Hub  Motion  in  Inertial  Space 

At  each  instant  in  time  the  hub  axis  (Section  5-2.2)  is  related  to  an 
inertial  coordinate  axis  system.  Inertial  accelerations  of  the  hub 
axis  system  are  defined  by  the  vector 


(5-27) 


where  the  quantities  represent  the  total  inertial  acceleration  of  the 
generalized  coordinates  of  the  vehicle  as  defined  by  motion  of  the 
principal  coordinate  axis  system. 

Orientation  of  this  system  relative  to  the  earth  is  specified  by  Euler 
angles  <f»E , 6^,,  and  as  seen  in  Figure  5-7.  The  sequence  and  defini- 

tion of  these  angles  is  i|>  (yaw),  0 (pitch),  <f>_  (roll).  Note  that  the 
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sequence  of  rotations  is  opposite  to  that  given  "by  Figure  5— 16 . The 


angular  rates,  p^,  q^,  r^,  of  the  huh  or  principal  reference  axis  system 


with  respect  to  the  inertial  coordinate  system  can  be  written  as 


0 f + 0 cos<J)E  sin<f>E 


0 -sin<j>E  cos<f>E 


cos0E  0 -sin0E  I 0 


sin0E  0 cos0E|  ^E 
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This  equation  can  be  rewritten  to  solve  for  <j>  , 0 and  ip„  as 
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The  Euler  angles  defining  orientation  of  the  principal  reference  axis 
system  with  respect  to  the  earth  is  next  obtained  by  integrating  the  rates 
with  respect  to  time,  or 
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tE=jo  *Edt 


(5-32) 


Angular  velocities  of  the  hub,  with  respect  to  the  inertial  axes  reference 


system,  p„,  q„,  rrr,  are  defined  in  terms  of  the  degrees  of  freedom  as 
n n rl 
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Linear  velocities  of  the  hub  or  principal  axis  system  are  now  determined. 
The  first  three  quantities  of  the  hub  axis  acceleration  vector  represent 
the  linear  inertial  acceleration  of  the  hub.  For  a system  in  motion,  the 

..  I 


inertial  acceleration,  aQ  , at  the  origin  of  the  system  is  defined. 


based  on  the  vector  algebra  of  Section  5 •*+■!>  as 


•■I  •• 

a = a + io  x v 
0 0 0 


(5-35) 
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where  a 


dV, 


0 1S  dt 


, the  rate  of  change  of  velocity,  V^,  of  the  origin 


of  the  moving  coordinate  system  and  co  is  the  rotational  velocity  of  the 
moving  coordinate  system,  both  relative  to  the  earth.  Now  making  the 
definition 
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gives 
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From  this  equation,  then,  the  rate  of  change  of  velocity  of  the  moving 
coordinate  reference  system  becomes 


X0H  " qHWH  + rHVH 


Y - r u + rj  v 
OH  HH  TH 


(5-38) 


W I H I X0H  " PHVH  + qHUH 


This  set  of  accelerations  and  the  time  integral  are  transformed  into  rotor 
and  fuselage  sets  and  represent  airflow  acceleration  and  velocity  incident 
on  the  helicopter. 


A separate  set  of  hub  accelerations  is  carried  through  the  analysis,  based 
on  the  hub  set,  which  contains  the  acceleration  due  to  gravity.  Ordinar- 
ily, gravity  is  treated  as  a force  of  mg  on  the  right-hand  side  of  the 
equations.  However,  the  gravitational  term  can  be  introduced  by  defining 


(5-39) 


where  gvlJ,  gv„,  g-.,  are  the  three  components  of  the  gravity  vector  to 
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be  defined.  The  hub  acceleration  on  the  left  may  be  defined  as  being  in 
"earth-inertial"  axes. 


The  logic  behind  this  substitution  is  as  follows.  For  a rigid  body  in 
motion,  the  equilibrium  equations  can  be  written  as 


mX  = m(u  + qw  - rv)  = F 


mY  = m(v  + ru  - pw)  = F 


(5-Uo) 


mZ  = m(w  + pv  - qu)  = Fr 


where 


Fyf  + 1 gy  rm 


(5-U1) 


F^,  Fy,  and  F^  represent  the  external  forces  acting  on  the  body,  exclusive 
of  gravitational  forces. 


Subtracting  the  gravitational  vector  from  each  side  of  the  previous 
equations  yields: 


m(X  - gx)  = m(u  + qw  - rv  - gx)  = Fx 


(5-1+2) 


m(Y  - gy)  = m( v + ru  - pw  - g ) = Fy 


(5-U3) 


m(Z  - gz)  = m(w  + pv  - qu  - g,,)  = Fz 


(5-1+1+) 


which  by  inspection  gives 


X - gx  = u + qw  - rv  - gy 


(5-1+5) 


Y - gy  = v + ru  - pw  - gy 


(5-1+6) 


Z - gz  = w + pv  - qu  - qz 


(5-1+7) 


Rearranging  these  equations  yields: 


u = (X  - gx)  - qw  + rv  + gy 


(5-U8) 


V = (Y  - gy)  - ru  + pv  + gy 


(5-49) 


w = '(Z  - gz)  - pv  + qu  + g 


(5-50) 


The  first  terms  on  the  right  side  of  the  equation  are  identified  with  the 
proposed  gravitational  acceleration  definition  of  equation  5—39  - 

Making  the  substitution: 


U X0H  ~ VH  + Vh  + gXH 
V 4 X0H  “ + PHWH  + gYH 
W H . ^0H  " PHVH  + qHUH  + gZH 
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In  this  equation,  the  accelerations  X.u,  Y_u,  and  ZOTy  are  the  degree- 

of-freedom  accelerations  of  the  hub  or  principal  reference  axis  system 
used  in  the  REXOR  analysis.  These  accelerations  represent  the  inertial 
accelerations  plus  the  equivalent  accelerations  of  the  reaction  force  to 
gravity.  Thus,  gravity  is  an  equivalent  acceleration  applied  to  the 
refeience  coordinate  axis  system.  Via  coordinate  system  referencing, 
every  mass  element  on  the  vehicle  is  therefore  acted  upon  by  this  accel- 
eration. This  avoids  including  gravitational  force  as  an  external  force 
individually  applied  to  each  mass  element. 

The  gravitational  vector  at  the  hub  is  simply  the  gravity  vector  in  earth 
axis  transformed  to  the  hub  axis  system  through  the  Euler  angle  rotations 

<t>g  > and  <f»E.  Or 
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5.5-2  Fuselage  Motion  in  Inertial  Space 

As  shown  in  Figure  5-3,  the  fuselage  is  rotated  from  the  huh  position  by 
bending  of  the  rotor  shaft.  The  origin  of  the  fuselage  coordinate  axis 
system  is  displaced  from  the  hub  axis  system  by  the  distances 


Yqp  V in  hub  coordinates 


(5-57) 


YqF  l in  fuselage  coordinates 


(5-58) 


As  discussed  earlier,  and  Y are  functions  of  shaft  or  pylon  bend- 

ing  and  hence  are  not  independent  degrees  of  freedom.  The  fuselage  rota- 
tions relative  to  the  hub  are  described  by  the  Euler  angles  -0g  and  -<j> 
taken  in  that  order.  The  minus  signs  and  reversed  order  are  taken  so  the 
transform  from  fuselage  to  hub  is  <J>S  and  then  0g. 

From  Figure  5-3,  the  location  of  the  fuselage  reference  in  hub  coordinates 
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Now,  defining  the  transformation  from  hub  to  fuselage  coordinates  as 


M 


"l 

0 

0 

COS0 

O 

0 

sin0s' 

0 

cos<j>s 

-sin^ 

0 

1 

0 

0 

sin<f>g 

COS<j>g 

-sin0s 

0 

C°S0S 

(5-60) 


The  inertial  velocities  of  the  fuselage  reference  system  in  fuselage 
coordinates  can  be  written  as 
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Differentiating  this  expression  once  more  with  respect  to  time  (Sec- 
tion 5.^.1)  yields  accelerations  of  the  fuselage  reference  point  in 
fuselage  axes. 
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These  two  equations  define  the  linear  velocities  and  accelerations  of  the 
fuselage  reference  point  in  fuselage  axes  in  terms  of  the  generalized 
coordinate  velocities  and  accelerations  of  the  hub  or  principal  reference 
axis  system,  and  in  terms  of  the  shaft  or  pylon  bending  generalized  coor- 
dinate displacements,  velocities,  and  accelerations. 


Keeping  in  mind  the  order  of  rotation  of  0 and  <j>  , and  the  rotation 

b b 

formulas  developed,  the  angular  velocities  (p,  q,  r)  , of  the  fuselage 

F 

reference  system  can  be  written  as: 
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Likewise,  the  angular  accelerations  of  the  fuselage  reference  become 
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Tha  above  equations  define  the  fuselage  reference  axis  system  angular 
velocities  and  accelerations  in  terms  of  hub  or  principal  reference 
axis  system  angular  velocities  and  accelerations  and  shaft  angular 
displacements. 

5.5-3  Motion  of  Rotor  Coordinate  Axis 

The  rotor  coordinate  axis  system  is  shown  in  Figure  5-5-  Note  that  the 

rotor  coordinate  axis  system  is  rotated  180  degrees  about  the  Y axis 

relative  to  the  hub  axis  system  at  the  time  when  the  rotor  is  at  azimuth 

position  zero.  That  is,  X and  Z change  directions.  The  rotor  coordinate 

system  then  rotates  through  the  angle  \p  from  this  position. 

R 

The  sequence  of  rotation  in  going  from  hub  to  rotor  coordinates  consists 

of  first  a 180-degree  0 rotation,  followed  by  the  \Jj  rotation. 

R 

Following  the  convention  established  in  Section  5-^.2  for  Euler  angles; 
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where  the  last  matrix  represents  the  180-degree  0 rotation.  The  next 

r t 


matrix  is  the  rotor  rotation,  i}> 


dt, 
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Since  the  origins  of  the  rotor  coordinate  system  and  the  principal  refer- 
ence axis  system  are  coincident,  the  linear  velocities  and  accelerations 
of  the  oi'igin  of  the  rotor  coordinate  system  can  be  directly  written  as: 


Yo  ' = [th-r]  I * 


= [th-r]  • * 
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Noting  gravity  has  been  treated  as  an  equivalent  acceleration  in  the  hub 
generalized  coordinate  accelerations.  This  same  equivalent  acceleration 

is  included  in  (Xn,  Y , Z ) , the  rotor  coordinate  accelerations. 

0 0 On 

The  angular  velocities,  p , q , r_ , and  accelerations,  p_,  q_ » of 

n n n n n n 

the  rotor  coordinate  system  are  determined;  again  noting  the  rotation 

order.  The  rotor  coordinate  system  angular  velocities  are: 
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Likewise,  accelerations  of  the  rotor  coordinate  system  are: 
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The  above  equations  then  define  the  coordinate  transformation  from  hub  to 
rotor  coordinates;  and  rotor  axis  system  linear  and  angular  velocities  and 
accelerations  in  terms  of  velocities  and  accelerations  of  hub  and  the  rotor 

degrees  of  freedom  tD- 

a 


5 - 5 . Blade  Coordinate  Relative  to  Rotor  Coordinates 


Interim  equations  to  describe  teetering  make  use  of  ZnT3T  , the  blade 

UdL 

undersling.  These  equations  are  approximations  given  in  Section  6.7. 

It  is  assumed  tne  blade  axes  and  rotor  axes  origins  coincide  and  ZnT)T  = 0 

UdL 

in  this  section. 

Referring  to  Figure  5—5 » observe  that  each  blade  has  its  own  coordinate 
axi s system. 

Since  each  blade  has  its  own  blade  reference  system,  the  and  Y__ 

BLn  BLn 

axes  are  rotated  with  respect  to  the  X_  and  Y axes  azimuthally  by  an 

n n 

angle  defined  by  the  equation 
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where  b is  the  number  of  blades  and  n is  the  blade  number.  This 
equation  states  that  the  XBL1  and  and  the  and  the  Yr  axes 

are  coincident. 

The  transformations  between  the  rotor  coordinate  axis  system  and  the  blade 
coordinate  axis  systems  are  defined  by  the  equation 
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Note  that  these  equations  define  blade  one  as  being  straight  aft  at  time 
zero. 
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In  the  blade  reference  axes,  the  velocities  and  accelerations  of  the 
origin  of  the  blade  reference  axis  system  become: 


and 


(5-72) 
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Likewise,  the  angular  velocities  and  accelerations  of  the  blade  reference 
axis  systems  become: 


5.5.5  Blade  Element  Motion 


The  following  blade  motion  description,  due  to  the  involved  nature  of  the 
geometry,  is  rather  lengthy.  First,  in  this  development,  the  motion  of 
the  blade  with  respect  to  the  relative  blade  coordinates  is  given.  This 
motion  is  the  sum  of  static  and  modal  deflections.  Then  the  relation  to 
freestream  coordinates  is  computed.  Partial  derivatives  are  extracted 
from  the  transformations  for  use  in  the  equations  of  motion  of  the  blade 
in  Section  6.6. 


The  blade  element  motions  for  the  nth  blade  are  defined  relative  to  the 
blade  (BLn)  coordinate  reference  axes  (Figure  5-5).  The  blade  element 
relative  motions  are  functions  of  the  static  shape,  of  blade  feathering 
and  torsional  deflection,  and  of  blade  bending  of  the  coupled  inplane  and 
flapping  modes . 


The  static  shape  includes  such  items  as  blade  twist,  <j>, 


TW’ 


hub  p 'econe 


angle,  3 , blade  droop  angle  relative  to  the  precone  angle,  y,  blade 


sweep  angle,  t 


O’ 


feathering  axis  precone,  3 , the  blade  feathering 

r A 


angle,  and  the  blade  element  center  of  gravity  location. 


The  blade  motions  about  this  static  shape  include  the  effects  of  the 
three  blade  bending  modes,  A^,  A^  and  A^,  Blade  feathering,  <j>^, 
and  blade  torsional  deflection,  <j>  . 


The  blade  element  motions  are  now  defined.  The  blade  static  position  in 
the  blade  reference  axis  system  is  first  developed.  The  blade  bending  and 
feathering  deflections  are  then  introduced.  Both  deflections  and  slopes 
are  developed  and  then  these  equations  are  differentiated  with  respect  to 
time  to  obtain  the  blade  element  linear  and  angular  velocities  and 
accelerations . 


The  blade  element  linear  motions  are  developed  in  blade  (BLn)  coordinates 
and  the  blade  element  angular  velocities  and  accelerations  are  developed 
in  blade  element  (BLE)  coordinates.  The  coordinate  transformation 


matrix 


[T  1 

L BLn-BLE 


is  also  defined  to  permit  the  transformation  of  the 


BLn-BLEj 

inertial  velocities  and  accelerations  from  one  axis  system  to  the  other. 
The  development  of  the  blade  relative  motion  equations  now  starts  with 
the  description  of  the  shape  of  the  blade. 
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5 . 5 . 5 • 1 Blade  Static  Shape 

Blade  elemental  motion  is  defined  as  motion  of  the  blade  element 
reference  axis  system  which  has  its  origin  at  the  blade  element 
center  of  gravity.  The  blade  aerodynamic  reference  axis  is  selected 
as  the  1/4  chord.  Likewise,  the  geometry  and  dynamics  are  referenced 
to  the  1/4  chord,  though  any  reference  line  could  have  been  used. 
Starting  with  the  straight  untwisted  blade  with  the  blade  1/4  chord 


lying  along  the  XT 


axis  as  in  Figure  5-18,  the  blade  element  eg 


and  blade  element  coordinate  axis  system  origin  are  defined  by  the 
coincident  point  defined  by  the  vector 


XCC-  ^ 1 ^ 
YCG(i)  ’ 
ZCG^  BLn 
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in  blade  coordinates.  The  dimension  X (i)  is  the  undeformed 

LG  BLn 

spanwise  location  of  the  cg/blade  element  origin.  The  dimension 

^G^BLr  is  the  chordwise  location  of  the  cg/blade  element  axis 

system  origin  forward  of  the  blade  1/4  chord  and  Zn_(iLT  is  any 

LG  BLn 

vertical  offset  of  the  cg/blade  element  origin  with  respect  to  the 
reference  chord  plane  of  the  blade . 

Now,  introducing  blade  twist  by  rotating  about  the  X axis 

BLn 

through  the  local  blade  twist  angles,  Figure  5-1*9 » results  in: 
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The  Roman  numeral  subscript  I denotes  the  first  of  a sequence  of 
static  line  transformations . 


Figure  5-18.  Blade  Element  CO/Origin  Location  in  Blade  Coordinates 


ZBLn'Zl 


Figure  5-19-  Effect  of  Blade  Tvist  on  Location  of  Blade 
Element  CG/Axis  System  Origin 
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At  this  point  the  subscripting,  BLE  will  be  dropped  to  simplify  the 
development.  Rewriting  the  above  equation,  we  have: 


(5-78) 
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Introducing  blade  coning,  f?  , results  in  the  location  of  the  blade 
as  shown  in  Figure  5-20.  This  results  in: 

fx(i)l  fcos£0  0 -sin30"|  fx(i)l 


sinf^  0 cos3q  z(i»  i 
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The  next  item  of  static  geometry  that  is  considered  is  blade  droop, 
Y,  and  then  blade  sweep,  xQ.  These  rotations  are  shown  in 
Figure  5-21.  Note  that  since  the  blade  sweep  and  droop  angles  are 
introduced  at  a distance  Xow  out  on  the  blade,  it  is  first  neces- 
sary  to  transfer  axes  to  this  location  before  making  the  rotations . 


Therefore,  the  blade  displacements  outboard  of  Station  Xg^  become: 


costq  -sinT0  0 cosy  0 siny  X(i) 

= sinTQ  cost q 0 0 1 0 < < Y(i) 


1 -siny  0 cosy 


xswoos6o 


XSWsinB0 
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At  this  same  station,  provisions  are  introduced  to  allow  for  offsets 
of  the  blade  in  both  the  vertical  and  horizontal  directions  by 

Z.  and  Y,  , respectively.  These  offsets  are  shown  in  Figure  5-22. 
jog  jog’ 

These  offsets  represent  displacement  of  the  blade  l/U  chord  with 
respect  to  the  blade  precone  line  at  blade  station  Xg^. 

Introducing  these  offsets,  then,  and  transferring  back  to  the  center 
of  rotation  through  Xg^  results  in  the  description  of  the  blade 
displacements  outboard  of  station  Xg^,  including  the  effects  of 
precone,  sweep,  droop,  and  offset  of  the  blade  from  the  precone  line. 


X ( i ) 

X ( i ) 

r i 

0 

xswcosBo 

Y(i) 

► = . 

Y(i) 

► + . 

Y. 

jog 

• + < 

0 

Z(i) 

IV 

Z(i) 

Ill 

Z. 

1 JOgJ 

__XSWSlnV 

(5-81) 


Figure  5-22.  Introduction  of  Blade  l/U  Chord  Offset,  Y and  Z 
With  Respect  to  Piecone  Line 
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At  this  point,  a reminder  that  the  prior  development  represents  the 
blade  displacement  inboard  of  Station  X and  the  above  equation 

D ft 

outboard  of  Station  XOII.  Therefore,  inboard  of  Station  X„lt: 

SW  SW 


082) 


Outboard  of  Station  XOTT: 

SW 


(5-83) 


With  this  in  mind,  the  remaining  developing  of  including  the  effects 
of  feathering  axis  static  precone  and  blade  reference  feather  angle 
in  describing  the  static  blade  position  continues.  No  distinction 
will  be  made  in  the  following  developments  between  inboard  of 

Station  X„I7  and  outboard  of  Station  X„,,. 

SW  SW 

Figure  5-23  shows  how  blade  feathering  is  introduced.  The  axis 
system  is  translated  to  a point  p which  is  located  at  the  inter- 
section of  the  precone  line  and  the  feathering  axis.  The  location  of 

this  point  is  a distance  2,^  along  the  cone  line,  as  shown  in  this 

figure.  The  blade  is  first  rotated  to  the  feather  axis;  then  rotated 

about  the  reference  feathering  angle,  $ , the  feathering  angle 

nilr 

for  which  the  blade  modes  are  defined.  Doing  this  results  in: 
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Figure  5-23.  Point  p and  Feathering  Axis  Precone  3 
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This  equation  defines  the  location  of  the  static  shape  of  the  blade 
in  an  axis  system  with  the  y-axis  horizontal  and  the  x-axis  aligned 
with  the  blade  static  feathering  axis.  Transforming  now  back  through 
the  feathering  axis  precone  angle  and  translating  back  to  the  rotor 
shaft  centerline  results  in  the  static  shape  of  the  blade  defined  in 
blade  coordinates , or 
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where  subscript  S refers  to  blade  static  or  undeformed  shape. 
Combining  equations  developed  so  far  results,  then, in  the  following 
two  equations  which  represent  the  static  shape  of  the  blade  for  both 
inboard  and  outboard  of  blade  station  X,. 
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Inboard  of  Station  X 
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Outboard,  of  Station  Xr 
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These  two  equations  then  define  completely  the  static  shape  of  the 
blade.  The  development  will  now  proceed  to  include  the  blade  bending 
or  elastic  deformation.  However,  before  proceeding  with  this,  the 
static  location  of  the  blade  feathering  bearing  is  defined  since 
these  will  be  used  in  the  development  that  follows. 


Referring  to  Figure  5-24,  it  can  be  seen  that  the  static  position  of 
the  inboard  feather  bearing  location  can  be  written  as : 


‘lBCOse0 


(5-88) 


*IBSin60  - 


(‘p  - lIB)(tan(eFA  - 60))COS60 


The  static  location  of  the  outboard  feather  bearing  is: 


‘obcos6o 
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W,inB0  + (l0B  - lp)(tanUFA  - 60))OOS60 


With  these  definitions,  the  analysis  will  now  proceed  to  include  the 
effects  of  blade  bending,  blade  feathering,  and  torsional  deflection. 
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Figure  5-24.  Static  Feather  Bearing  Geometry 
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5. 5. 5.°  Blade  Shape  - Elastic  Deformation 

In  the  foregoing  development,  the  analysis  has  proceeded  in  a 
completely  rigorous  fashion.  At  this  point,  though,  a departure 
from  a completely  rigorous  simulation  of  the  elemental  blade  motions 
will  be  made.  It  will  be  assumed,  as  far  as  blade  elastic  deforma- 
tion is  concerned,  that  the  cosine  of  angles,  like  precone  less  droop, 
blade  sweep,  elastic  flapping,  and  elastic  inplane  slopes,  but  not 
blade  feathering  is  approximately  equal  to  1,  and  therefore,  the 
blade  elastic  deflections,  y and  z,  in  blade  coordinates,  will  be 
assumed  to  be  equal  to  those  in  the  static  blade  element  coordinates. 
This  assumption  is  a reasonably  valid  assumption  and  is  completely 
consistent  with  standard  practice  in  the  mathematical  representation 
of  blade  element  motions . 

Additionally,  as  far  as  the  effect  on  structural  axis  reorientation 
due  to  blade  <j>  rotation,  the  effect  due  to  blade  elastic  twist  is 
considered  to  be  small  compared  to  that  due  to  blade  cyclic  and 
collective  feathering.  Also  it  will  be  assumed  that  the  contribu- 
tions to  blade  Y and  Z motion  are  small  due  to  blade  torsional 
motion,  other  than  that  due  to  local  center  of  gravity  offset. 

With  these  assumptions  in  mind,  blade  elastic  bending  will  now  be 
introduced.  The  contribution  to  elastic  blade  bending  is  simply 
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Note  that  X or  spanwise  motions  are  not  included  in  this  equation. 
Blade  spanwise  motion  will  be  determined  separately  by  utilizing 
blade  slope  data  to  determine  the  change  in  the  projected  blade 
length  upon  the  blade  X axis.  With  this  in  mind,  the  total  Y and  Z 
blade  motions  including  blade  bending,  but  not  yet  including  blade 
feathering  or  blade  elastic  twist,  is  strictly  the  sum  of  the  pre- 
vious static  line  expressions  and  the  modal  deflection.  Blade 
torsional  deflection  is  treated  as  an  independent  degree  of  freedom, 
and  therefore  is  not  included  as  part  of  these  blade  modes.  Combining 
the  previous  static  deflection  with  the  modal  deflections  gives: 


5 - 5 - 5 • 3 Blade  Feathering 


Blade  feathering  is  relative  to  the  reference  feathering  angle  <L. 


The  feather  angle,  then,  as  far  as  blade  motion  is  concerned,  is  due 


to  the  difference  in  the  total  feather  angle  ^ and  the  reference 


feather  angle  (L™  • 
nJcir 


The  blade  feathering  motion  is  introduced  similarly  to  the  way  the 
blade  reference  feathering  angle  was  introduced,  except  that  the 
feather  axis  slopes  are  due  to  the  static  position  as  well  as  due  to 
elastic  deformation  in  both  the  flapwise  and  inplane  deflection. 


If  we  let  Z*  . and  Y1 

FA  FA 


represent  the  instantaneous  vertical  and 


inplane  slopes  of  the  feathering  axis,  then  transferring  to  the 


inboard  feathering  bearing,  making  the  rotations  through  Z'  and 

r A 


Y'  to  the  feathering  axis,  rotating  through  the  delta  feather 
r A 


angle  -( <f,F-<l>REF)  cr  rotating  back  through  -Y'  and  -Z'  , 


and  then  transfer.  ;.ng  back  to  the  BLn  axis  system  results  in  the 
definition  of  the  displacements  in  blade  axis  coordinates. 


However,  before  proceeding  with  this,  the  feathering  axis  slopes 


Y'fa  and  Z'  are  defined.  The  slopes  are  simply  defined  as  the 


difference  in  the  total  static  and  elastic  deflection  of  the  outboard 
and  inboard  feather  bearings  divided  by  the  spanwise  distance 
between  the  bearings.  Then  from  Figure  5-2h  and  the  bearing  static 
location  equation: 
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where  in  terras  of  the  static  and  modal  deflections 
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In  the  development  that  follows,  the  time  derivatives  of  Y’  and 

FA 


Z'  are  required,  so  therefore,  they  are  now  defined.  Taking  the 


first  and  second  time  derivatives  of  the  slope  equations  yields: 
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Transferring  the  blade  displacements  as  indicated  above  to  the 
inboard  feather  bearing,  transforming  to  the  feathering  axis,  and 


performing  the  feathering  rotation  as  discussed  earlier,  yields  the 
following  equation  which  defines  the  blade  displacements  in  blade 
axis  coordinates : 
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This  equation  then  gives  the  blade  displacement  in  blade  coordinates, 
including  the  effects  of  the  static  shape,  blade  bending,  and  blade 
static  twist.  The  effect  of  blade  elastic  twist  is  now  considered. 


5*5*5.^  Blade  Elastic  Twist 


Blade  motion  due  to  blade  elastic  twist  is  accounted  for  by  going 
back  to  the  static  twist  equation.  Blade  elastic  twist,  <J>  is 
assumed  to  be  directly  superpositionable  with  blade  static  or  blade 
pretwist,  <J>  except  that  the  static  pretwist  takes  place  about 
the  l/U  chord, and  the  blade  elastic  twist  takes  place  about  the 


Figure 


a)  BLADE  PRETWIST,  <p (i)^  ABOUT  BLADE  REFERENCE  AXIS 


b)  BLADE  ELASTIC  TWIST,  0{i)T  ABOUT  BLADE  SHEAR  CENTER 

-25.  Blade  Static  Pretwist,  ^ and  Elastic  Twist,  <)> 
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Note  that  for  convenience  of  using  the  condensed  matrix  notation 

:tors  for  such  terms  as  £ , 

P 

have  been  used.  As  can  be  seen  in  this  equation. 


discussed  above,  the  most  general  vectors  for  such  terms  as  £ , 


XSW’  Zjog’  and  YSC 


these  have  all  been  treated  as  full  vectors.  Making  the  appropriate 
substitutions  of  course  will  result  in  the  expressions  previously 
obtained. 


It  .is  noted  that  the  equation  is  written  for  the  relative  displace- 


ment of  points  on  the  blade  outboard  of  Station  X ...  Inboard  of 

oW 


that  station,  the  displacements  are  determined  from  the  previous 
inboard  equation  or  simply  by  zeroing  out  such  terms  as  rj0g}  and 


r^j  and  substituting  unit  diagonal  transformations  for  £t^  j 


and 


Y 
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T in  the  full  equation.  Following  either  approach  yields  the 


blade  displacement  equation  for  points  inboard  of  Station  Xc. or 
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The  ith.  station  blade  displacements,  Y and  Z,  in  blade  coordinates 
for  points  on  the  blade  both  outboard  and  inboard  of  station  X, 
are  then  defined. 


SW 


5. 5- 5. 6 Blade  Element  Y and  Z Relative  Velocities  and  Accelerations 


The  blade  element  coordinate  axis  system  linear  Y and  Z velocities 
relative  to  the  blade  reference  axis  system  can  be  found  by  differ- 
entiating the  position  equation  with  respect  to  time.  Note  that  no 
distinction  will  be  made  at  this  ooint  between  outboard  or  inboard  of 


station  XgW,  but  using  t.ie  equation  for  displacements  outboard  of 
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Jei^!aJn°?hand  as,discufed  earlier,  zeroing  out  certain  terns, 
results  in  the  equations  for  velocities  or  accelerations  of  points 
inboard  of  that  station.  points 
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Note  in  the  above  equation  that  the  T matrices  are  not  time 

r 1 L y 

derivatives  of  the  IT  I matrices  but  are  derivatives  of  the  trans- 


formation matrices  with  respejt  to  the  transformation  angle  £. 
This  is  arrived  at  by  making  the  substitution  that: 
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and 
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Taking  the  time  derivative  again  of  Equation  5-111  yields  the  blade 
element  Y and  Z linear  accelerations  relative  to  the  blade  reference 
axis  system. 
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These  equations  define  the  blade  element  relative  displacement 
velocities  and  accelerations,  respectively,  required  by  the  blade 
inertial  velocity  equations  developed  shortly.  Note  that  in  the 
preceding  development  these  equations  are  written  for  the  nth  blade, 


and  with  the  exception  of  the  [g-],  [t„  J . 

ho]’  W’  (rscl’  W’  hoe)’  h}* ana  h0} matrices’ th 


matrices,  the 


terms  are  all  blade  dependent.  Remember,  also,  that  inboard  of 
the  T and  T matrices  are  unit  diagonal. 
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5 . 5 • 5 • 7 Blade  Element  Slopes 


The  blade  element  Y1  and  Z'  slopes  are  determined  by  differentiating 
the  deflection  equation  with  respect  to  the  nth  blade  radial  distance, 

These  formulations  are  used  for  quasi-static  torsion  formula- 
tion and  output.  Performing  the  required  differentiation  for  points 
along  the  blade  reference  line: 
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The  transformation  matrix. 


5 . 5 • 5 • 8 Transformation  from  Blade -to-»Blade  Element  Coordinates 

In  this  section,  the  transformation  matrix  from  blade  (root)  to  ith 
blade  element  motions  will  be  developed.  Each  blade  will  have  its  own 
transformation  matrix  for  each  ith  blade  station.  The  transformation 
matrix  will  initially  be  developed  as  the  transform  from  blade 

element  to  blade  coordinates  , {^BLE-BLnJ  ' 

The  transformation  matrix,  [^TgLE  BLnJ  > can  developed  by  referring 

to  the  development  of  the  deflection  equations.  The  first  rotation 
from  blade  element  to  blade  coordinates  is  through  the  combined  twist 

angle,  'the  second  rotation  is  through  the  negative  of  the 

precone,  B^;  the  thi’-’d  through  the  negative  of  the  sweep  and  droop 

angles,  t_,  y;  the  fourth  through  the  feathering  axis  angle,  B„.  ; 

the  fifth  through  the  negative  of  the  reference  feathering  angles 

and  the  sixth  back  through  the  negative  of  the  feathering  axis 

Khr 

precone  angle,  B . 

These  rotations  then  define  the  transformation  from  blade  element  to 
blade  coordinates,  including  the  effects  of  the  static  shape  of  the 
blade,  pretwist,  precone,  sweep,  droop,  etc.  Also  included  is  the 
effect  of  blade  elastic  twist.  Again  note  that  for  stations  inboard 

of  Station  Xgw,  “the  sweep  and  droop  angles,  and  y,  respec- 
tively, must  be  set  to  zero  in  the  formulation  of  the  transformation 
matrix  as  in  the  definition  of  the  blade  displacements  and  blade 
slopes.  This  portion  of  the  transformation  matrix  which  includes 
the  static  blade  share  and  combined  twist  is  defined  as  follows: 
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The  next  two  rotations  from  blade  element  to  blade  coordinates  are 

due  to  the  elastic  blade  bending  slopes.  Since  Y ' and  Z ' 

BEND  BEND 

are  motions  of  the  blade  elements  with  respect  to  the  blade,  then  to 
transform  from  blade  element  to  blade  coordinates  requires  negative 


rotations  of 


y'bend  and  z'bend 


to  be  included.  Finally,  the  blade 


feathering  rotation  from  the  reference  feather  angle  must  be  included. 
The  final  transformation  then,  from  blade  element  to  blade  coordi- 
nates, is  defined  by  the  following  equation: 
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where 
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and  again  where 
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inboard  of  Station  X^TT. 

bW 


Also: 


r 

Y' 

BEND 

■ 

► s 

Z' 

BEND 

l_ 

Y V 


Y 


V V 


V 


In 


2n 


3n  J 


(5-122) 


The  inverse  or  transpose  of  this  equation  yields  the  transformation 
from  blade  to  blade  element  coordinates,  or: 
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again  where 


W ■ [■] 


(5-124) 


inboard  of  station  XOIT. 

o w 


5. 5. 5- 9 Blade  Element  Angular  Velocities  and  Accelerations 


From  the  foregoing  discussion,  the  blade  element  angular  velocity 
vector  can  be  determined.  Starting  with  the  angular  velocities 
(p,  q,  r)  of  the  blade  reference  axis  system  and  systematically 

and  progressively  transforming  these  velocities  through  each  axis 
rotation  and  adding  the  respective  angular  velocity  associated  with 
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each  of  the  indicated  angular  rotations,  results  in  the  following 
equation  for  the  blade  element  angular  velocities. 
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Note  that  in  this  equation,  starting  on  the  right-hand  side  with  the 
quantities  in  the  innermost  brackets , the  blade  reference  system 
angular  velocities  are  first  transformed  through  the  increment  of 


feathering  axis  flapping  slope  due  to  bending,  Z'  , and  then  the 

FA 


feathering  axis  flapping  angular  velocity,  -Z'  , is  added.  Minus 

r A 


is  used  since  Z'  is  a negative  8 rotation.  Next,  the  resultant 


W vector  is  transformed  through  Y'  and  Y'  is  added.  This  is 

r A r A 


then  transformed  through  the  delta  feathering  angle,  A<j>  , and  the 

r 
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feathering  angular  velocity,  <j>  , is  added.  This  is  then  transformed 

r 

back  through  the  increments  of  feathering  axis  slopes  due  to  blade 
bending,  giving  the  vector: 
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which  represents  the  blade  element  angular  velocities  due  to  combined 
blade  feathering  and  blade  reference  axis  system  angular  velocities. 


Next,  the  effects  of  blade  bending  at  each  blade  station  are  intro- 
duced. The  above  vector  is  first  transformed  through  the  local 

blade  element  flapwise  bending  slope,  Z'  , and  then  the  angular 

velocity,  -Z'BEND,  is  adcie<i*  This  result  is  transformed  through 

the  blade  element  inplane  bending  slope,  Y'  and  the  inplane 

angular  velocity  due  to  blade  bending,  Y'  , is  added,  resulting 

in  the  total  vector  less  the  initial  transformation  string. 

This  vector  then  represents  the  blade  element  angular  velocities  due 
to  the  combined  effects  of  the  blade  reference  axis  system  angular 
velocities  of  the  blade  feathering  angle  and  of  the  blade  angular 
velocities  due  to  blade  elastic  bending.  The  remaining  transforma- 
tions then  include  the  static  effects  of  the  blade  feathering  axis 

precone,  8 . , the  blade  reference  feathering  angle,  blade 

r A Khr 

sweep,  Tq,  blade  droop,  y,  and  blade  or  hub  precone,  Bq,  and 

the  combined  effect  of  blade  static  and  elastic  twist,  represented 

by  Finally,  the  blade  elastic  twist  angular  velocity,  , is 


added,  giving  the  total  blade  element  angular  velocities, 


a 


BLE 
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Also  note,  as  indicated  before,  the  matrix 


N M 


has  the  value 


calculateu  if  X is  greater  than  Y.^  and  has  the  value  of  unity 

if  X is  inboard  of  station  X„,,. 

bW 

At  this  point  it  has  been  assumed  that  the  contributions  of  Y' 


and  Z'pA  are  sraall  compared  to  the  other  contributions  to  -I  q 


This  assumption  is  supported  by  referring  to  the  final  form  of  the 
above  development.  First  of  all,  both  of  these  vectors  are  small 


compared  to  •fir  5 which  is  fundamentally  the  rotational  speed 
lrj  BLn 

of  the  rotor.  Also,  both  of  the  feathering  axis  flapping  and  inplane 
angular  velocities  are  first  added  and  then  transformed  through  the 
delta  feathering  angle  and  then  subtracted,  meaning  that  fundamentally 

the  pxzncipal  magnitude  or  component  contributions  due  to  and 

Z'  are  self-cancelling. 

r A 

With  the  above  assumption: 
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The  blade  element  angular  accelerations  can  now  be  determined  by 
differentiating  this  equation  with  respect  to  time.  Again,  as  in 
the  case  of  the  angular  velocities,  the  contributions  due  to  time 
derivatives  of  the  feathering  axis  flapping  and  inplane  slope  changes 
due  to  bending  are  neglected.  With  this  assumption,  the  time 
derivative  is: 
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5.5.5-10  Blade  Element  X Motions 


previous  development,  the  equations  did  not  account  for  the 
ment  displacement,  velocity,  and  acceleration  in  the  span- 
direction.  The  method  used  to  define  these  is  one  of 
.it.,:  neutral  axis  as  the  axis  of  no  stretch  and  determining  the 

projecuiv ..  of  this  axis  onto  the  X-axis  as  the  blade  be-ds.  This 
projection,  then,  is  the  spanwise  cr  X 1 cation  of  the  neutral  axis 
in  blade  coordinates.  The  rate  of  change  of  this  projection  is  the 
spanwise  relative  velocity  and  the  second  rate  of  change  is  the 
spanwise  relative  acceleration  or  the  blade  element  neut.ral  axis 
location  or  point.  The  motions  are  then  transformed  to  the  center  of 
gravity  to  obtain  the  spanwise  motion  of  the  origin  of  the  blade 
element  reference  axis. 


in  Figure  5 -26,  the  deflected  neutral  axis  is  shown  as  a function  of 
blade  radius.  The  ( i — 1 ) and  ith  station  are  shown.  It  can  be  seen 


from  this  figure  that  as  X^li-l)  approaches  X (i),  then  the 


delta  length  of  the  blade  ( i ) - S,^(i-l)j,  can  be  written  as: 


_ bHA(l  1)^  = (XNA(l)  " XNA(l-1^BLn  + ^YNA(l)  YhA(l_l)^ 

+ (ZNAU)  - ZM(i-l}) 
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Rearranging  this  equation  and  summing  from  the  blade  root  to  the  kth 
blade  station  yields: 
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where  for  i=l, 
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Likewise, 


v*>  - v;1  • V3>  = ° 

BLn  BLn  oLn 


(5-135 


S (i)  is  simply  the  blade  length  to  the  ith  station  measured  along 

the  neutral  axis  and  Y (i)  and  Z (i)  are  the  Y and  Z locations 

N BLn  ^ BLn 

of  the  neutral  axis  in  the  blade  coordinate  axis  system  for  the  nth 
blade.  These  displacements,  along  with  their  derivatives,  will  be 
defined  later.  First,  however,  by  taking  the  first  and  cond  time 
derivative  of  X equation,  the  spanwise  velocities  and  accelerations 
of  the  blade  element  neutral  axis  point  are  determined  and  are  given 
by  the  following  two  equations. 

-/y  (i)  - Y„.(i-l)\  li„Ai)  - Y„.(i-1)\ 


X„. ( i-1 ) 


Z...  i-1 


XNAU)  - ^NA  1-1 


Y„.(i-1)\  - ( Z ( i ) - Z„.(i-l) 

’BLn  1 


X...U-1) 


Z„. (i-1) 


Z...U-1) 


XNA(i)  - XNA(i"l} 


If  Y NA(i)  is  the  distance  along  the  ith_ blade  element  chord  line 

from  the  blade  element  reference  axis  origin  or  center  of  gravity  to 
the  blade  element  neutral  axis,  then  the  blade  element  neutral  axis 
motions  can  be  written  in  terms  of  the  blade  element  motions  as: 
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Referring  to  Section  5-5,  the  time  derivative  of  the  above  equation 
is : 
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and  likewise,  the  second  time  derivative  is: 
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These  three  equations,  then,  define  the  Y and  Z displacements, 
velocities,  and  accelerations  of  the  neutral  axis  point  used  in  the 
X equations  and  time  derivatives.  Also,  the  increments  of  spanwise 
motions  due  to  the  offset  between  the  center  of  gravity  and  neutral 
axis  are  defined  by  these  same  three  equations.  This  increment 
represents  the  motion  of  the  neutral  axis  relative  to  the  blade 
reference  axis  origin,  therefore,  the  span  motion  at  the  center 


of  gravity  is  determined  by  subtracting  AX  (i)  from  the  spanwise 

NABLn 


motion  of  the  neutral  axis,  or: 
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(5-11*3) 


These  equations,  then,  along  with  the  previous  expressions  for 
X and  Z,  define  the  blade  element  relative  displacement,  velocity, 
and  acceleration  vectors  required  for  the  total  inertial  vectors 
which  follow. 
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5.5.5.11  Blade  Motion  in  Absolute  Coordinates 

To  this  point  the  blade  element  motion  has  been  defined  in  terms  of 
the  blade  axis  or  relative  coordinates.  The  elements  defined  are: 


blade  element  relative  displacements 
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blade  element  relative  velocities 
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and  blade  element  relative  accelerations 
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Using  the  method  of  Section  5-l.l»  expressions  in  freestream 
(absolute)  coordinates  can  be  written  for  use  in  the  equations  of 
motion.  The  blade  element  velocity  becomes: 
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The  blade  element  accelerations  are: 
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and  matching  rotation  terms  are  defined  in  Section  5.5-4  in  terms  of 
rotor  axis  terms  which  are  in  turn  related  to  the  principal  (hub) 
reference  axis. 
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5.5.6  Swashplate  Motion 

As  shown  in  Figure  5-9,  the  swashplate  reference  axis  system  is  defined 
with  the  Z-axis  down.  The  motion  of  the  swashplate  reference  system  is 

defined  by  three  generalized  coordinate  displacements,  Z , <j>  , and  0 , 

or  or  or* 

which  move  relative  to  the  hub  or  principal  reference  axis  system. 

The  rotations  4>gp  and  Q^p  are  taken  in  the  same  order  as  shown  in 
Figure  5-17  and  therefore,  from  Section  5-^.3  the  angular  velocities  are: 
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Where  ^SP  is  the  rotational  speed  of  the  swashplate,  and 


'SP  “ ^ R 


(5-151) 


where  is  the  rotational  speed  of  the  rotor. 

As  indicated  before,  the  swashplate  axes  do  not  rotate  at  the  rotational 
speed  ifigp.  However,  the  total  angular  velocities  reflect  the  rotational 
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rate  ^gp*  Therefore,  the  total  angular  rates  of  the  swashplate 
in  swashplate  axes  are  obtained  with  ^ = 0.  This  gives 
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Nonrotating  swashplate  angular  velocities,  subscripted  SP  , are 

„ * NR 

obtained  by  deleting  <(/  above. 

Of 

The  swashplate  angular  accelerations  can  be  similarly  determined  by 

evaluating  the  general  expression  at  ij;  = ij ; = 0.  This  yields: 
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The  vertical  velocities  and  accelerations  of  the  swashplate  are  simply 
defined  as: 
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and 
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It  is  noted  in  these  equations  that  the  Z-axis  motion  is  assumed  to 
remain  parallel  to  the  hub  Z axis. 


The  swashplate  angular  displacements  are  obtained  by  integrating  the 
angular  velocities,  or: 


^p  y ^spdt 
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and 
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Likewise,  the  vertical  displacement  of  the  swashplate  relative  to  the  hub 
or  principal  reference  axis  system  is: 


= J Z_dt  + Z, 
0 


(5-158) 


Where  ZQgp  is  the  location  of  the  swashplate  reference  point  with  respect 


to  hub  or  principle  axis  reference  points  along  the  Z-axis. 


It  is  noted  that  in  a feathering  moment  feedback-gyro  control  system,  the 
swashplate  is  used  for  the  control  gyro. 
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5.5.7  Direct  Feedback  Control  System  Gyro  Motion 


As  indicated  before,  the  equations  of  motion  incorporate  a description  of 
a control  gyro  that  is  used  in  a direct  flapping  moment  feedback  control 
system  (DFCS),  which  is  an  option  in  the  program.  The  gyro  motion  is 

defined  by  the  Euler  angle  rotations  <j>  and  0 relative  to  the  fuselage 

G ( / 

reference  aris  system  as  shown  in  Figure  5-9*  The  order  of  rotation  has 
been  selected  to  be  first  <p„  and  then  0„ , so  therefore,  the  trans- 
formation  from  the  fuselage  reference  axis  system  to  the  gyro  system  is: 
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In  the  same  manner  as  the  swashplate  equations,  the  angular  velocities  are 


written  with  iJj  = ij/  = 0. 
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The  p and  q terms  of  this  equation  give  the  angular  velocities  of  the 
nonrotating  portion  on  the  gyro,  i.e.,  its  gimbal  frame. 
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The  angular  accelerations  rotating  gyro  mass  or  inertia  are  likewise 


evaluated  from  Section  5.^.3  by  setting  \ji  = \p  = 0. 
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The  accelerations,  then,  of  the  nonrotating  part,  subscripted  G , are 


simply  obtained  by  deleting  the  first  term  of  the  previous  equation. 
Or: 
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5-5-8  Blade  Feathering  Motion 


The  feathering  occurring  at  the  feather  bearings,  Figure  5-lU,  is  taken  to 
be  the  sum  of  the  motions  of  the  following  dynamic  and  kinematic  elements: 


• Swashplate  - collective  command 


« 


i 


• Blade  bending  to  feathering  couplings 


Elastic  pitch  horn  and  associated  components 


The  total  feathering  response  is: 


*Fn  = % - Ais  cos^BLn  + V ' B1S  sin(*BLn  + V 


81 'fa 

p — — a 
°1  3A  In 


3z'fa  3z'fa  3*f 

°2  3A2  A2n  U3  3A *3n  3BpH  PHn 
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Velocities  and  accelerations  are  formed  by  differentiation.  The  desired 
relations  are: 


e0  - A1S  cos(*BI,n  + V - Bis  sin(*BLn  + V 


+ A1S  sln(*BLn  + V - B1S  cos(Vn  + *R>  *R 


3Y'fa  • 32 'fa  • 32 'fa  • 3*f  • 

_ p ill  A _ p a _ P £iL  a + £_  a 
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and  for  accelerations: 


eo  - Ais  COslt,BLn  + ^R>  - Bis  sin(*BLn  + V 


* 2 A1S  si"(*BLn  + V “ B1S  cos (^BLn 
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The  commanded  cyclic  blade  angles  are: 

fsintk.,. 


■(!) 


sin^pH  cos*pH  *sp 

C0S*PH  -sin*pH  ' 6sp 


(5-166) 


where  the  angle  *p^  is  the  pitch  horn-swashplate  connection  lead  to 

feather  axis.  See  Figure  5-27.  Thir  angle  is  computed  as  a static  value. 

It  should  be  noted  that  some  hub  configurations  carry  the  pitch  horn 

toward  the  blade  trailing  edge.  These  configurations  are  entered  in 

REXOR  by  forming  the  supplement  of  | . 

PH 


^PH  = 180  “ ^PH  (degrees) 


(5-167) 


This  angle  gives  the  correct  modeling  of  the  sense  of  rotation 
reversal  with  the  trailing  pitch  horn  geometry. 

The  velocities  and  accelei ations  of  the  command  cyclic  are  obtained  by 
differentiation. 
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Figure  5-27.  Pitch  Horn  Blade  Feathering  Phase  Angle 
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The  overall  coupling  (swashpuate  to  feathering)  gear  ratio,  d/e,  is 
expressed  as  a static  terra  plus  a first-order  collective  correction. 


(!)  ■ (!)„ * (IX  e< 
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The  commanded  collective  is: 


eo  - -ZSP/e 


(5-171) 


The  swashplate  vertical  motion,  Z 


SP’ 


is  developed  in  Section  5*5.6.  The 


value  e is  the  static  effective  crank  (pitch  horn)  arm  about  the  blade 
feather  axis.  This  crank  length  is  inputted  as  a negative  number  for  a 
trailing  pitch  horn  geometry  to  give  the  proper  sense  of  collective  for 
swashplate  vertical  translation. 


Taking  time  derivatives: 


eo  "zsp/e 


(5-172) 


and: 


6o  ■ ~ZSP/e 


(5-173) 


The  blade  bending  to  feathering  coupling  factors  are 


C^,  C2,  and  C3 


for 


the  first,  second,  and  third  blade  models.  The  blade  bending  modes  are 
described  without  a torsion  component;  this  allows  freedom  in  varying  the 
blade  sweep,  droop,  jog,  or  oLher  geometric  parameters  without  new  input 
data  for  the  blade  mode  shape.  The  torsion  either  is  calculated  sepa- 
rately along  the  blade  proper  or  as  a blade  root  component  by  pitch  horn 


bending.  The  ...  C3  factors  are  intended  to  add  a feathering  component 


to  the  blade  mode  which  would  exist  even  with  no  torsion  or  feathering 
moments.  As  such,  they  are  in  effect  the  63,  a2»  etc.,  coupling  usually 
described  in  the  literature.  These  couplings  are  usually  determined  as 
a function  of  the  distance  from  the  flap  or  inplane  mechanical  or  verti- 
cal hinge  to  a pitch  horn  projection. 
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The  factor  C^  is  defined  positive  as  nose  up  feathering  (radians)  per 

radian  aft  inplane  deflection  of  the  feather  axis  slope.  The  factors 

and  are  defined  positive  as  nose  up  feathering  (radians)  per  radian  of 

up  flapping  feather  axis  slope  for  the  first  and  second  flap  modes, 
respectively.  These  factors  require  the  inplane  feathering  axis  slope 

aY’FA  ?Z'fA  3Z’fA 

— — for  mode  one  and  the  out-of-plane  slopes  — — and  -r- — for 


the  other  two  modes.  The  slopes  are  defined  as 
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(5-175) 


(5-176) 


The  mode  partial 


, etc .,  describe  the  inplane  or  outplane  component 


displacement  at  the  inboard  or  outboard  bearing  location  due  to  mode  A 


Ao  . K • 
2n  3n 


The  length  l is  the  distance  between  bearings.  These 
B 


partials  are  the  same  type  as  described  in  Section  5-5.5  where  the  blade 
motion  is  built  up  from  the  bending,  feathering,  and  twist  components. 
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6.  EQUATIONS  OF  MOTION 


INTRODUCTION 


With  the  coordinate  systems  and  transformation  between  systems  well  in 
hand,  the  development  can  proceed  to  the  equations  of  motion.  The 
development  yields  a set  of  second-order  differential  equations  with  time 
varying  coefficients.  These  equations  are  formulated  using  the  energy 
approach  in  a form  credited  to  Lagrange . The  solution  to  the  system  of 
equations  is  in  the  time  domain  by  numerical  integration.  The  result  is  a 
time  history  of  the  displacements,  velocities,  accelerations,  and  loads  of 
the  components  of  the  helicopter  modeled.  Extra  attention  is  given  the 
main  rotor  where  the  blade  geometry  is  modeled  in  detail,  and  the  program 
treats  each  blade  separately. 

Following  the  development  of  the  equation  methodology  used,  the  math 
modeling  of  the  vehicle  component  parts  is  carried  out. 


6.2  ENERGY  APPROACH  TO  DEVELOPMENT  OF  EQUATIONS  OF  MOTION 

There  are  two  basic  approaches  to  developing  the  equations  of  motion  for 
a physical  system.  These  are: 

• Vector  summation  of  forces 

• Energy  approach. 

Given  an  equal  set  cf  conditions,  limitations,  and  assumptions,  both  pro- 
cedures should  result  in  equivalent  sets  of  equations.  The  difference  is 
in  the  ease  of  arriving  at  a complete  set  of  equations.  Note  that  force  is 
a vector,  whereas  energy  is  a scalar  quantity.  Therefore,  in  dealing  in 
terms  of  energy,  less  information  regarding  direction  needs  to  be  handled. 
Also  the  systematic  nature  of  the  energy  approach  reduces  the  risk  of 
error.  As  stated  by  Lagrange  (Mecanique  Analytique,  1788),  "The  methods 
which  I present  here  do  not  require  either  constructions  or  reasonings  of 
geometrical  or  mechanical  nature,  but  only  algebraic  operations  proceeding 
after  a regular  and  uniform  plan". 


The  starting  point  of  this  development  is  Lagrange's  equation.  It  may  be 
derived  by  postulating  Newton's  second  law,  or  from  Hamilton's  principle. 
Lagrange's  equation  may  be  written  in  the  following  form: 


d /3T  \ 8T  t 9B  , 3U 
dt  V^r/  " 9qr  8<ir  8qr  1 


where 


T is  kinetic  energy 
q is  a generalized  coordinate 
B is  dissipation  function 
U is  potential  energy  function 


is  the  generalized  force,  derived  from  the  virtual  work,  6W,  and 


is  defined  by  the  equation 


96W 


9a 


(6-2) 


~r 


Equation  6-1  will  now  be  developed  into  the  form  as  applied  in  REXOR. 

This  form  bears  a close  resemblance  to  a force  balance  equation,  but  is 
derived  from  energy  considerations.  For  clarity,  the  development  is  first 
shown  for  a set  of  discrete  mass  particles,  then,  in  the  section  that 
follows,  is  extended  to  the  distributed  elementa.1  masses  of  the  REXOR 
modeling  and  to  the  iterative  solution  scheme  used. 


In  a conventional  manner  the  equation  is  formulated  in  terms  of  generalized 
coordinates.  These  coordinates  are  a function  of  time,  and  completely 
define  the  system.  They  are  generally  not  directly  identifiable  as  a 
physical  quantity. 


The  physical,  parameters  or  elemental  coordinates  are  defined  to  be  func- 
tions of  the  generalized  coordinates  and  in  turn  a function  of  time. 
Consider  a system  to  be  composed  of  particles  whose  physical  coordinates 
are  a function  of  n generalized  coordinates.  For  the  ith  particle: 


= xi('l1>  *2*  ' ' 

• ■ V 4) 

(6-3) 

= • • 

• , q^;  t) 

(6-1*) 

■ Vv  V • • 

• . v tj 

(6-5) 

Note:  a Cartesian  coordinate  set  is  selected,  and  used  in  REXOR.  However, 

the  argument  is  true  for  an  arbitrary  coordinate  set. 
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The  functional  relationship  of  the  physical  or  constrained  coordinates  and 
generalized  coordinates  yields: 


3x.  3x.  3x. 

<5x.  = t ~ 6q  + — — 6q  + ■ • ■ + — — 6q 
1 3q  U1  3q_  H2  3q  Hn 

l n 


(6-6) 


3y . 3y . 

6yi  = 3^  S + 3^  6*2  + 


+ ~ 6q 
3qn  n 


(6-7) 


6z.  = T-i  *5q  + — - 5q  + • • • + t— - 6q 
l 3q  *1  3q  ^2  3q„  ^n 


(6-8) 


The  time  dependence  is  implicit  in  the  increments  of  the  generalized 
coordinates.  The  equation  is  strictly  true  for  infinitesimal  increments. 

In  REXOR  the  generalized  coordinates  are  distinct  from  physical  coordinates 
in  the  main  rotor  blade  descriptions.  Here  the  generalized  coordinates 
are  blade  modal  variables . The  modal  variables  represent  tangible  deflec- 
tions of  the  blade  from  a reference  position,  and  as  such  are  small  but 
not  infinitesimal  variables . 

As  the  variables  are  a function  of  time: 


3x.  3x. 
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z.  = — - q_  + T— - q_  + 
1 3q1  1 3q2  ^2 


3z. 
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(6-n) 


In  terms  of  the  ith  particle  the  kinetic  energy  for  the  system  may  be 
identified  as: 


= (*i2  + *1  + 62) 


Toward  the  particular  formulation  of  Lagrange's  equation  used  in  REXOR, 
the  first  two  terms  of  the  previously  stated  form.  Equation  6-1,  are 
developed: 


d_  /3T_\  3T_ 

dt  \3qr/  " 3qr 


(6-13) 


Performing  these  operations  for  the  ith  particle  case  and  the  rth  gen- 
eralized coordinate  and  summing  over  the  system  yields: 


/3T  \ 3T  V^/l  d 3 //■  2 -2  . 2\ 

VV  ~ 3qr  ~ frf\2  mi  dt  34r\\Xi  + Yi  + Zi  / 
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( 6-lU ) 


A useful  math  operation  of  cancellation  of  the  dots  is  developed  prior  to 
proceeding.  Recall: 


3x.  3x. 
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Then  also 


3x.  3x. 

x.  = -j—  q,  + 40  + 

l Bq^  1 9q2  2 


3x . 

* • + 7 — 4 

\ n 


(6-16) 


9x.  3x. 
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This  is  also  true  for  y and  z and  for  the  double  dot  terms  in  x,  y 
and  z. 


An  operation  to  reverse  the  order  of  spacial  and  temporal  differentiation 
is  required.  To  show  this  the  time  derivative  of  a partial  is  taken  as 


3 /3Xi 


. /3x.\  . /3x.\  . /3x.\ 
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Next  the  spacial  derivative  of  x^  is  given  as 


3x.  „ /3x.  3x. 

1 _ _3 _/ 1 + 1 • + 
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Now  since 


xi  = VV  V ' ' ' ’ 9„> 


(6-20) 


the  order  of  spacial  differentiation  is  reversible 
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and  hence 


d h i 


3x, 


dt\Dq  I 3q 


(6-22) 


Similar  relations  exist  for  y.  and  z. . 

i i 

Proceeding  on  with  the  kinetic  energy  terms: 
! N 
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Using  the  relationship  for  cancelling  dots  in  partials , reversing  the  order 
of  differentiation  and  cancelling  terms  gives 
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Then  from  Equation  6-1,  Lagrange's  Equation  in  constrained  coordinates 
with  point  masses  becomes 


v r- 9Xi + •• 3yi + ••  3zii  sb  3u 

2^mi  xi  + yi  aT  + zi  aT  + aT  + aT  = 

i=1  L r r rj  ^r  ^r 


(6-25) 


Also,  in  the  same  vein  of  defining  the  generalized  coordinates,  the 
relationship  between  the  elemental  and  generalized  forces  can  be  developed. 
This  relationship  is  developed  from  the  definition  of  virtual  work  on  a 
particle  as  the  scalar  product  of  the  applied  force  and  an  infinitesimal 
displacement.  Therefore  for  the  total  system  of  N elements, 


N i- 

= / F fix.  + F fiy.  + F 6z 
L-J  x i y.  l z.  i 

i=l  L 1 1 1 J 


(6-26) 


Using  the  definition  of  from  Equation  6-2  gives: 


3x.  3y.  3z. 

F — — + F — — + F — — 

*i  3qr  yi  3,Jr  2i  3«r 


(6-27) 


Substituting  Equation  6-27  into  Equation  6-25  yields  the  final  form  of  the 
Lagrange  energy  equation  in  constrained  coordinates  for  point  masses,  which 
is  in  the  form  from  which  the  REXOR  Equations  of  motion  are  developed. 
Making  this  substitution  and  rearranging  the  equation  yields 


t [(Vi  - V)  5 + (Vi  - ry)  % + (Vi  ' \)  ft] 


+ |?-  + |!L.  0 

3^  9tlr 


(6-28) 


The  above  equation  is  the  basis  for  the  entire  derivation  of  the  equations 
of  motion  of  REXOR.  Note  that  this  equation  is  written  for  discrete  ele- 
ment masses  and  discrete  forces.  Also,  at  any  instant  in  time  all  of  the 


ingredients  required  to  define  the  elemental  accelerations,  , y^,  z^, 


are  not  known.  Specifically,  the  generalized  coordinate  displacements  and 
velocities,  q^  and  q^ , are  known  at  any  instant  in  time  but  the  generalized 
coordinate  accelerations,  q^,  remain  to  be  determined  at  the  time  the 
elemental  accelerations  are  computed. 


The  following  section  presents  the  manner  in  which  the  foregoing  equation 
set  is  adapted  to  the  REXOR  numerical  solution  to  solve  the  equilibrium 
equations  or  equations  of  motions  for  the  generalized  coordinate  accelera- 
tions. This  development  is  first  presented  in  the  simpler  form,  for 
clarity,  for  discrete  mass  elements  and  forces  and  then  in  expanded  form 
to  include  elemental  distributed  masses  and  applied  moments. 


6.3  ITERATIVE  CONCEPT  AND  EQUATION  SET  SOLUTION  METHOD 


Given  a set  of  equations  as  developed  in  the  previous  section,  the  next 
step  is  to  establish  a method  of  solution.  The  solution  process  is 
defined  as  solving  the  equation  set  for  the  accelerations,  integrating 
the  accelerations  for  updated  velocity,  and  position;  then  substituting 
the  integrands  back  into  the  equations  to  determine  new  values  of 
accelerations. 


The  first  step  of  the  process  is  to  define  explicitly  the  accelerations 
from  the  equation  set.  In  the  process  of  implementing  the  REXOR  equations, 
it  is  desirable  to  handle  the  accelerations  as  an  estimated  plus  a correc- 
tive term.  In  generalized  coordinates  then  we  can  write 


= < 


J NEW 


+ < 


J CORR . 


(6-29) 


-'OLD 


This  operation  proceeds  on  a sequential  time  basis.  For  each  increment 
advance  in  time,  the  previous  'NEW'  becomes  the  'OLD'.  In  REXOR,  the  time 
increment  corresponds  with  a step  azimuthal  advance  of  the  main  rotor 
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blades.  However,  this  need  not  be  the  case.  The  'NEW'  accelerations  must 
be  used  in  the  numerical  integration  process  to  define  the  generalized 
coordinate  velocities  and  displacements.  But  if  some  form  of  a predictor 
on  accelerations  is  used  then  the  'OLD'  would  be  this  predicted  value  and 
in  this  case  it  would  be  an  estimated,  'EST',  value. 


Using  the  notation  'OLD'  and  'EST'  interchangeably  the  linear  elemental 
accelerations  can  be  written  at  time  t as 


I 

x. 

1 


where  the  estimated  accelerations  are  determined  using  the  generalized 


displacements  and  velocities,  q^  and  q^,  at  time  t,  and  the  generalized 


coordinate  accelerations  q^,  either  estimated  or  from  one  previous  time 
step  in  the  numerical  integration  process. 


Then,  at  any  given  instant  in  time  where  the  'EST'  elemental  accelerations 
are  thusly  determined,  it  can  easily  be  shown  that  the  corrective  ele- 
mental accelerations,  (x,  y,  z)^C0RR,  are  a function  only  of  the  gen- 
eralized coordinate  corrective  accelerations. 
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Now  making  the  substitution  of  Equations  6-27  and  6-31  into  Equation  6-25 
from  the  previous  section  and  rearranging  terms  yields  the  Lagrange 

equation  for  the  q^  coordinate  in  terms  of  the  estimated  elemental 

accelerations  and  the  corrective  generalized  coordinate  accelerations. 
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The  equations  of  motion  for  the  system  can  now  be  combined  and  presented 
in  matr.'1':-.  form. 
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where  the  matrices,  M^»  C ^ and  K will  be  defined  in  the  following 


discussion.  However  oefore  proceeding  with  this,  Equation  6-3b  is  now 
rearranged  into  the  form  actually  used  in  the  numerical  process  in  REXOR 
The  equation  is  solved  in  terms  of  the  corrective  accelerations. 


The  correction  terms  come  from  an  inversion  (or  simultaneous  equation, 
Cholesky  method)  operation  on  the  model  equation  set. 
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As  indicated  before  estimated  accelerations  in  physical  coordinates  come 
from  the  'EST'  or  'OLD'  generalized  coordinate  accelerations  and  the  cur- 
rent generalized  coordinate  velocities  and  displacements.  The  integration 
part  of  the  solution  operation  supplies  the  (q)  and  (q)  data. 


The  whole  package  operates  in  a cyclical  fashion,  as  shown  in  Figure  6-1. 
Arranging  the  solution  sequence  as  such  gives  it  some  important  attributes 
and  advantages. 

First,  to  determine  the  corrective  acceleration,  the  inverted  mass  matrix 
premultiplies  the  difference  of  applied  and  estimated  reactive  forces 
represented  by  the  quantity  in  the  large  brackets  on  the  right-hand  side 
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Figure  6-1.  Equation  Solution  Loop 
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Second,  as  will  be  shown  in  the  Section  6. it,  (blade  equations  section), 
carrying  the  running  acceleration  in  elemental  coordinates  allows  for  the 
simple  separation  cf  the  centrifugal  and  structural  stiffness  of  the  rotor 
blades  which  has  important  advantages  which  have  been  discussed.  Also, 
the  aerodynamic  loading  terms,  already  by  nature  in  physical  coordinates, 
are  easily  accounted  fo^. 

In  the  actual  application  of  Equation  6-35  to  REXOR,  distributed  elemental 
rigid  body  masses  are  associated  with  each  coordinate  point  and  applied 
moments  in  addition  to  forces  at  each  coordinate  point  are  accounted  for. 

Referring  back  to  Equation  6-27  the  generalized  force,  Q , from  virtual 

work  can  be  simply  written  in  the  following  form  to  account  for  applied 
moments  at  each  of  the  ith  grid  points  as 


°r  - 2_  F: 


Xi  3qr  + Fyi  3<1r  + \ 3<ir 


j-—  + M t-=-  + M ~ (6-37) 

\ yi  Zi  % 


The  terms  of  Equation  6-2h  in  Equation  6-28  can  be  developed  for  the  dis- 
tributed masses  by  going  back  to  the  elemental  acceleration  equation, 
Equation  5-11  of  Section  5.^.1  which  is  repeated  here,  in  a rearranged 
Torm,  for  clarity  of  this  development. 


of  Equation  6-35-  With  the  usual,  small,  integration  steps  these 
differences  will  be  relatively  small.  Therefore,  inaccuracies  in  the 
mass  matrix  or  its  inversion  process  only  slightly  affect  the  total 
acceleration  determination.  This  means  approximations  and  simplifications 
to  the  mass  matrix  are  acceptable.  In  some  instances,  a diagonal  mass 
matrix  will  give  convergence  to  the  required  solution. 


The  total  partial  derivatives  relating  the  motion  of  the  coordinate  point 
and  set  imbedded  within  each  elemental  body  and  the  motion  of  the 
generalized  coordinate  becomes 


(6-47) 
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(6-48) 


(6-49) 


where  on  the  right  side  these  equations,  x,  y,  and  z represent  the 
location  of  the  distributed  masses  within  the  rigid  body  elemental  mass, 

and  Xq,  y^  and  zQ  represent  the  motion  of  the  mass  element  reference 

point . 

For  each  jth  coordinate  of  the  system,  the  elements  of 

Equation  6-24  can  be  written  by  substitution  of  Equations  6-4o,  6-47, 

6-48  and  6-49.  This  gives 
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and  finally  collecting  and  grouping  terms  yields  the  final  and  complete 
definition  of  the  terms  of  Equation  6-35  for  the  estimated  elemental 
accelerations . 
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where  in  this  case  the  summation 


represents  summation  over  the  jth 


rigid  body  element.  With  this  in  mind,  substituting  Equations  6-52  and 
6-37  back  into  Equation  6-25  yields  the  complete  form  of  the  Lagrange 
energy  equation  in  constrained  coordinates  with  distributed  elemental 
masses  and  forces  from  which  the  REXOR  equations  of  motion  are  all  devel- 
oped. Also  including  these  terms  as  well  as  the  moment  terms  of  Equa- 
tion 6-37  in  Equation  6-35  yields  the  final  form  of  the  equation  as  used 
in  REXOR.  This  form  will  be  presented  following  the  development  of  the 
generalized  mass  damping  and  stiffness  matrices. 

From  Equation  6-33  it  is  easily  seen,  by  examining  the  coefficients  of 
the  corrective  accelerations  that  generalized  mass  matrix  elements, 

M , , can  be  written  as 
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This  equation  is  for  point  masses.  Actually,  as  discussed  earlier,  the 
REXOR  equations  model  a set  of  distributed  masses  characterized  by  an 
overall  mass,  center  of  gravity,  and  moment  of  inertia  values.  As  shown 
in  the  previous  section,  extension  to  the  distributed  mass  form  is  made 
by  describing  the  particle  absolute  coordinates  in  terms  of  the  position 
of  a relative  coordinate  set  in  inertial  space  ana  the  particle  position  in 
terms  of  this  relative  set  as  developed  in  Section  For  a rigid  body 

the  associated  relative  set  and  the  particle  associated  with  the  body  main- 
tain a fixed  relationship.  The  summing  over  the  particles  of  the  system 
then  becomes  a sum  over  products  of  masses  and  lengths  yielding  mass 
moment  end  moment  of  inertia  terms. 

The  mass  elements  can  be  developed  by  substituting  the  partial  derivatives 
developed  in  the  preceding  discussion.  These  partials  describe  both  the 
motion  of  the  mass  element  reference  and  also  the  distributed  masses  within 
the  rigid  body  elemental  masses. 


Substituting  these  partials,  Equations  6-47,  6-48  and  6-49  in  the 
generalized  mass  expression,  Equation  6-53  yields: 
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and  using  moment  of  inertia  and  mass  moment  definitions: 
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and  is  identified  as  a generalized  mass.  For  orthogonal  systems  M is 

rk 


zero  except  for  r - k.  The  development  of  REXOR  is  mostly  nonorthogonal 
coordinates,  therefore,  the  generalized  mass  matrix  has  many  off-diagonal 
terms. 


Similarly,  terms  can  be  developed  for  the  strain  (potential)  energy  and 
damping  functions. 
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The  stiffness  and  damping  matrix  terms  in  REXOR  are  defined  with  reference 
to  relative  coordinates;  which  parallels  the  physical,  configuration.  The 
coordinates  used  with  these  terms  then  should  be  on  a relative  basis. 

This  statement  at  first  appears  to  be  contradictory  to  the  premise  of  the 
equation  development.  However,  if  these  matrix  terms  were  defined  on  an 
absolute  basis  the  terms  other  than  those  associated  with  a relative 
motion  would  be  identically  zero.  The  integration  of  the  accelerations 
produces  changes  in  velocity  and  position.  These  changes  with  the  proper 
starting  reference  are  the  relative  coordinates  and  velocities. 

Equation  6-35  is  now  repeated  here  in  a slightly  expanded  form  to  include 
the  effect  of  applied  elemental  moments.  Equation  6— 37 » and  distributed 
elemental  masses,  Equation  6-52. 


- 

-1 

“ N 

■ 

V 

'/  . 3x.  9<J>. 

q. 

y 

|ra. x.  - F Y-—  + M — — — + . . . 

XL 

\ 1 1est  Xij8ql  Xi  3*i 

i=l 

- 

. 

M 

. 

rk 

• 

► = _ 

MASS 

< 

• 

MATRIX 

N 

7 ..  _ \^i  9*i 

*n 

CORR 

L 

fm.x.  - F hr — + M -r — + • • • 

_\  1 1EST  xir\  xi  9“l 

k.  - 

/ 

_i=l 

m 

(Equation  6-55) 


(Sum  of  Equations  6-37  and  6-52) 
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Even  though  Q was  defined  as  the  generalized  forces  of  the  system,  for 
r 

the  purpose  of  further  development  and  of  the  application  of  the  above 
equation  in  the  REXOR  analysis,  each  line  in  the  large  brackets  on  the 
right  side  of  Equation  6-60  will  hereafter  be  referred  to  as  a generalized 
force  or  a generalized  delta  force  and  will  be  referred  to  by  the  symbol, 
FR,  in  the  following  development. 
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6.1+  OVERVIEW  OF  ROTOR-BLADE  MODEL 


Many  elements  of  the  rotorcraft  can  be  directly  modeled  along  the  lines 
developed  in  Sections  6 through  6.2  and  systematized  in  Section  6.5  with- 
out further  ado.  However,  there  are  enough  special  considerations  and 
concepts  involved  in  modeling  the  individual  blades  and  combined  rotor  to 
justify  a separate  section  to  address  these  topics. 

6.1+.1  Concept  of  Modes 

The  basic  textbook  principles  governing  solutions  for  eigenvalues  (natural 
frequencies)  and  eigenvectors  (mode  shapes)  for  systems  of  several  degrees 
of  freedom  can  be  applied  to  those  of  many  degrees  of  freedom.  For  each 
independent  degree  of  freedom  there  is  an  additional  natural  frequency  and 
mode  shape. 


Free  vibrations  of  continuous  systems  such  as  beams,  or  for  example  the 
helicopter  fuselage,  or  rotor  blades,  are  generally  analyzed  mathematically 
by  reducing  the  system  to  a system  of  discrete  masses  and  elastic 
constraints . 


6.4.2  Blade  Bending  - Modal  Variable 

The  blade  is  a twisted  rotating  beam  and  its  analysis  requires  considering 
the  coupled  flapwise-chordwise-torsional  response  of  the  blade.  For  the 
REXOR  analysis,  coupled  flapwise-chordwise  mode  shapes  are  used,  upon 
which  is  superimposed  one  of  a number  of  torsional  response  representations 
of  varying  complexity  (Sections  5-?-4,  6.4.8,  6.6. 5,  and  6.6.6). 


If  one  applies  generalized  coordinates , each  blade  mode  in  the  analysis 
may  be  treated  as  a single  degree  of  freeuom.  The  generalized  coordinates 
are  called  normal  coordinates  for  the  special  case  when  the  modes  are 
orthogonal,  in  which  case  the  generalized  mass  matrix  reduces  to  a diagonal 
matrix,  as  does  the  generalized  stiffness  matrix. 


The  REXOR  analysis  uses  blade  modes  calculated  for  the  blade  at  a fixed 
rpm,  fixed  collective,  and  in  an  unswept,  unconed  orientation.  Since  the 
program  allows  for  variation  of  all  of  these  parameters,  which  is  accounted 
for  in  the  overall  REXOR  analysis,  the  predetermined  modes  become  non- 
orthogonal  as  used  in  the  program.  Thus,  blade  motion  is  effectively 
described  by  a set  of  modal  variables,  each  representing  a characteristic 
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frequency,  and  a set  of  modal  coefficients  that  describe  the  relative 
amplitude  of  oscillation  for  each  blade  segment  and  each  frequency. 

Since  the  modes  are  nonorthogonal,  ve  will  find  in  REXOR,  as  would  be 
expected  in  such  a case,  off-diagonal  coupling  terms  in  both  the  gen- 
eralized mass  and  stiffness  matrices.  It  can  readily  be  shown  in  cases 
where  generalized  or  normal  coordinates  are  applied,  that  relatively  few 
modes  need  be  taken  to  define  accurately  the  time-history  of  blade  deflec- 
tion. This  assumes  that  the  primary  frequencies  of  excitation  fall  within 
the  range  of  modes  considered. 

Some  caution  should  bo  applied,  however,  in  interpreting  time-histories  of 
moment  (stress)  or  shear.  These  variables  generally  represent  the  second 
and  third  spanwise  derivatives  of  the  deflection  curve.  The  higher  the 
degree  of  spanwise  derivative,  the  greater  the  number  of  modes  required  to 
define  it. 


However,  another  unique  feature  of  REXOR  is  that  it  actually  makes  spanwise 
integration  of  the  blade  element  aerodynamic  and  inertial  distributed  load- 
ing functions  to  compute  the  moments  and  shears  rather  than  using  the 
second  and  third  spanwise  derivatives  of  the  deflection  curve.  This 
approach  greatly  enhances  the  accuracy  of  the  internal  loads  for  a given 
number  of  modes. 

6.4.3  Adapting  Modal  Description  to  Variable  Geometry 

It  is  noted  that  the  selection  of  the  modal  description  is  such  that  the 
periodic  reorientation  of  the  structural  axis , due  to  blade  feathering  and 
torsional  deflection,  with  respect  to  the  centrifugal  force  field  stiff- 
ness is  accounted  for.  This  is  accomplished  by  generating  blade  modes  as 
rotating  blade  modes  for  input  into  the  program,  generating  structural 
only  generalized  stiffness  terms  for  each  mode,  and  the  couplings  between 
modes.  The  centrifugal  stiffness  effects  are  then  included  separately 
within  the  REXOR  program.  Thus  accounting  for  the  effects  of  blade 
feathering,  etc.,  as  indicated  above. 

The  blade  equations  developed  in  Section  6.6  permit  the  easy  separation  of 
structural  and  centrifugal  stiffness  type  terms.  Including  the  centrifugal 
terms  in  the  blade  element  accelerations  and  the  structural  stiffness  terms 
. . 3U  ....  . . , . . „ „ . 


in  the  expression 


permits  accurate  simulation  of  the  effect  of  rotor 


speed.  Secondly,  due  to  the  blade  element  motion  description  of 

au 

Section  5.5.5,  the  modal  structural  stiffness,  ; — , remains  independent 

&q 

r 

of  blade  feathering  angle  in  each  mode. 
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In  addition  to  the  above  terms,  a considerable  number  of  nonlinear  terms 
are  included  in  REXOR.  These  include  blade  spanwise  acceleration,  non- 
linear Coriolis  type  terms  and  feathering  moments  due  to  products  of 
flapping  loads  times  inplane  deflection  with  respect  to  the  feathering 
axis  and  inplane  loads  times  flapping  displacements  relative  to  the 
feathering  axis.  Additionally,  £ as  used  in  REXOR  is  time  dependent 
and  includes  the  effect  of  blade  pretvist  and  the  instantaneous  blade 
feathering  angle  due  to  collective  and  cyclic  including  pitch-flap  and 
pitch-lag  coupling  effects. 


The  effects  of  static  precone,  blade  sweep,  blade  droop,  jogs  in  the 
structural  axis,  feathering  axis  precone,  etc.,  are  also  included. 

6.k.h  Blade  Mode  Generation 

The  blade  modes  can  be  determined  by  any  appropriate  classical  method  of 
analysis  for  coupled  flapwis e-inplane  bending  beams.  The  only  require- 
ment is  a cantilever  (hinge  or  hingeless)  boundary  condition  for  the  modes 
and  that  the  terms  included  in  the  homogeneous  part  of  equations  28  and  29 
of  Reference  9 be  accounted  for.  These  equations  are  repeated  here  for 
convenience.  Flapwise: 

cos2£  + EI2  sin2g)w"  + (Elg  - EI^)  sing  cosg  v"J " 


- (Tv')'  - fi  mv  + mv  = 0 


(6-65) 


and  inplanewise: 


j(EI2  - EI1 ) sing  cosg  v"  + (E^  sin2g  + EIg  cos2g)v"J" 


- (Tv' ) ' - n mv  + rav  = 0 


(6-66) 





6.4. 5 Modal  Coefficients 

Several  additional  points  need  to  be  made  regarding  modes  in  order  that 
the  equation  development  be  properly  understood.  First,  the  same  modal 
coefficients  apply  to  the  first  and  second  time  derivatives  of  the 
function,  since 


= f(x) 
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Second,  the  motion  is  not  necessarily  confined  to  one  direction.  A given 
modal  frequency  may  excite  or  couple  with  motions  in  other  directions. 

For  example: 
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6.4.6  Independent  Blades 

In  REXOR  the  blade  motions  are  computed  and  tracked  individually.  One  set 
of  equations  operates  on  a blade  in  BLE  coordinates  as  explained  in 
Section  6.4.11.  The  result  for  a time  step  is  stored  in  BLn  coordinates 
for  that  blade.  The  operating  set  in  BLE  coordinates  then  performs  the 
computations  for  the  next  blade  in  turn. 
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6.4.7  Blade  Element  Aerodynamic  Forces  - Overview 

The  functions  F , F , F , and  moment  terms  from  Section  6.2  are 

X . I . Li  . 

Ill 

primarily  aerodynamic  loads  for  the  blade  equations.  These  loads  are 
derived  from  blade  inertial  velocity  (equivalent  to  air  velocity)  and 
table  lookup  aerodynamic  coefficients  as  given  in  Section  7. 

6.4.8  Blade  Torsional  Response 

6. 4. 8.1  Pitch  Horn  Bending 

Several  alternate  approaches  to  modeling  blade  feathering  dynamics 
exist  in  REXOR.  The  simpliest  is  to  assume  the  blade  is  torsionally 
rigid,  and  that  the  flexibility  is  in  the  pitch  horn. 

6. 4. 8. 2 Quasi-Static  Blade  Torsion 


The  blade  pitch  horn  bending  description  is  improved  by  the  addition 
of  a blade  twist  dependent  on  the  moment  loading.  This  quasi-static 
torsion  is  computed  by  integrating  the  blade  pitching  moment  times 
the  torsional  flexibility  from  tip  into  the  root.  (Developed  in 
Sections  5-3.4  and  6.6.5.) 

6. 4. 8. 3  Dynamic  Blade  Torsion 

A third  approach  to  blade  torsional  response  in  REXOR  is  an  uncoupled 
torsional  mode  which  operates  as  additional  blade  twist.  This 
material  is  developed  in  Section  6.6. 

6.4.9  Radial  Integration 

For  each  element  of  a rotor  blade  the  equations  of  motion  are  formed  per 
Section  6.2.9.  As  briefly  touched  on  in  Section  6.4.6  these  data  are 
formed  in  BLE  axis.  These  elements  are  summed  to  total  equations  for 
each  blade  in  BLn  coordinates  at  the  blade  root.  This  is  explained  in 
Section  6.6.  These  blade  root  summations  are  also  used  in  the  hub  axis 
(Section  6.9) . 


6.5  EQUATION  SYSTEM  DEVELOPMENT 

6.5.1  Reference  to  Base  Operation  Matrix 

The  equation  of  motion,  as  developed  in  Sections  6.2  and  6.3  and  as  pre- 
sented in  most  general  form  by  Equation  6-60,  may  be  given  in  abbreviated 
form  as 
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' = Mrk 
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The  M's  are  the  generalized  mass  matrix  elements,  the  F^'s  are  the 


generalized  forces,  and  the  q ' s the  generalized  coordinates  or  degrees 


of  freedom.  As  explained  previously,  the  F^'s  are  the  complete  set  of 


external  forces  and  internal  reactions  computed  with  estimated  values  of 


the  accelerations,  q 's,  at  the  next  time  point.  The  q 's  are 
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then  corrections  to  the  estimated  values . 


The  generalized  mass,  M , is  developed  in  Section  6.3.  The  generalized 

ric 


force  may  be  expanded  as  (using  the  point  mass  form)  : 
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The  inertia,  damping,  and  elastic  terms  are  developed  further  in  Sec- 


tion 6.3  (see  Equation  6-60).  The  friction  force  F__,  , the  aerodynamic 

fn 
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external  forces  FA  , and  the  pilot  control  forces  F£ 


are  described  as 


needed.  Note  that  the  potential  energy  and  dissipation  terms  have  been 
directly  included  in  the  force  expression.  Where  the  stiffness  and  damping 
matrices  are  simple  diagonals,  this  is  done.  In  the  case  of  the  blade 
equations  the  distinct  stiffness  and  damping  matrix  form  (Section  6.3)  is 
computed  before  combining  all  the  applied  forces.,  internal  reactions, 
stiffness  and  damping  terms  into  an  overall  force. 


6.5.2  Organization  by  Degrees  of  Freedom 


In  developing  the  equations  of  motion  there  are  three  types  of  ingredients 
needed: 


• Generalized  masses 

• Generalized  forces 

• Partial  derivatives  (used  in  both  of  the  above  items) 
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The  equation  development  can  then  proceed  with  these  ingredients  along  one 
of  two  lines  of  organization. 

• For  every  major  rotorcraft  piece  (fuselage,  rotor,  etc.),  compute 
all  the  ingredients  and  sort  according  to  degree  of  freedom  for 
equation  use. 


o For  every  degree-of- freedom  group,  sort  through  the  rotorcraft 
pieces  for  applicable  ingredients.  Sorting  is  minimal  because  of 
close  association  of  degrees  of  freedom  and  component  parts. 


The  latter  development  is  used  here, 
given  in  Figure  5-10. 


The  degrees  of  freedom  modeled  are 


The  following  subsections  will  describe  the  appropriate  partial  differen- 
tiations, the  generalized  masses,  and  the  generalized  forces  in  detail. 

Each  generalized  mass  couples  the  inertia  of  one  generalized  coordinate 
with  another  or  itself.  The  algebraic  equations  for  each  generalized  mass 
will  be  given  only  once.  If  the  reader  cannot  find  a particular  mass 
element  under  one  subsection,  he  should  look  into  the  other  subsection 
relating  to  the  coupled  generalized  coordinate. 

6-5-3  Partial  Derivatives 

The  generalized  masses  and  forces  use  partial  derivatives  which  describe 
the  variational  motion  of  each  physical  mass  element  in  rectangular  coordi- 
nates relative  to  the  motion  of  each  generalized  coordinate.  The  partial 
derivatives  required  are  determined  from  the  generalized  mass  and  force 
expressions  for  distributed  masses  of  Section  6.2.  The  partial  derivatives 
are  easily  constructed  from  the  coordinate  transformations  which  have  been 
developed . 

In  developing  the  motions  of  a physical  mass  element  relative  to  a gen- 
eralized coordinate,  a number  of  transforms  may  be  used.  These  can  be 
categorized  as  either  linear  or  Euler  axes  transforms  which  either  displace 
without  rotation  or  rotate  without  displacement.  The  overall  partial  will 
be  the  product  of  partials  associated  with  each  of  these  transforms.  The 
typical  form  of  these  partials  will  now  be  illustrated. 

To  obtain  the  partials,  the  equations  relating  the  velocities  are  obtained 
first.  Reviewing  Sections  5-^-1  and  5-^-2,  the  velocity  relations  of 
interest  are  restated.  For  linear  transforms: 
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and  for  Euler  transforms 
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The  partials  of  interest  are  conveniently  organized  into  3 by  3 matrices. 
They  are  for  the  linear  transform: 
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For  Euler  angles  defined  in  reversed  order  or  reversed  sign  the  last  matrix 
will  differ.  Note  by  inspection  that  rotary  to  linear  derivatives  such 


as 


30 


3x 


are  all  zero.  The  derivatives  can  be  strung  together  to  get  motion 


in  a third  axis  c relative  to  motion  in  axis  a.  Abbreviating  the 
matrices : 
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assuming  in  general 
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(6-89) 


(6-90) 


(6-91) 


(6-92) 


(6-93) 
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The  abbreviations  used  are 


r = {X,  Y,  Z| 


C = {*,  8,  <l»( 


6.5.*+  Generalized  Masses 


As  discussed  before,  the  helicopter  is  assumed  composed  of  a finite  number 
of  mass  elements.  They  are  the 

• fuselage 

• tail  rotor 

• propeller 


engine  rotor 


• swashplate 


• fixed  hub 


k mass  elements  on  each  of 


blades . 


The  reader  should  realize  the  mass  matrix  is  symmetric  from  the  definition 
of  equation  (6-33)  and  interchange  of  the  order  of  differentiation. 

\r  - Mrk  '6-9T> 

Only  the  elements  in  the  diagonal  and  the  upper  right  triangle  will  be 
given  in  the  following  sections . 

Each  of  these  mass  elements  must  be  summed  for  each  of  the  generalized 

mass  matrix  elements.  Each  mass  is  handled  with  the  distributed  mass  M . 

rk 

relation  of  Section  6.2.  Fortunately,  only  the  fuselage  requires  the  full 
equation.  The  center-of- gravity  terms  drop  out  if  the  mass  motion  is 
determined  at  the  center  of  gravity.  This  situation  is  true  for  the  blade 
line  which  passes  through  the  center  of  gravity  of  the  blade  section  mass 
elements.  Only  the  fuselage  and  swashplate  have  reference  axis  origin  off 
the  center  of  gravity.  Another  simplification  is  that  cross  products  of 
inertia  exist  only  for  the  fuselage.  Each  blade  mass  element  is  con- 
sidered to  be  in  the  shape  of  a rod  lying  along  the  chord  at  the  blade 
station  in  question. 
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Certain  small  terms  and  factors  are  dropped  from  the  generalized  masses. 

As  discussed  in  Section  6.3,  the  equations  of  motion  are  solved  for  small 
incremental  corrections  to  the  accelerations.  With  this  formulation  the 
masses  can  tolerate  approximations  as  contrasted  to  the  generalized  forces. 
For  example,  some  partials  are  assumed  to  be  zero  or  one  such  as  the  hub 
to  fuselage  partials,  which  in  effect,  says  the  shaft  bending  deflections 
are  small. 


6.5-5  Generalized  Forces 


The  equation  formulation.  Section  6.3,  requires  that  precision  be  used  in 
compiling  the  generalized  forces  per  Equation  6-7^,  expanded  per  Sec- 
tion 6.3,  Equation  6-60,  to  include  rigid  body  distributed  mass  elements. 
This  formulation  includes  for  each  degree  of  freedom: 


Summation  over  all  mass  elements  of  the  mass  times  inertial 
acceleration  times  partial  derivative.  (Section  6.3  expression 
for  distributed  masses.) 


External  (aerodynamic)  loadings  times  a partial  derivative. 


• Potential  energy  and  damping  terms  or  assembled  stiffness  and 
damping  terms  with  partial  derivatives  (Section  6.3). 


For  some  degrees  of  freedom  the  applicable  mass  elements  and  the  total 
integration  are  directly  written  as  final  results  which  can  be  verified 
by  inspection.  Degrees  of  freedom  that  properly  include  summation  over 
the  main  rotor  blades  involve  some  extensive  numerical  integrations  and 
complicated  coordinate  transformations . 


6.6  BLADE  BENDING  AND  TORSION  EQUATIONS 


6.6.1  Blade  Radial  Summation 


The  contribution  from  all  the  individual  blade  sections  are  summed  to  give 
the  blade  generalized  masses  and  forces.  These  are  given  for  blade  root, 
bending,  feathering,  and  torsion  motions.  The  blade  root  values  are  then 
transformed  to  the  final  degree  of  freedom  variables  by  partial  derivatives. 
The  summation  is  carried  out  over  all  elements  of  the  rotorcraft,  including 


the  independent  blades.  Due  to  the  relative  isolation  of  one  blade's 


modes  from  another,  only  the  by  *+  submass  matrices  along  the  diagonal  of 
the  1+b  by  Ub  rotor  matrix  are  filled,  where  b is  the  number  of  blades. 


6.6.2  Partial  Derivatives 


The  generalized  masses  and  forces  utilize  partials  relating  the  X,  Y,  Z, 
<j>,  0,  and  ij)  linear  and  rotary  motion  of  each  blade  element  to  the  blade 
bending,  blade  torsion,  body,  rotor,  and  swashplate  degrees  of  freedom. 
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Only  the  blade  bending,  torsion,  and  feathering  partials  are  derived  in 
this  section;  the  blade  partials  for  the  body,  rotor,  and  swashplate  are 
to  be  found  in  their  respective  sections. 


As  developed  in  Section  5.3.^,  the  blade  torsion  may  be  modeled  either  as 
a pitch  horn  bending  or  an  uncoupled  dynamic  torsion  mode.  For  the  former 


case  the  partial 


is  a blade  spanwise  constant  multiplier  to  summa- 


"trtr 

tions  which  couple  in  the  feather  angle.  In  the  latter  case,  is  a 

3BPHn 

function  of  span  and  blade  number. 

The  first  partials  to  be  considered  are  those  relating  motions  at  any 
point  i on  the  blade  to  the  rigid  body  motion  of  the  blade  root.  These 
partials  are: 
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Note  that  X_T_,  Y_Tr,  and  ZnTr,  are  expressed  in  blade  root  coordinates; 
dLE  BEE  dLE 


while  4>_  , 6 and  <J/  are  in  terms  of  blade  element  axes  aligned 

BLE  BLE  BLE 


with  the  blade  element  principal  axes. 


Next  consider  the  blade  Y and  Z bending  response  with  respect  to  the 
blade  bending  modes.  A number  of  equations  can  be  used  to  develop  the 
required  expressions.  The  velocity  equations  from  Section  5-5-5  are 
selected  for  ease  of  analysis.  Using  cancellation  of  the  dots 
(Section  6.2) : 
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Defined  Section  5- 5-8. 
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The  input  coefficient  is  a pitch  lag  coupling.  The  coefficients, 


C2  and  C , are  the  pitch  flap  couplings  for  first  and  second  flap  modes. 


The  existence  of  two  coefficients  accounts  for  the  difference  in  position 
of  the  flap  (virtual)  hinge  for  the  tvro  flap  modes.  These  coefficients  are 
developed  further  in  Section  5- 5-8. 


Even  though  these  coefficients  are  referenced  to  the  inplane  or  flapping 
slope  in  the  inplane  or  flapping  modes,  they  are  the  total  pitch-  flap- 
lag  coupling  in  each  mode.  Therefore,  they  must  account  for  the  effect  of 
both  the  inplane  and  flapping  feathering  axis  slope  changes  in  each  mode. 
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The  angular  derivatives  with  respect  to  the  blade  bending  modes  are 
constructed  in  the  velocity  form. 
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Note  that  the  angular  derivatives,  being  applied  to  local  segments, 
presented  in  BLE  axis.  Referring  to  Equation  5-125: 


also 
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are 


Also  note  that  in  the  same  context  and  argument  of  Section  5-5-5  the 

feathering  axis  slopes,  Y'  and  Z'  have  been  neglected  in  the  above 

rA  j.  a 

angular  partials. 

Derivatives  with  respect  to  the  blade  feathering  are  also  constructed 
using  cancellation  of  the  dots. 
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Similarly,  for  angular  motion: 
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The  partials  developed  with  respect  to  are  used  directly  in  swash- 
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plate  and  rotor  summations  as  well  as  some  of  the  following  mass  and  force 


terms . 


Some  terms  require  a further  compounding  derivative, 
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for 


the  case  of  the  pitch  horn  bending  torsion  option.  Taking  the  degree  of 
freedom  to  be  pitch  horn  angular  deflection  about  the  feather  axis,  the 
constant  is  approximately  1/e.  For  ^he  dynamic  torsion  option  a different 
set  of  mass  formulations  is  used  in  terms  of  BLE  axis,  obviating  the 
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need  for  the  compounded  derivative.  As  indicated  in  Section  5*5-5* 

expressions  inboard  of  X equal  expressions  outboard  of  X . with 
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Blade  X motions  must  now  be  accounted  for.  The  assumption  of  the  neutral 
axis  as  the  axis  of  no  stretch  is  discussed  in  Section  5*5-5  and  the 
derivation  of  the  X motions  shown.  The  equation  for  the  partials  in 
BLn  axis  for  a point  on  the  neutral  axis  is  taken  from  the  formulation  for 
the  X velocities: 
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The  program  data,  however,  is  at  the  blp.de  center -of -gravity  axis.  The 
transfer  is : 
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where  the  right-hand  side  elements  are  from  Section  5-5-5  and  the  previous 


development  of  this  section.  The  distance  Y(i)„„,  = Y(i)„„  - Y(i) 
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is  the  distance  the  neutral  axis  is  from  the  center-of-gravity  axis. 
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positive  forward.  The  partials 
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are  used  in  the 
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■preceding  equation  for 
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NA  is  then  the  difference  in  X 
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motions  between  the  reference  and  the  neutral  axis,  and  is  subtracted  from 
the  neutral  axis  motions : 
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to  obtain  a center-of-gravity  value. 
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The  spanwise  variation  of  X with  feathering, 
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, can  be  derived 
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in  a manner  similar  to 
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The  formulations  are : 
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and 
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The  program  assumes 
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to  be  zero  for  generalized  mass  calcula- 
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tions.  This  is  done  because  of  the  latitude  possible  in  the  generalized 
masses  and  the  second  order  nature  of  the  term.  In  contrast,  the  dei  iva- 
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tive  |— — — 1 is  retained.  The  partial  is  set  to  unity  for 
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the  generalized  mass  terms.  The  full  equations  are  used  for  all  these 
terms  in  determining  the  generalized  forces. 

A simple  partial  derivative  is  also  needed  when  defined  as 

dynamic  torsion.  Since  torsion  occurs  along  the  bent  and  twisted  blade 
line,  in  blade  element  axis  BLE,  only  the  vertical  or  normal  to  chord 
motion  of  the  shear  center  is  of  interest,  hence, 
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is  program  input  for  the  torsion  mode  shape. 
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6.6.3  Generalized  Masses 

The  blade  generalized  masses  in  conjunction  with  partial  derivatives 
couple  blade  feathering,  blade  torsion,  blade  bending,  body  motions,  and 
rotor  tilt.  All  the  blade  generalized  masses  that  are  assumed  to  exist 

9xble 

are  given  in  the  following  table.  As  mentioned  before,  — is  assumed 


94> 


Fn 


94. 


zero  and 


BLE 


94> 


is  assumed  one  in  the  program,  although  it  is  given  in 


Fn 


the  table.  The  blade  has  a rotary  inertia  about  the  center  of  gravity 

axis  I . The  blade  also  has  inertia  I about  a vertical  axis, 

XXBLE  ZZBLE 

but  its  contribution  is  considered  small  to  the  generalized  mass  and  is 
neglected.  This  contribution  is  not  neglected  in  the  generalized  forces. 

The  coupling  of  the  dynamic  blade  torsion  with  the  body  is  assumed  zero 
(but  not  for  pitch  horn  bending).  The  dynamic  torsion  is  taken  as  a high 
frequency  mode,  and  is  assumed  effectively  decoupled  from  rigid  body 
motions  of  the  hub  or  rotor. 

Table  6-1  lists  all  the  terms  coupling  rotary  motion  at  the  blade  root, 

M,  , and  similar  terms.  However,  not  all  listed  are  used, as  certain 
*BLn  ^BLn 

approximations  are  made  in  developing  the  principal  axis  generalized 
masses  which  reduce  the  number  of  blade  coupling  generalized  masses  needed. 
Since  the  mass  matrix  operates  on  the  acceleration  error  term  rather  than 
the  total  acceleration,  these  approximations  do  not  detract  from  the 
validity  of  the  results  produced. 
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TABLE  6-1.  BLADE  GENERALIZED  MASSES 
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TABLE  6-1  - Continued 
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TABLE  6-1  - Continued 


Blade  Bending  Coupling 
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B Defined  As  Dynamic  Pitch  Horn  Bending 
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6.6.U  Generalized  Forces 


The  development  herein  proceeds  by  first  deriving  the  equations  for  the 
loads  on  an  individual  blade  element.  The  blade  element  loads  are  com- 
posed of  aerodynamic  and  inertial  components  conveniently  found  in  either 
the  blade  root  axes  BLn  or  the  blade  element  axes  BLE.  The  loads  will 
be  summed  in  BLn  axes  with  the  appropriate  transformation.  The  desired 
equations  are: 
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The  aerodynamic  loads  are  in  BLE  axes  alignment  about  the  blade  reference 
datum  line  which  is  the  quarter  chord.  A transfer  through  the  distance 


Ynn  is  made  to  the  aerodynamic  moment.  To  put  the  data  on  a common 

cgble 


basis  with  dynamic  terms.  The  blade  aerodynamics  is  detailed  in 
Section  7*  Since  only  the  blade  section  pitching  moment  is  considered. 


i = = Note  the  blade  element  is  assumed  configured  as  a 
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chordwise  rod  for  inertia;  hence 
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A number  of  blade  summations  are  desired.  All  will  be  made  in  BLn  axes 
along  the  center-of-gravity  axis.  The  loads  at  the  principal  reference 
axes  and  for  rotor  tilt  make  use  of  the  blade  root  shears  and  moment. 
These  are  simply  the  stun  of  the  k total  blade  elements. 
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and  likewise  for  root  moments. 
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The  summations  illustrated  above  are  for  the  total  inertial  and  aerodynamic 
components.  In  a similar  manner,  the  blade  root  aerodynamic  loads  are 
derived.  The  blade  root  loads  are  summed  over  all  the  blade  to  give  main 
aerodynamic  loads  for  downwash  computations  (Section  7.2.2)  in  the  manner 
the  total  main  rotor  loads  are  found  in  Section  6.9. 
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Feathering  moments  are  used  by  the  swashplate  equations  of  motion.  These 
moments  are: 
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Only  the  X component  is  used  by  the  program.  The  equations  above  trans- 
fer the  summed  blade  loads  to  the  inboard  bearing,  then  transform  them  to 
feathering  axes.  Using  the  blade  root  loads  is  correct  when  one  recalls 
that  the  blade  is  defined  as  those  portions  that  are  feathered;  the  fixed 
hub  is  excluded. 


The  blade  bending  generalized  forces  are  now  presented.  They  are: 
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The  potential  energy  is  given  as: 
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are  inputs  calculated  external  to  the  program  from  a bending  beam  model. 


El  and  El  are  the  flapping  and  chord  stiffness  about  axes  aligned 
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with  the  blade  element  principal  axes.  The  chord  and  flapping  moments, 


My  and  reflect  the  contribution  of  the  bending  moment  from  the  i 

i i 


(or  j)  mode.  The  integration  goes  from  root  to  tip.  The  K's  are 
evaluated  for  whatever  normalized  modes  are  used  as  program  input. 


The  last  equation  can  be  derived  from  the  Bernoulli-Euler  law  for  bending 
beams : 
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where  r is  a radius  of  curvature.  The  strain  energy  is 
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Substituting  in  from  the  Bernoulli-Euler  law  and  noting  that  dS  = r^dG  = 
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Partial  differentiation  gives 
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Considering  the  moments  as  a linear  sum  of  components  from  each  bending 
mode , one  has : 
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and  likewise  for  M . Then,  by  substitution,  the  desired  equation  is 

Lx 

obtained. 

The  damping  factor  is  assumed  proportional  to  the  spring  rate.  This 
assumption  in  conjunction  with  a mechanical  lead-lag  damper  operating  at 
the  feather  axis  with  respect  to  the  fixed  hub  leads  to  : 
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The  coefficients  K . used  above  can  be  directly  identified  with  the 

mj 

coefficients  K ^ developed  in  Section  6.3.  Here  the  partial  derivatives 

associated  with  the  elemental  springs  are  embodied  in  the  mode  shape 
determinations.  The  lag  damper  is  restricted  to  respond  only  to  inplane 
motion. 

The  generalized  force  is  developed  for  pitch  horn  bending  and  dynamic 
torsion.  For  pitch  horn  bending, 
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where  Mpn  is  the  total  feather  moment  as  derived  in  Section  6.10.  For 
the  uncoupled  dynamic  torsion  option. 
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is  needed.  The  spring  constant  can  be  interpreted  as 
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where  the  generalized  mass  is  computed  continuously  and  w is  the 

PPH 


nautral  frequency  of  the  uncoupled  torsion  mode,  a program  input 
constant . 


6.6.5  Quasi-Static  Blade  Torsion 


To  improve  the  pitch  horn  bending  blade  feathering  representation  a quasi- 
static blade  torsion  distribution  is  introduced.  Quasi-static  torsion  is 


computed  from  the  structural  stiffness,  GJ^,  each  station  and  the 


torque  MX  at  the  shear  center.  The  torque  is  summed  from  the  tip  to 

oL 


the  blade  station  in  question  as  shown  in  Section  5* 3-1.  The  increment  of 
twist  produced  at  a blade  station  j can  be  displayed  as: 
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assuming  a first-order  lr.g  represents  the  torsional  dynamics.  The  time  con- 
stant i.’>  chosen  to  be  representative  of  the  blade  first  torsional  mode 

frequency. 

To  obtain  this  result,  the  available  computation  elements  require  some 
further  operations.  First,  REXOR  conducts  blade  integrations  from  root 
to  tip,  in  BLn  axes.  To  obtain  tip  to  root  values: 
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Note  the  summation  is  conducted  from  root  to  the  station  j in  question. 
Thus  the  j represents  a summation  whereas  the  i represents  a blade 
station. 

Second,  these  data  are  used  to  form  the  required  torque  at  the  shear  center. 
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etc . , act  along  the 


etc . , are  employed  to 


Small  angles  are  assumed.  The  moments  / M \ etc.,  act  along  the 

\ ^LE/BLn 

BLn  axes  and  hence  the  matrix  of  lengths  /X \ etc . , are  employed  to 

\ SC/BLn 

obtain  moments  at  the  shear  center  which  are  then  transformed  into  shear 
center  axes,  subscripted  SC,  parallel  to  blade  element  center- of- gravity 
axes,  subscripted  BLE. 

The  blade  deflections  and  slope  ir.  BLr  are  also  needed  for  the  above 
expressions. 
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6.6.6  Quasi-Static  Pitch  Horn  Bending 

To  facilitate  troubleshooting  numerical  instability  problems  an  optional 
quasi-static  pitch  horn  bending  degree  of  freedom  is  available.  The  com- 
putation elements  are  the  same  as  developed  in  Section  6.6.U  except  that 
the  solution  does  not  use  generalized  masses,  is  therefore  an  uncoupled 
mode,  and  is  calculated  externally  to  the  main  computation  flow.  The 
formulation  used  is: 
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The  dynamics  are  assumed  represented  by  a first-order  lag  with  x_,,  as  the 

PH 

time  constant.  The  variable  <()__„  is  used  to  distinguish  this  formula- 
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tion  from  the  usual  8™,  symbology. 
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6.7  ROTOR  TILT  EQUATIONS 

As  an  interim  measure,  the  capability  to  model  teetering  rotor  systems  has 
been  included  in  REXOR  by  the  partial  redefinition  of  the  available  blade 
degrees  of  freedom.  The  methodology  is  to  utilize  blade  positions  n=l  and 
n=3  to  represent  the  rotor  and  suppress  blade  matrix  positions  n=2  and  n=U. 

The  four  degrees  of  freedom  per  blade  are  used  as  follows: 

• tfl  First  Inplane  (no  change) 

• #2  Coning  (first  flap-collective) 

• if 3 Teetering  (first  flap-cyclic) 

• ffk  Blade  Torsion  (currently  not  used) 

The  inplane  calculations  proceed  as  usual.  Currently  the  hard  swashplate 
option  is  specified  with  the  interim  teetering  modeling.  This  means 
kinematic  relations  are  used  between  command  functions  and  blade  angles 


produced.  The  fourth  mode,  3 


therefore  is  not  used. 


Positions  2 and  3 are  programmed  to  be  blade  flapping  modes,  and  proper 
generalized  forces  and  masses  will  be  generated  assuming  the  modal  data  is 
entered.  A fixed  root  restraint  mode  shape,  coning,  is  used  for  mode  §2 , 
and  a pinned  connection  pair  mode  shape  (teetering)  is  used  for  mode  #3. 

A rigid  flapping  mode  (straight  line)  will  also  work  for  the  teetering 
mode. 

To  generate  coning  preferentially  in  mode  tt2  and  teetering  in  mode  ff 3 the 
opposite  type  is  nulled  from  each  mode.  This  is  done  as  follows: 

• Coning  Blade  ft 1 = 1/2  (Blade  ffl  + #3) 

• Coning  Blade  ft 3 = 1/2  (Blade  ffl  + #3) 

• Teetering  Blade  ffl  = 1/2  (Blade  tfl  - ff3 ) 

• Teetering  Blade  ff3  = -1/2  (Blade  ffl  - ff3 ) 

This  procedure  is  used  for  the  blade  generalized  forces  and  accelerations. 

Modifications  are  made  to  the  generalized  mass  matrix.  As  shown  in 
Figure  6-2,  the  elements  involved  with  the  unused  blade  positions  2 and  4 
are  zeroed  except  for  diagonal  ones  to  prevent  singularities.  Note  that 
the  teetering  is  associated  with  roll,  pitch,  longitudinal,  and  lateral 
hub  velocities;  whereas  coning  gives  rise  to  yaw  and  vertical  velocities. 
The  remaining  zeroes  shown  in  Figure  6-2  remove  the  teetering  coupling 
with  coning  associated  hub  motions  and  vice  versa. 


The  remaining  element  of  the  teetering  representation  is  the  huh  undersling 

(Zqbl)j  which  is  taken  positive  above  the  teeter  axis.  The  effect  is 

approximately  modeled  by  shifting  the  final  blade  X coordinate  by  the 
undersling  times  the  teeter  mode  root  angle. 
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6.8  SHAFT  BENDING  EQUATIONS 

6.8.1  Prime  Contributions  From  Fuselage 

The  shaft  equations  of  motion  involve  only  the  fuselage  and  shaft  forces 
and  masses.  This  simplification  results  from  the  hub  axes  being  the 
principal  axes.  Therefore,  only  derivatives  of  elements  physically  con- 
nected to  the  shaft  branch  of  the  system,  i.e.,  the  fuselage,  need  be 
considered. 

6.8.2  Partial  Derivatives 

The  needed  partial  derivatives  are  those  relating  fuselage  motions  to 
shaft  bending  motions.  The  reader  should  recall  that  shaft  roll  right  and 
nose  up  are  positive  for  hub  motions  in  fuselage  axes.  Here  the  interest 
is  in  fuselage  motions  in  hub  axes  where  fuselage  roll  left  and  nose  down 


w i-r: 


are  positive.  The  desired  derivatives  are  obtained  from  the  fuselage 
equations  of  Section  5 -5- 2. 
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where  the  elements  are  given  in  Section  5.5.2. 

6.8.3  Generalized  Masses 

Table  6-2  lists  the  generalized  masses.  It  includes  cross-coupling  terms 
with  the  principal  axis.  Note  that  only  upper  diagonal  terms  are  given. 
The  lower  diagonal  terms  are  obtained  by  symmetry.  Certain  terms  are 
Judged  negligible  and  are  dropped. 
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6.9  PRINCIPAL  REFERENCE  AXIS  EQUATIONS 

6.9*1  Nonzero  Contributions  From  Most  Vehicle  Mass  Elements 

The  principal  reference  axis  equations  of  motion  consider  contributions 
from  most  of  the  physical  elements  of  the  rotorcraft.  The  elements 
involved  are: 

• Main  rotor  - defined  as  all  portions  that  can  be  feathered 

• Rotor  hub  - includes  all  portions  of  the  main  rotor  assembly  that 
cannot  be  feathered,  and  is  treated  as  a rigid  body 

• Swashplate 

• Tail  rotor 

• Propeller  (if  used) 

• Fuselage 

• Engine 

The  internal  control  gyro  (if  used)  contribution  is  considered  negligible. 

The  six  rigid  degrees  of  freedom:  X,  Y,  Z,  <f>,  6,  ^ are  taken  with 

respect  to  the  stationary  hub  axes  which  are  also  the  principal  axes.  The 
origin  lies  on  the  main  rotor  shaft  where  it  is  intercepted  by  the  rotor 
reference  line. 

The  other  elements  considered  are  then  referenced  to  the  hub  axes.  The 
fuselage  is  subject  to  shaft  bending  motions  relative  to  the  principal 
axes.  The  tail  rotor,  the  propeller,  and  the  engine  are  installed  on  the 
fuselage  and  rotate  at  the  main  rotor  speed  times  the  appropriate  gear 
ratio.  Positive  rotations  are  defined  as: 


• Hub 


Swashplate 


same  as  main  rotor 


• Propeller  - Counterclockwise  looking  forward 

• Tail  rotor  - Clockwise  looking  right 

• Engine  - Counterclockwise  looking  forward 

The  propeller  and  engine  are  treated  as  rigid  rotating  bodies;  the  tail 
rotor  is  also  allowed  to  flap  (teetering  hinge,  etc.).  This  flapping  is 
considered  secondary  and  enters  only  into  the  aerodynamic  computations. 
The  main  rotor  is  allowed  a variety  of  motions;  teetering,  feathering, 
bending,  and  twisting. 


6.9-2  Partial  Derivatives 


A prime  ingredient  for  the  equations  of  motion,  Section  6.3,  as  outlined 
in  Section  6.5,  are  the  partial  derivatives.  First  considering  partial 
derivatives  concerning  the  fuselage  set  with  respect  to  the  hub  set. 
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As  usual,  partials  relating  rotary  motions  to  linear  motions  are  zero. 

For  certain  applications,  the  shaft  bending  is  small  and  some  small  terms 
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in  the  partial  equation  can  be  neglected.  Note  in  the  above,  the  axis 
notation  is  left  off  where  the  axis  notation  is  the  same  as  the  mass 
element  notation.  For  example, 
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Next  the  hub*to- blade  root  partial  derivatives  are  given. 
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Note  that  all  the  rotary-to-linear  partials  like 
deriving  these  derivatives  the  relationship 
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is  helpful  in  moving  the  rotary  quantities  from  a matrix  to  a vector. 


Derivatives  properly  ought  to  he  taken  for  the  other  physical  elements 
that  have  teen  considered  for  the  principal  axis  degrees  of  freedom. 
Omission  of  the  swashplate  partials  could  he  serious  for  swashplate  mass 
and  inertia  that  are  significant  compared  to  the  fuselage  and  rotor. 
Derivatives  for  the  other  rotating  elements  are  covered  for  the  constant 
rotational  speed  in  the  fuselage  derivatives.  The  variable  rotational 
speed  case  is  estimated  to  develop  third  or  fourth  order  chain  rule 
derivatives,  and  therefore  can  he  safely  neglected. 
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6.9*3  Generalized  Masses 

Table  6-U  lists  the  generalized  masses.  Note  that  small  terms  are  dropped 
from  the  generalized  mass  equations  but  not  from  the  generalized  forces. 
This  is  justified  because  as  discussed  in  Section  6.3,  the  equations  of 
motion  are  solved  for  small  increment  corrections  to  the  accelerations. 
With  such  a formulation,  the  generalized  masses  operate  only  on  accelera- 
tion errors,  and  can  tolerate  approximations. 

The  swashplate  and  the  hub  motions  are  given  in  the  principal  axes . The 
hub  motions  apply  to  the  hub  center  of  gravity  which  is  assumed  to  be  at 
the  origin  of  the  principal  axes.  The  swashplate,  however,  can  be  dis- 
placed vertically  from  that  origin.  The  fuselage  reference  axis  origin 
is  taken  to  be  directly  below  the  principal  axes.  Note  the  fuselage  mass 
and  center  of  gravity  include  the  propeller,  tail  rotor,  and  engine,  but 
the  moment  of  inertia  excludes  these  items. 

The  terms  can  be  checked  by  inspection  from  the  overall  generalized  mass 
formulation  of  Section  6.3  and  the  preceding  partial  derivatives. 

6.9.1+  Generalized  Forces 

As  pointed  out  in  Section  6.5,  the  generalized  forces  contain  both  the 

external  forcing  terms  (F  , etc.)  and  the  internal  reaction  terms 

i 

such  as  nmx^.  The  latter  is  entered  with  a negative  sign  so  that  terms 

such  as  m.x.-F  from  Section  6.3  are  properly  compiled, 
i 

The  loads  associated  with  the  six  principal  reference  axis  degrees  of 
freedom  are  listed  in  Table  6-5*  The  propeller,  tail  rotor,  main  rotor, 
and  engine  are  assumed  to  have  shafts  parallel  to  the  fuselage  reference 
axes.  The  transfer  of  the  aerodynamic  loads  from  tail  rotor  axes  with 
origin  at  hub  center  and  parallel  to  the  fuselage  reference  axes  is  shown 
in  the  table.  The  propeller  is  similar.  The  fuselage  aerodynamic  loads 
are  assumed  to  already  be  in  the  fuselage  reference  axes.  Further 
development  of  the  main  rotor  blade  component  loads  is  in  Section  6.6.1 
and  the  aerodynamics  for  all  rotors  and  fixed  surfaces  is  left  to 
Section  7. 
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TABLE  6-1+  - Continued 
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where  the  following  weight  and  balance  relations  hold. 
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TABLE  6-5-  PRINCIPAL  AXIS  GENERALIZED  FORCES 
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TABLE  6-5  - Continued 
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The  main  rotor  shaft  loads  are  simply  found  from  the  blade  root  loads . 
The  approximate  equations  are: 
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matrix  accounts  for  rotor  gimballing,  the  transfer  between 


stationary  and  rotating  axes,  and  the  switch  between  Z up  and  Z down 


axes . The 


tr- 


BLn 


matrix  accounts  for  the  individual  blade  azimuth. 


These  matrices  were  developed  in  Section  5- 5- 5-  Z is  the  amount  of 

OBLn 


undersling  (negative  value). 


6.10  SWASHPLATE  EQUATIONS 


6.10.1  Use  as  a Pure  Swashplate  or  a Feathering  Feedback  Control  Gyro 


The  swashplate  converts  the  nonrotating  system  (hub  axis)  collective  and 
cyclic  blade  commands  into  first-harmonic  Fourier  components  in  the 
rotating  system  (rotor  coordinates).  For  a conventional  swashplate  system 
the  command  signal  operates  open  loop  from  input  to  blade  feathering. 

The  swashplate  formulation  used  in  REXOR  handles  this  case,  but  this 
original,  and  prime,  use  is  for  the  external  control  gyro  with  feathering 
feedback.  This  is  a closed-loop  system  where  the  command  input  torques 
the  control  gyro,  which  in  turn  feathers  the  blades.  Response  rates  of 
the  principal  axes  and  rotor  moment  via  hub  geometry  form  a control  system 
loop  closure.  Both . systems , as  well  as  the  swashplate  element  of  the 
isolated  control  gyro,  are  handled  merely  by  changing  the  system  constants 
input . 
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6.10.2  Partial  Derivatives 

The  swashplate  partial  derivatives  are  readily  obtained  from  Section  5*3-5* 
Using  matrix  notation 
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Also : 
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Since  the  swashplate  axis  is  directly  referenced  to  the  principal  (hub) 
set,  the  above  derivatives  are  complete.  The  lack  of  translation  to 
angular  derivatives  is  explained  by  the  parallelogram  linkages  used  with 
swashplates  to  isolate  the  collective  and  cyclic  inputs.  The  terms  left 
out  of  the  matrix  indicate  that  the  swashplate  does  not  have  a yaw  degree 
of  freedom. 

The  reader  should  be  aware  that  the  angular  notai  ion  <J> , 0,  and  \ t>  have 
two  meanings,  depending  on  whether  they  are  in  the  numerator  or  the 
denominator  of  the  partial.  The  numerator  is  the  displacement  of  the 
mass  element  with  respect  to  the  principal  reference  axis,  whereas  the 
denominator  is  the  degree  of  freedom  incremental  variable. 

Swashplate  motions  pick  up  large  inertia  loads  from  the  rotor  due  to  blade 
feathering.  Partials  relating  feathering  to  swashplate  motions  are 
assembled  by  first  relating  the  feathering  motion  in  the  rotating  system 
with  feathering  in  the  stationary  system: 
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From  Section  5-5.8,  equations  relating  swashplate  motions  to  the  stationary 
feather  angles  give 


3A1S 

3A1S 

3*sp 

39SP 

3B1S 

3B1S 

3*SP 

39SP 

= — T 
e tp 


■at 


^PH  |9SP 


Sln*PH 

COS*PH 

.C0S,,PH 

-3iD*w. 

(6-3014 ) 


(6-305) 


(6-306) 


(6-307) 


The 


K] 


matrix  does  not  follow  the  conventional  Euler  angle  notation 


since  a desire  existed  to  define  ^ as  the  angle  the  pitch  horn  to 

rH 

pitch  link  attachment  point  leads  the  blade.  The  overall  derivatives  can 
be  put  together  as: 


3<j> 


Fn 


dp 


3*Fn  3A1S  + 3*Fn  3B1S 


SP 


3A1S  3*sp 


SB  Sdi 
IS  VSP 


(6-308) 


^Fn  8*Fn  3A1S  , 3*Fn  3B1S 
36SP  ' 3Ais  39sp  3Bis  39SP 


(6-309) 


3<f> 


Fn 


SP 


9 


!!ln  3^S  + !lFn  ^is\^0_ 

30o  3Ais  38o  3Bis36o/3Zpp 


(6-310) 


6.10.3  Generalized  Masses 

Table  6-6  presents  the  generalized  masses  which  couple  the  swashplate 
motions  with  one  and  another,  with  the  rotor,  blade,  and  body  degrees  of 
freedom.  The  table  uses  summations  of  the  blade  that  are  described  in 
detail  in  Section  6.6.1. 
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The  feathering  moment,  M^  , is  taken  to  be  composed  of  blade  loads, 
friction,  and  tension-torsion  pack  loads. 
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The  detailing  of  M^  , feathering  moments  due  to  blade  loads,  is  accom- 
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plished  in  Section  6.6.4.  The  friction  load,  M , follows  the  func- 
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tion  shown  in  Figure  6-3.  By  reducing  ‘(j>  to  near  zero,  stiction 

Fn » BK 


Fn  BK 

is  obtained.  Otherwise,  if  Vf>  is  large,  the  ratio  — 2 

Fn  ,BK  $ 


Fn  ,BK 


determines  the  amount  of  viscous  friction. 


The  first  part  of  the  tension-torsion  pack,  M , has  a simple  spring 

TTFn 


contribution  of 


KTT^Fn  * ^Fn^) 


(6-350) 


from  the  torsion  spring  K_m;  <{> 

TT  Fn„ 


being  the  blade  ar gle  for  zero  spring 


load.  The  tension  contribution  from  the  tension-torsion  pack  requires  a 
knowledge  of  the  locations  of  the  inboard  and  outboard  feather  bearings 
and  the  inboard  and  outboard  pin  ends  of  the  tension-torsion  pack.  The 
locations  of  the  bearings  are  described  in  Section  5-5-5>  The  inboard  end 
of  the  pack  is  the  sum  of  the  static  and  blade  bending  contributions: 
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he  static  position  formulation  follows  the  development  of  Section  5.5.5 

inboard  of  Xgw  (blade  attack  point  to  movable  hub)  except  that  geometric 
twist  is  not  included. 
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Knowing  the  location  of  the  tension-torsion  pack  and  the  feather  bearings, 
the  total  contribution  from  the  pack  is 
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where  for  each  blade: 
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If  no  stretch  is  assumed,  the  £'s  can  be  treated  as  constants.  Direc- 
tional cosines  are  used  for  convenience.  The  negative  value  of  the  tension 
in  the  feather  moment  equation  indicates  a tension-torsion  pack  reaction 


force  acting  through  the  small  arms  X^  - X^,  - Y^,  and  - ZJB 


to  produce  a feathering  moment. 


The  tension  is  equal  to  the  component  of  the  root  shear  lying  along  the 
tension-torsion  pack  axis.  The  lateral  components  of  the  shears  reacted 
at  the  feather  bearing  are  ignored.  Using  directional  cosines 
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The  remaining  portions  of  the  generalized  force  are  the  potential  energy 
and  dissipation  functions.  First  consider  the  angular  potential  energy 
terms  which  model  the  swashplate  tilt  spring  rate.  This  spring  rate  has 
a center  dead-hand,  an  operating  range  spring  rate,  and  a high  spring  rate 
to  simulate  a travel  limit  stop. 


Consider  the  normal  operating  range  spring  rate  first.  The  swashplate 


springs  are  defined  in  control  axes  (Figure  6-5)  as  K,  and  KQ  (can 

*SP  9SP 


be  unequal  in  size).  To  find  the  elastic  spring  loads,  the  swashplate 
motions  are  first  found  in  control  axes  as: 
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The  geometric  interpretation  of  tpr  is  shown  in  Figure  6-5. 


A shaft  bending  to  swashplate  coupling  ratio,  R , is  allowed.  This 

orb 


coupling  exists  for  the  unloaded  swashplate  in  much  the  same  way  as  the 
blade  bending  to  feathering  coupling  exists  for  the  unloaded  pitch  horn 
and  defines  a modal  component. 


Taking  the  swashplate  deflections  i;  the  control  axes,  subtracting  control 
inputs  <J>£  and  6^,  and  using  the  inverse  transform,  the  swashplate  spring 

terms  in  swashplate  axes  become: 
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where  the  subscript  (l)  is  used  to  distinguish  these  values  (used  in  sub- 
sequent logic  calculations)  from  the  final  expressions  developed  below. 
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Substituting  for  the  swashplate  motions  in  terms  of  the  swashplate  axes 
and  rearranging: 
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where 
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Note : 


: [ksp  1 


is  a symmetric  matrix  of  constants. 


The  spring  constants,  K,  , KQ  , exist  as  springs  in  the  usual  sense 

^SP  8SP 

when  the  swashplate  serves  strictly  as  a swashplate.  When  the  swashplate 
is  used  as  a control  gyro,  the  springs  may  be  composed  of  springs  operating 
in  the  usual  sense  plus  a preloaded  on  center  spring, which  counteracts  the 
other  spring  (positive  spring)  when  the  system  is  viewed  from  the  swash- 
plate. The  preloaded  spring  is  called  the  negative  spring.  With  this 
arrangement,  it  is  possible  to  have  zero  spring  rate  of  the  swashplate 
motion  while  having  a very  stiff  input  spring  rate. 


to  account  for  a limit  travel  stop.  The  limit  deflection  for  the 
swashplate  is 


/ 2 2 v'"'^ 

(^SP  + 9SP  ) “ 6s, 


SP 


( 6-37** ) 


where  <5^  is  the  circular  stop  swashplate  deflection  limit. 


The  angular  damping  term  is  analogous  to  the  spring  load : 
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where 


M 


has  the  same  formulation  as 


[Ksp] 


Control  friction  is  treated  as  having  rotating  and  nonrotating  components . 
The  rotating  component  has  already  been  discussed  as  part  of  the  feathering 
moment.  The  nonrotating  component  is  applied  to  the  swashplate.  It  has 

the  formulation  shown  in  Figure  6-3  with  a change  ir  labels  such  that 


Fn 


is  either  *sp  or  e,.p,  and  is  either  or 


The  vertical  potential  energy  term  is  described  as: 
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F is  a constant  to  center  the  gyro  springs. 
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The  spring  rate  is  taken  to  be  K_  out  to  deflection  Zn  and  K 

^ZSP  1SP  ^ZSP 

beyond . 

A simple  coupling  from  the  rotary  dampers  gives  the  vertical  dissipation 
function . 


: z-rc  d>d)-Rc  e e 

zsp  sp  nz<t>  ^gp  9SP  ?SP  ze  egp  sp  sp 


(6-378) 


where  R^,  R_,0  are  coupling  ratios.  Note  the  effect  of  vertical 
motion  on  the  swashplate  tilt  loads  through  the  rotary  dampers  is  assumed 


To  correlate  with  flight  test  records  and/or  to  force  the  swashplate 
vertical  response  to  cross  the  spring  rate  changeover  a force  offset  con- 
stant is  used.  Introducing  this  constant  into  the  swashplate  vertical 
degree  of  freedom  equation  line,  causes  the  variables,  primarily  the 
swashplate  vertical  motion,  to  shift  and  rebalance  the  equation. 

6.10.5  Control  Inputs 

The  swashplate  angular  input  is  controlled  by  the  pilot's  lateral  stick  for 
either  a pure  swashplate  or  use  as  a control  gyro.  The  input  torques  are: 


Vc 


Vc 


(6-379) 


The  inputs  are  aligned  with  the  control  axis  (Figure  6-5). 

Note  the  equivalence  of  forms  in  terms  of  angular  commands  ®c’  or 

longitudinal  stick  (aft)  Xg  and  lateral  stick  (right)  Y^. 

The  controls  are  frequently  linked  to  the  swashplate  through  actuators 
which,  as  a first-order  approximation,  can  be  simulated  by  a first-order 
lag.  For  illustration,  consider  the  longitudinal  stick.  Its  rate  becomes: 


<v-v 


(6-380) 
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The  position  is  found  by  a time  integration.  The  pilot's  input  is  Xp  and 
is  the  time  constant.  The  actuator  may  be  rate  limited: 


X„  s X„ 
C C, 


(6-381) 


or  displacement  limited: 


X„  2 X„ 
C C, 


(6-382) 


Similar  formulations  can  be  applied  to  Y , <j>Q,  9q.  For  the  vertical  or 
collective  system,  the  command  input  is  = -0^  e and  the  actual,  true 
value  is  Z = -0  e in  terms  of  the  collective  angle  and  the  pitch  horn 

Ox  O 

arm.  The  vertical  system  is  typically  fully  powered,  and  0 is  modeled 

UG 

similarly  to  the  technique  just  described  for  Xc< 

6.11  CONTROL  GYRO  E^r  -u'IOHS 
6.11.1  Introduction 

As  developed  briefly  in  Section  6.10.1,  REXOR  models  three  cyclic  control 
configurations.  The  first  two  type^  are  the  conventional  swashplate  and 
the  feathering  feedback  external  control  gyro.  Both  of  these  systems  use 
the  same  modeling  in  REXOR;  the  type  being  set  by  the  selection  of  the 
system  constants.  The  third  configuration  is  an  isolated  direct  flap 
feedback  control  gyro  which  operates  a swashplate,  as  modeled  in 
Section  6.10,  via  irreversible  actuators. 

This  configuration,  and  the  associated  degrees  of  freedom,  <()_,,  0„,  are 

G G 

defined  somewhat  differently  in  that  no  generalized  mass  couplings  are 
assumed  with  the  other  degrees  of  freedom.  The  generalized  masses  involved 
are  simply  the  diametral  moments  of  inertia  of  the  control  gyro. 

This  simplification  is  possible  because  the  inertia  of  the  gyro  is  so  small 
the  principal  axis  loads  do  not  include  the  control  gyro  contributions. 

This  leads  to  zero  generalized  mass  coupling,  and  advantage  is  taken  of 
this  fact  to  decouple  the  gyro  degrees  of  freedom.  Being  decoupled,  the 
control  gyro  can  be  handled  apart  from  the  main  stream  of  the  computations. 
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The  following  development  also  differs  from  other  sections  in  that  the 
control  gyro  interface  with  the  other  rotorcraft  elements  (feedback  and 
phasing)  are  treated  as  separate  subsections. 

6.11.2  Partial  Derivatives 

Partials  are  required  to  relate  the  motions  in  the  instantaneous  gyro  axes 
to  the  gyro  Euler  angle  degrees  of  freedom.  The  desired  partials  are 
extracted  from  the  rotary  velocities  relationship  of  Section  5- 5-7  and  the 
geometry  shown  in  Figure  5-9*  They  are: 
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The  values  of  the  above  partials  are  the  same  for  axes  rotating  with  the 
rotor  or  nonrotating  axes  attached  to  the  gimbal.  As  shown  in  Sec- 
tion 5*5.7»  this  is  due  to  the  rotor  axes  being  inertial  axes  defined  to 

be  momentarily  attached  body  axes  at  i|>_  = 0 and  therefore  coincident  with 

gimbal  axes. 

The  difference  between  the  inner  and  outer  gimbals  is  ignored  and  the  gim- 
bal supporting  the  gyro  rotor  is  taken  to  contain  all  the  inertia. 

6.11.3  Generalized  Forces 


The  control  gyro  loads  result  from  inertia,  springs  and  dampers,  pilot 
input  and  rotor  flap  feedback.  The  detailing  of  the  feedback  is  extensive 
enough  for  a separate  section.  The  net  loads  are: 
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where  subscript  R is  rotating  component  only,  and  NR  is  r.’nrotating 
component  only. 
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Stiction  is  allowed  in  the  feedback  load,  and  viscous  friction  can  be 
equivalenced  to  damping. 

The  equations  can  be  solved  for  "<f>  and  0 by  substituting  in  the  expres- 

G G 

sion  for  angular  accelerations  derived  in  Section  5 - 5 - T • These  are: 


where : 
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Substituting  the  values  of  the  partials  and  the  above  into  the  gyro 
equations  of  motion  and  rearranging  yields: 
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Expanding  the  inertia  load,  the  general  formulation  for  the  moments  gen- 
erated by  a body  along  principal  axes  are: 
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where  the  moments,  the  velocities,  and  the  accelerations  are  all  along  the 
principal  axes.  If  the  body  is  axially  symmetric  about  the  Z axis, 


My  = IXX(P  + 9.r)  = Izz(p  + qr)/2 
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My  = Iyy(q  - rp)  = IZ2(q  - rp)/2 
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Next,  the  motions  are  conveniently  described  as  composed  of  a nonrotating 
component  plus  a spin  \p. 


In  addition,  the  body  motions  are  assumed  referenced  to  inertia  axes  which 
are  instantaneously  attached  to  the  body  at  0=0.  Specifically  for  the 
gyro  rotor  the  relations  are: 
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Therefore,  for  the  gyro  gimbals  which  do  not  rotate,  the  moments'  are: 

V 


Next,  recalling  that  all  the  gimbal  inertia  is  in  the  gimbals  that  hold 
the  rotor , we  have  for  the  gyro  rotor : 
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The  contribution  from  spin  thrust  appear  as  add-on  terms.  Collecting 
terms : 
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The  NR  notation  was  dropped  from  p.  and  q since  p_  = p_  and 
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In  = 'In  iro*  The  control  gyro  rotation  is  assumed  to  be  independent  of 
G G , NR 


rotor-engine  speed  variations,  and  driven  through  a constant  speed 
coupling.  These  assumptions  are  reasonable  even  with  rotor  shaft  driven 
configurations  because  the  coupling  of  drive  torque  equation  into  the 
required  roll  and  pitch  equations  is  a second-order  effect.  The  yaw 
moment  reduces  to 
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The  spring  and  dampers  are  taken  to  be  defined  in  control  axes.  The 
formulation  follows  the  corresponding  formulation  for  the  swashplate: 
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The  development  of  the 


follows  precisely  the  formulation 


presented  for  the  swashplate  in  Section  6.10.5  or  can  he  directly  sup- 
plied as  input.  Note  the  cross  coupling  terms  only  exist  if  the  spring 

and  damper  are  not  symmetric . The  stick  positions  X and  Y are  the 

c e 

equivalent  position  output  of  first-order  servos  modeled  in  the  manner  of 
Section  6.10.5. 

The  gimbals  are  allowed  to  he  unbalanced.  A simple  formulation  couples 
huh  vertical  motion  into  the  gyro,  and  is  included  in  the  previous 
generalized  force  expressions. 
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Gyro  tilt  transforms  are  ignored.  The  distances  X^  and  are  positive 

if  the  imbalance  mass  is  to  the  right  or  forward  of  the  gyro  spin  axis. 


6.11.4  Control  Gyro  Feedback 


The  feedback  is 
concept,  and  is 
not  precessing, 
the  gyro  spring 
proportional  to 
is  proportional 


highly  important  to  the  direct  feedback  gyro  control 
detailed  with  care.  For  steady  flight,  when  the  gyro  is 
the  pilot  input  counters  the  moments  from  the  feedback  and 
, The  spring  contribution  is  minor.  Since  the  feedback  is 
blade  flap  displacements,  the  resulting  rotor  shaft  moment 
to  the  pilot's  stick  displacement. 


Typically,  to  achieve  a feedback  proportional  to  the  rotor  shaft  moment, 
a feedback  lever  is  mounted  part  way  out  on  the  fixed  hub.  Displacement 
and  tilt  of  the  lever  mount  with  blade  bending  causes  the  tip  of  the  lever 
to  move,  driving  a rod  that  is  attached  to  the  gyro  feedback  spring.  The 
displacement  of  the  tip  of  the  lever  is : 


ZJ  ZRM  “ Z'rm 


(xrm-xj) 


Where  the  subscript  RM  and  J refer  to  the  lever  moment  and  tip,  respec- 


tively. The  lever  spanwise  length  is  then  - Xj.  The  linear  and 


slope  displacements  of  the  lever  mount  are  given  as: 


Z = y*  8"RM  . 

n ^ ^ SA  mn 
m=l  mn 


3 3Z\ 


-E 


for  blade  n due  to  modes  1,  2 and  3.  The  displacement  does  not  depend 
on  feathering,  and  the  mount  is  taken  to  be  at  the  same  waterline  as  the 
hub  origin.  Polarity  of  the  motion  needs  to  be  observed.  For  instance, 
the  construction  of  the  lever  for  the  AH-56A  helicopter  caused  the  slope 


displacement  contribution  to  Zj  to  be  greater  than  Z^  . Therefore, 

n n 


the  tip  of  the  lever  went  down  as  the  blade  flapped  up. 


A feedback  spring  preload  is  used  to  preclude  deadband.  Friction  exists 


and  is  modeled  by  letting  the  point  reverse  direction  an  incremental 


amount,  AZ. 


LIMIT 


LIMIT 


(6-1+15) 


'Phis  acts  without  changing  the  applied  gyro  moment. 


Mpg  is  the  gyro 
FR 


stiction  moment  and 


the  feedback  spring.  The  feedback  spring 


displacement  becomes: 


. {f  \ 


(6-1+16) 
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where : 


(fi, 

\Jo  n 


= SIGN 


LIMIT 
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dt  < AZ 

J . 
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Jo  n 


dt  (6-1+17) 


The  action  of  stiction  on  the  feedback  displacement  during  a reversal  of 
velocity  is  shown  in  Figure  6-6. 


Point  J 

r 

the  blade  reference  and  a distance  from  the  rotor  shaft  center.  The 


Jn  on  the  feedback  lever  is  located  at  an  azimuth  ip  ahead  of 

FB 


point  in  fuselage  axis  becomes 
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(6-1*18) 


where 
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(6-1+19) 


The  other  end  of  the  feedback  spring  is  grounded  to  the  gyro  at  a radius 
°f  XFBG  whlch  leads  the  blade  by  the  angle  ik  . Hence, 

r .Dtr 


FBG 


FBG 


V = 


"FBG 


m "It  r It  T "It  I”  "]r 
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(6-1+20) 
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The  feedback  spring  force  is  assembled  as: 


(FFBn)p  ~ KFB(ZFBn  ~ *Fb) 


(6-423) 


The  length  £,  serves  as  the  preload  constant.  It  does  not  affect  the 

I*  JD 


cyclic  spring  feedback  to  the  gyro.  The  vertical  displacement  Z of 

r BG 

n 


the  lower  end  of  the  spring  is  excluded.  It  produces  a force  proportional 
to  gyro  tilt  and  this  force  is  best  considered  to  be  a component  of  the 
gyro  spring  value. 


The  feedback  moment  is  the  spring  force  times  the  gyro  attachment  arm. 
The  moment  becomes : 


\ ^ 'l 

r Y F -V 

FBG  FB 
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(6-424) 


which  is  taken  about  the  gyro  Euler  axes . 
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To  close  the  gyro-swashplate-rotor  "blade  loop,  a relation  is  needed 
between  gyro  motions  and  swashplate  motions.  Assuming  a power  actuator 
between  the  gyro  and  swashplate  gives : 


GSP  YSP, 


(6-H25) 


tgsp  6SP, 


(6-1*26) 


The  subscript  C refers  to  the  command  or  control  value. 


The  actuator  variables  can  have  rate  and  displacement  limits  in  the  manner 
described  in  Section  6.10.5. 


6.11.5  Summary  Description  of  Phase  Angles 


A number  of  phase  angles  were  introduced  in  this  section  and  Section  6.10 
on  the  swashplate.  Figure  6-7  shows  the  relationship  between  all  these 
phase  angles  as  they  apply  to  the  Direct-Flap  Feedback  Control  System. 

The  values  indicated  are  for  the  AH-56A  Cheyenne  helicopter.  The  pilot  is 
assumed  to  have  pulled  the  stick  back,  precessing  the  gyro  until  the  flap 
feedback  was  sufficient  to  react  against  the  pilot  input.  The  vehicle  is 
assumed  fixed,  otherwise  it  would  pitch  up  and  relieve  the  flapping.  The 
gyro  spring  load  is  minor  and  is  not  shown.  Swashplate  to  blade  linkage 

is  made  with  a trailing  pitch  horn,  hence  is  180  degrees  minus  the 


actual  pitch  horn  angle  to  account  for  the  reversal  in  feathering  going 
from  a leading  to  a trailing  pitch  horn.  The  pitch  horn  phasing  is  chosen 
to  place  feathering  at  an  azimuth  that  will  produce  flapping  at  the  desired 
direction,  in  this  case  up  flap  forward  for  a nose  up  pitching  moment. 


Note  the  gyro  precesses  up  in  an  axis  perpendicular  to  the  applied  moment 


F^  , but  that  the  swashplate  displacement  is  in  the  same  axis  as  the 


applied  moment  F^  . The  swashplate  control  acts  through  very  stiff 


spring  rate  which  prevents  the  rotating  swashplate  from  responding  in  a 
low-frequency  precession  mode. 
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Figure  6-7.  Phasing  For  Direct-Flap  Feedback  Control  System 


6.12  ENGINE  EQUATIONS 

6.12.1  Rotor  Azimuth  and  Rotation  Rate 

The  program  allows  a variation  of  rotor  speed  in  maneuvers  due  to 
variations  in  the  torque  required  by  the  various  rotors  and  in  the  torque 
supplied  by  the  engine.  The  dynamic  system  rotates  as  a rigid,  geared 
unit.  That  is,  the  shafts  are  not  allowed  elastic  windup.  The  main 

rotor  speed,  ^ and  hence  the  engine  speed,  is  referenced  to  the  fuse- 
R 

lage,  and  not  to  inertial  space.  The  displacement  i|j  is  the  azimuth  of 

R 

the  number  one  blade. 

6.12.2  Engine  Model 

Figure  6-8  illustrates  the  engine  model  used  in  the  program.  The  figure 
also  plots  typical  engine  torque  characteristics.  The  model  represents 
the  first-order  lag  power  response  characteristics  of  the  free  turbine 
powerplants  commonly  used  in  rotorcraft  applications . 

Being  a perturbation  model,  the  engine  is  referenced  to  its  trim  position. 
The  change  in  engine  torque  in  a maneuver  is 


M - M 

^ENG  ** 


F,NG,TRIM 


9MENG  r 3MENG  / • 

9*gen  gb<  ' 8W,  ( 0,0 


ENG  y ENG, TRIM 


(6-427) 


where  0 < < M . The  zero  limit  occurs  if  the  overrunning 

^ENG  iSNGjMAX 

clutch  disconnects  the  engine  in  the  'transition  to  autorotation.  The 
maximum  value  corresponds  to  the  engine  shaft  torque  limit. 

The  gas  generator,  speed,  , is  a degree  of  freedom.  It  is  considered 

a "secondary"  degree  of  freedom  in  that  the  coupling  through  the  general- 
ized masses  with  the  "primary"  degrees  of  freedom  can  be  neglected.  An 
equation  for  the  generation  speed  can  be  supplied  from  its  torque 
characteristics : 

TGEN  ^GEN  + CGEN  ^GEN  ~ " KENG1  ^ENG  " KENG2(’^ENG  “ ^ENG,TRIm)  ( 6-1+28 ) 


5 VI*? 
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... 


The  terms  on  the  left  represent  acceleration  inertia  torque  and  steady- 
state  torque.  On  the  right,  the  fuel  control  causes  torque  to  he  added  if 


the  engine  speed  drops  below  the  trim  value.  The  ipz 


tern  exists  since 


the  control  is  modeled  with  simple  lag.  Restating  this  equation,  using 
rotor  speed  and  a generator  time  constant,  gives: 


^GEN  ~ Si  ~ ^^R  ~ ^R,TRIM^ 


(6-U29) 


where  ij> 


GEN  r. 


is  the  order  of  a second. 


The  engine  droop  characteristic  can  be  used  to  size  the  engine  constants. 
With  4^™.,  = <L  = 0,  substituting  the  generator  equation  into  the  engine 

equation  and  rearranging. 


( XArac;)  ° (%<;)_  3 ('wig). 


8*oe»  *** 


(6-1*30) 


Only  A incremental  changes  are  of  interest.  The  bracket  subscripted  R 
indicates  the  torque  is  determined  at  the  rotor  speed  and  includes  the 
engine  gear  ratio.  The  term  on  the  right  is  the  static  droop  line  shown 
in  Figure  6-8.  This  plot  also  geometrically  interprets  the  partial 
derivatives  on  the  left. 

The  generator  speed  is  not  given  a reference.  Its  value  is  zero 

when  trim  is  completed. 

6.12.3  Partial  Derivatives 

Partial  derivatives  relating  linear  motions  of  the  blade  elements  to  shaft 
rotations  are  needed: 


!Ible 

l S*R  J 


^BLeI  / _^Fn\ 

13*ri  L W ) 


(6-1*31) 
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Shaft  rotation  not  only  involves  blade  root  rotation  ij>  , but  also 


feathering  motions.  The  feathering  partial  is  obtained  by  differentiating 
the  feather  angle  equation  in  Section  5-5.8: 


IT1  - A1S  sln(*BLn  + V - B1S  °°s(*BLn  * *R> 


(6-1+32) 


6.12.1+  Generalized  Masse? 


The  engine  degree  of  freedom  couples  in  with  every  other  degree  of  freedom 
except  those  for  shaft  bending.  Equations  for  the  engine  generalized 
masses  are  given  in  Table  6-7- 


TABLE  6-7  - Continued 


M,  . = M, 


3<j> 

TFn 

+ M. 


<P nA  1.  A 9i}>  4>  A 

R mn  BLn  mn  R Fn  ran 


(6-itUx) 
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(6-443) 


6.12.5  Generalized  Forces 

Only  one  generalized  force  is  needed: 


(6-444) 


The  inertia  component  of  the  propeller,  the  tail  rotor  and  the  engine  are 
formulated  in  Section  6.9.4.  The  aerodynamic  description  follows  in 
Section  J.  The  positive  direction  of  these  moments  are  right  roll,  nose 
up,  and  yaw  right  per  fuselage  axes. 


7.  AERODYNAMICS 


7.1  INTRODUCTION 

Other  than  gravity,  the  external  loadings  acting  on  the  REXOR  equations  of 
motion  can  he  traced  to  aerodynamic  sources.  The  following  subsections 
trace  the  source,  nature  and  use  of  these  aerodynamic  loads. 

7.1.1  Aerodynamic  Forces  Producing  Surfaces  Considered 

The  aerodynamic  loads  considered  in  REXOR  are  divided  into  the  categories 
of  (l)  associated  with  the  main  rotor,  or  (2)  the  rest  of  the  rotorcraft 
(nonrotating  surfaces,  tail  rotor  and  propeller).  In  view  of  the  stated 
objectives  of  REXOR  the  program  development  emphasis  is  on  the  main  rotor. 

7.1.2  Use  of  Forces  Generated 

As  mentioned,  the  aerodynamic  loads  are  in  essence  the  external  forcing 
functions  of  the  equations  of  motion.  Generally  the  developed  loads  are  in 
the  axis  of  the  apparent  air  velocity  of  the  loaded  element.  Thus  trans- 
formations are  required  to  put  the  loads  into  the  reference  axes  of  the 
equation  of  motion  considered. 


7.2  MAIN  ROTOR 
7.2.1  Overview 

To  generate  a main  rotor  model  with  sufficient  detail  to  do  dynamic  inves- 
tigations, a reasonably  good  quality  aerodynamics  presentation  is  re- 
quired. To  this  end  a table  lookup  of  blade  section  properties,  multifunc- 
tion inflow  model,  quasi-steady  aerodynamics,  and  dynamic  stall  are  used 
in  REXOR. 

7. 2. 1.1  Blade  Flow  Field 

As  developed  in  the  following  subsections,  the  instantaneous  blade  air- 
flow is  the  inertial  velocity  of  the  blade  element.  This  velocity  in- 
cludes the  motion  of  the  principal  reference  set  and  the  motion  of  the 
blade  element  with  respect  to  the  principal  reference  set.  The  cal- 
culation assumes  the  airmass  is  at  rest,  which  is  reasonable  for 
dynamics  investigations. 

7.2.1. 2 Air  Pressure  and  Angle  of  Attack 

The  dynamic  pressure  used  for  these  calculations  is  based  on  sea 
level  standard  density.  The  loads  are  ratioed  to  the  actual  air  density. 


The  angle  of  attack  is  the  sum  of  geometric  pitch  angle  and  the  instan- 
taneous air  velocity.  The  rate  of  angle  of  attack  is  also  calculated 
and  used  for  the  transient  blade  aero  loads.  Sections  7*2. 3. 5 and  7. 2. 3. 4. 


7.2.1. 3 Forces  and  Moments  Produced 


The  steady  blade  loads  are  produced  from  the  air  velocity  components 
of  Section  7-2. 3.1  and  the  coefficient  data  (C^,  Cp,  Cj^)  of  Section 
7.2.4.  The  transient  lift  and  moment  effects  are  developed  in  Sec- 
tions 7*2.3. 3 (quasi-steady  aerodynamics)  and  7*2. 3.4  (dynamic  stall). 


7.2.2  Concept  of  Rotor  Inflow  Model 


The  main  rotor  inflow  model  used  in  REXOR  is  based  on  the  air  flow  incident 
upon  the  rotor  disc  plus  the  air  velocity  imparted  due  to  momentum  exchange 
due  to  integrated  blade  span  loading.  This  is  to  be  contrasted  with  a 
formulation  which  tracks  the  rotor  blade  positions  and  the  attendant  trail- 
ing vorticies 


The  incident  air  flow  is  the  inertial  velocity  of  the  rotor  coordinates, 
and  is  directly  available  from  the  preceeding  mechanical  development.  How- 
ever a number  of  assumptions  need  to  be  stated  and  utilized  to  arrive  at 
the  induced  velocity  component  of  the  inflow  model. 


7. 2. 2.1  Induced  Velocity  Assumptions 


1.  Only  the  vertical  downwash  and  its  variations  radially  and 
azimuthally  over  the  rotor  disk  are  considered.  Induced 
swirl,  and  lateral  downwash  components  are  considered  minor 
and  therefore  neglected. 


2.  Downwash  effects  due  to  unsteady  aerodynamics  are  not  treated 
here  as  an  overall  effect,  but  as  a blade  segment  condition 
in  Section  7* 2. 3- 3. 


Rotor-induced  flow  distribution  in  hover  and  forward  flight 
is  patterned  after  Reference  10.  This  reference  assumes  a 
uniform  loading  in  hover.  Figure  7-1,  from  Referenc°  11, 
shows  this  distribution  compared  with  typical  loading  and  a 
triangular  loading  model.  Figure  7-2  from  Reference  10  shows 
the  theoretical  induced  velocity  distribution  in  forward 
flight  as  a consequence  of  a uniform  hover  distribution. 

This  data  is  fitted  to  slopes  or  a longitudinal  skew  as  a 
function  of  speed  in  REXOR.  Lateral  distribution  remains 
uniform  in  accord  with  Reference  12, which  corrects  the 
lateral  distribution  work  of  Reference  10. 


A variation  in  lateral  and  longitudinal  induced  velocity  is 
included  to  account  for  roll  and  pitch  aerodynamic  shaft 
moments. 


LOCAL  DISK  LOADING 


5.  Lifting  line  theory  correction,  rotor  tip  loss,  and  root  cut- 
out effects  are  ignored. 

6.  Transient  effects  are  simulated  by  a single  time  lag. 

7. 2. 2. 2 Steady  State  Values 

The  starting  point  for  determining  the  do\mwash  is  momentum  theory  as 
applied  to  an  elementary  area  dA: 


dT  = (f)  2Mi  - (PViMR  “) 


= P\/UH2  + VH2  + (VH  - Wi>2  “ *1 


(7-1) 


The  thrust  increment  is  dT,  dm/dt  is  the  flow  of  air  through  the  rotor 
disk  with  resultant  velocity  V3-MR,  p is  the  air  density  and  w^^  the 
downwash  velocity.  The  velocities  are  taken  in  hub  coordinates  and  no 
effort  is  made  to  account  for  rotor  tilt. 

The  thrust  expression  above  is  used  to  define  the  following  induced 
velocity  components. 

• Average  component,  w^ 

• Longitudinal  variation  with  pitching  aerodynamic  moment, 

<HMR 

• Lateral  variation  with  roll  aerodynamic  moment,  p 

,blMR 

The  downwash  velocity  becomes 


”i  WiMR  + r ?iMR  C0S  *8  + r PiMR  sin  *R 


(7-2) 


The  coefficients  can  be  evaluated  by  equating  the  thrust  and  moment 
values  for  the  main  rotor  equations  to  the  integrals  of  the  momentum 
expressions  at  hand.  First  consider  the  thrust  expression.  The 
evaluating  task  can  be  reduced  by  employing  some  boundary  conditions. 
For  rotor  thrust  only  (no  moment),  q^R  and  piMR  = 0.  A convenient  ex- 
pression for  the  elementary  area,  dA,  is  shown  in  Figure  7-3.  While 
radial  anulii  would  serve  for  thrust  integration,  the  form  selected  is 
particularly  suited  for  the  moment  expressions. 


■ 


For  the  average  rotor  thrust, 


R 

,,H  = T ' / "V""2  + VH2  * H - \m)  2 (2Vr2  - r2  dr  ) 2 ”i 


(7-3) 


A further  assumption  is  required  to  solve  the  square  root  of  this 
pressxon  and  the  corollary  momentum  equations. 


For  forward  flight 


Wi  <<  ViMR  “\/Uh  + \ + (WH  - wiMR)  = (constant) 


(7-1*) 


Completing  the  integration  gives 


\/PlR  ViMR2wiMR 


(7-5) 


exo  consider  the  case  of  no  rolling  moment;  ie.,  only  thrust  and 
pitching  moment.  Figure  7-3  is  used  with  the  incremental  strip  con- 
sidered to  he  right-left  oriented  so  that  all  equal  values  of  qiMR  are 
integrated  at  once.  lMR 


Then,. 


'He,! 


r P\l\  + VH2  + (VH  - "if 


- r q 


imr) 


(2  vR2  - r2  dr)  2 


(ViMR  + r Hmr) 


(7-6) 


‘ e *•  ViME  2 qiMR 


(7-7) 


Likewise  for  rolling  moment,  and  using  fore-aft  increment  strips 
gives 


V„ = 0 v 2 Pil 


(7-8) 


Note  the  subscript  A on  F_.  , M . , and  My.  denot 

\R,H  ,H 

aerodynamic  component  only  of  main  rotor  loads  in  hub  axes. 


denotes  the 


The  foregoing  expressions  are  now  developed  for  hovering  and  low- speed 
flight.  In  this  condition, 


2 2 
u + v << 
H H 


("h  - ViMn)‘ 


(7-9) 


Integrating  gives 
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(7-10) 
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) (7-3 


Vh=  ^ 2 


WiMR  ^WH  ~ WiMR^ 


(w  - w )‘ 
' WH  iMR' 


) (7-1 


The  consequence  of  cyclic,  first- harmonic  downwash  was  explored  in 
Reference  13.  Their  conclusion,  which  parallels  Lockheed's  exper- 
ience, is  that  the  phase  and  magnitude  of  the  flap  response  of  a hinge- 
less blade  to  cyclic  feathering  is  markedly  affected  by  cyclic  down- 
wash.  The  shaft  moments  variation  with  feathering  angle  and  the 
phase  angle  between  flap  and  feathering  are  both  reduced  with  cyclic 
downwash,  the  effect  being  greater  in  hover  than  in  forward  flight. 


A physical  interpretation  can  be  rationalized  for  the  formula  above, 
at  least  in  hover,  in  that  the  aerodynamic  thrust  and  moment  produces 
a flow  of  linear  and  angular  momentum.  Imagine  the  flow  as  a contin- 
uous stack  of  disks  having  a mass  per  unit  thickness  pirR^  and  dimetral 
inertia  per  unit  thickness  p(irR^/U).  2 w^r,  2 p^MR  and  2 qp^R  are 
the  final,  for  downstream  position,  values  of  induced  velocities  ob- 
tain by  these  disks  orientated  with  the  flow.  The  terms  pirR^v^  and 
PttrVL  are  the  mass  flow  per  unit  time,  and  the  moment  of  inertia 


flow  per  unit  time  through  the  "actuator"  disk,  which  times  2w. 


^iMR  or  ^iMR  nioraentum. 


For  programming  purposes,  an  empirical  blend  of  the  forward  flight  and 
hovering  sets  of  expressions  is  used.  The  limiting  cases  of  the 
empirical  set  give  the  derived  cases.  The  expressions  used  are: 


-F  = pirR  V 2 w 

"Vh  a®  iMR 


(7-13) 
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(7-15) 


7. 2. 2. 3 Variations  in  Forward  Flight  and  In  Ground  Effects 


The  previous  development  can  be  assembled  and  combined  with  linearized 
forward  flight  distribution  and  ground  effect  factors. 


^BLE^  " WiMR  fiMR  1 + KiMR  R cos^r  + ^RLn  + V 

DW,BLn  L 


+ rpiMR  Sin(*R  * W + r“iMR  COs(*R  + W I7’16* 


The  ground  effect  factor,  fj_MR>  and  the  longitudinal  linear  gradient 
factor,  KiMR,  accounts  for  forward  flight. 


Note  this  formula  is  in  rotating  coordinates,  and  that  the  forward 
flight  distribution  actually  is  applied  along  the  line  of  the  apparent 
airflow,  ik,.  In  doing  this  the  distribution  is  valid  for  forward 
flight,  sideward  flight  and  sideslip  conditions.  The  angle  is 


4'..  = Tan 
W V ^ 


(7-17) 


as  shown  in  Figure  5-^* 


The  aerodynamic  moment  factors,  and  remain  attached  to  the 

hub  axis. 


The  downwash  factor  K-yvffi,  as  explained  in  assumption  3,  is  given  as 
a function  of  the  wake  angle  defined  as 


X iMR  Tan 


-1 


2 2 
H VH 


w.,,„  - wIT 
iMR  H 
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which  is  zero  in  hover  and  near  90  degrees  in  high-speed  flight.  The 
function  can  be  constrained  by  a number  of  factors.  In  hover,  the 
value  is  zero.  A 90- degree  value  of  about  1.6  can  be  read  from  Fig- 
ure 7-2.  Also  from  this  figure  a set  of  linearized  distributions  is 
read,  and  plotted  as  Figure  7-1*. 


The  ground  effect  factor. 


f 

iMR 


= 1 


-Ml)1 


(7-19) 
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is  taken  from  Reference  lU.  Where  h = -(Z0,,),,  from  Section  5. 5-1. 
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7. 2. 2.1+  Downwash  Transients 


II 

f 
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The  downwash  velocities  are  assumed  subject  to  simple  time  lags: 


wiMR  ^ wiMR , N “ wiMR^TiMR 


(7-20) 


and  similarly  for  and  q^j^.  The  value  w^j^  is  the  new  value 

obtained  from  the  steady-state  formulation  of  the  previous  sections. 

The  time  lag  is  difficult  to  estimate.  It  probably  should  be 

related  to  a length  such  that  iiMR  = JiiMR  ViMR.  Suppose  there  is  a 
sudden  change  in  collective  or  cyclic.  The  characteristic  length 
at  which  the  old  trailing  vortices  are  having  a minor  effect  on  the 
rotor  would  be  expected  to  be  of  the  order  of  a small  fraction, 
roughly  1/3,  of  the  rotor  diameter. 

Also,  de  facto,  the  time  lag  serves  as  an  averager  which  has  some 
effect  in  assisting  the  convergence  of  the  solution. 

7.2. 2. 5 Iteration  of  Downwash  Solution 

As  is  the  case  with  any  rotary-wing  loading  calculation,  there  is  an 
interplay  between  the  downwash  variance  from  calculating  the  loading 
and  a variance  in  the  loading  from  recomputing  the  downwash.  A common 
practice  is  to  solve  an  iterative  loop  to  satisfy  both  equations  (i.e., 
lift  and  momentum) . In  REXOR  the  iteration  does  not  take  place  in- 
dependently, but  proceeds  stepwise  with  the  rotor  azimuthal  advance. 
With  the  normal,  rapid  convergence  of  the  iteration  the  solution  will 
essentially  be  complete  with  the  step  advance.  However,  large  step 
sizes  will  incur  an  additional  downwash  time  lag. 

7.2.3  Blade  Element  Velocity  Components 

In  the  following  subsections  the  blade  aerodynamic  loading  is  categorized 
and  developed  along  two  lines.  They  are: 

• Steady-state  aerodynamics 

• Transient  phenomena  consisting  of  quasi-steady  aerodynamics  and 
dynamic  stall. 

7. 2. 3.1  Sources  and  Resolution  from  Blade  Motion 


The  steady  aerodynamics  are  based  on  the  air  velocities  while  the 
quasi-steady  aerodynamics  (from  flutter  theory)  and  dynamic  stall  de- 
pend on  accelerations. 
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The  air  velocity  is  the  blade  mechanical  velocities  summed  with  a com- 
ponent due  to  downwash.  In  a similar  manner,  the  air  acceleration  is 
taken  to  be  the  mechanical  blade  accelerations  minus  the  downwash  accel- 
erations. The  downwash  formulation  as  developed  in  Section  7.2.2  allows 
for  lags,  and  it  is  these  lags  that  result  in  downwash  acceleration 
terms . 

7- 2. 3. 2 Steady  Aerodynamics 

The  air  velocities  relative  to  a blade  section  are  desired  for  an  axis 
system  with  origin  at  the  quarter  chord  to  match  the  airfoil  table 
data.  From  Section  5 *5* 5,  the  mechanical  blade  velocities  relative 
to  the  free  stream  or  earth  axes  are  available  as  {&  , Y ~ , Z l1  • 

The  desired  relative  air  velocities  at  the  quarter  chord  (orblaae  BLn 
reference  axis)  are: 


BLn  - BLE 


0 

0 

“BLEJ  DW,] 


(7-21) 


Where  the  second  vector  on  the  right  is  the  downwash  velocity  de- 
veloped in  Section  7-2.2,  and  the  third  transfers  the  velocity  from 
the  BLE  reference  point  at  the  blade  center  of  gravity  back  to  the 
quarter  chord.  The  distance  Yrr  is  positive  with  the  center  of 
gravity  ahead  of  the  quarter  chord.  For  notational  convenience. 


US  ~Al/4  c 
UC  ‘ = « *1/4  c 


(7-22) 


The  angle  of  attack  is  defined  as 


a 1/4 


sin_1  (UK  / A/  + ) 


(7-23) 


Airflow  aspects  of  quasi-steady  aerodynamic  formulation  are  developed  at 
this  point  for  convenience.  The  quasi-steady  aerodynamic  contribution  is 
conceived  as  composed  of  circulatory  and  noncirculatory  components.  The 
circulatory  components  are  taken  to  be  equivalent  to  finding  the  aerodyna- 
mic force,  and  moment  coefficients  are  based  on  an  angle  of  attack  at  the 
three-quarter  chord: 


ct  . . — ct  _ . i + 

3A  c 1/h  c 


PBLE  2 

^7 


( l-2h ) 


As  such,  the  effect  of  angular  rates  is  included  in  deriving  the 
steady  aerodynamic  coefficients.  The  formulation  above  does  not 
attempt  to  account  for  local  downwash  rotation  or  curvature  and  its 
chordwise  variation.  The  net  result  is  that  aerodynamic  coefficients 
determined  in  Section  7.2.U  are  computed  with  c • 

A number  of  quantities  used  in  the  dynamic  stall  computations.  Section 
7.2.3*^,  are  also  available  from  the  previous  mechanical  development. 
They  are  also  defined  here  for  convenience.  First,  the  angle  of  side- 
slip appears  only  in  the  dynamic  stall  formulation.  For  this  purpose 
it  is  defined  as 


A = Tan 


© 
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Also,  dynamic  stall  is  based  on  the  time  derivative  of  the  angle  of 
attack  at  the  three-quarter  chord: 


1 + (UN/UCr 


UC  UC 
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The  above  equation  requires  U^  and  U . 


3Ln  - BLE 


+ J 0 


DW,BLE 


( r-pq)Y 


-( p+qr )Y 
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The  subscript  I has  been  used  to  define  motions  to  earth  axes.  The 
additive  of  the  gravity  term  places  these  accelerations  in  a true  in- 
ertial axis  system,  not  earth  inertial  axes,  as  appropriate  for  aero- 
dynamic calculations.  See  Section  5 - 5 - 1 • Gravity  does  cause  bouyancy 
forces , but  these  can  be  ignored.  The  gravity  vector  can  be  obtained 
from  hub  values  as 


1 BLn  i 


[£] 
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By  differentiating  the  downwash  velocities  , the  downwash  accelerations 
are  obtained: 


!,BLn)  Vi' 


(Wble)  DW.l 


- w f K 
iMR  iMR  iMR 
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7. 2. 3. 3 Quasi-Steady  Aerodynamics 


Quasi-steady  aerodynamics  is  accounted  for  in  REXOR  by  incorporating 
the  terms  from  the  two-dimensional  flutter  theory  of  Theodorse'.  (ref- 
erence 15).  In  the  REXOR  analysis,  Theodorsen's  lift  deficiency  func- 
tion C(k)  is  taken  as  unity.  This  means  that  the  flutter  theory 
presently  incorporated  neglects  shed  wake  effects,  or  in  physical 
terms  does  not  account  for  the  phase  change  between  blade  element  lift 
(or  pitching  moment)  and  angle  of  attack,  due  to  shed  vorticity,  or 
the  assumption  of  quasi-steady  aerodynamics  is  expressed  by  C(k)  = 1. 

Referring  to  the  classic  text  on  aeroelasticity  by  Bisplinghof f , 
Ashley,  and  Halfman  (Reference  1 6),  the  expressions  for  lift  and 
pitching  moment  are  given  as: 


L = it  pb^[h  + U a - ba  a]  + 27ipU  bC(k)  [h  + U a + b - aj  a] 

(7-30) 
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and 


M = Tipb2  [bah  - Ub  (-i-  - a)  i - b2  (4  + a2)  d] 


+ 27TPU  b2  (a  + |)  C(k)  [h  + U«+  b (1-  a)  a]  (7-  3l) 


In  REXOR  the  blade  aerodynamics  and  quasi-steady  aerodynamics  are 
referenced  to  the  local  section  quarter-chord  properties.  This  is 
done  because  the  majority  of  available  airfoil  data  uses  this  refer- 
ence.  Note  that  the  final  aerodynamic  loads  are  translated  to  the 
local  BLE  axis  (eg  location)  for  use  in  the  equations  of  motion. 


Reviewing  the  above  expressions,  and  referencing  the  rotation  point  to 
the  quarter  caord  gives  a = -1/2.  If  we  take  C(k)  as  unity,  replace 

v,  f»C1+ru  t0ry  lm  by  {dCLR/^)>  and  substitute  c/2  for  the  semi- 
chord b,  these  equations  become 


L = 


■nPc 


[Ktui*va] 
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and: 


M = 
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Note  that  the  entire  last  term  in  the  moment  equation  vanishes  with 
a - -1/2.  Referring  to  the  lift  expression,  noncirculatory  aerody- 
namic lift  is  accounted  for  in  REXOR  by  the  first  term  in  which  h + Ua 
are  combined  into  UN  in  blade  element  coordinates.  The  second  term 
results  from  table  lookup  where 


dCT 


PU2  c C 
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in  which  the  angle  of  attack  is  previously  computed  from 


‘3/Uc 


_ j.  u.  c°ti 

Lu  “ 2uJ 


The  a within  the  brackets  is  identified  as  0,  the  actual  physical 
angle  of  the  blade  with  respect  to  the  freestream  direction.  The  a 
on  the  left  hand  is  that  due  to  the  air  velocities  which  include  the 
plunging  velocity  h and  rotation  component  c/2a.  Hence  Pblic  = <*• 


The  total  aerodynamic  pitching  moment  is  the  sum  of  the  quasi-steady 
loads  computed  above  and  the  table  lookup  blade  section  properties 
(Section  7.2.4). 


7. 2. 3. 4 Dynamic  Stall 


Dynamic  stall  is  included  in  REXOR  based  upon  the  Boeing- Vertol  formu- 
lation set  forth  in  References  17,  18,  and  19-  It  is  similar  to  the 
treatment  of  dynamic  stall  in  the  Bell  C-8l  program.  A comparison  of 
REXOR  with  the  C-8l  program  is  given  in  Appendix  IV,  pages  393-404,  of 
Reference  4.  Dynamic  stall  is  specifically  addressed  with  respect  to 
the  two  programs  beginning  on  page  395  of  that  report.  A significant 
point  of  difference  between  the  treatment  of  dynamic  stall  in  the  two 
programs  is  that  C-8l  puts  a 20-percent  limit  on  the  angle-of-attack 
overshoot  in  obtaining  the  dynamic  maximum  lift  coefficient,  whereas 
REXOR  has  no  limit.  The  correctness  of  the  treatment  of  dynamic  stall 
in  either  program  is  difficult  to  assess  since  the  concensus  of  re- 
searchers in  this  area  is  that  current  methods  are  empirical  at  best, 
and  much  research  still  remains  to  be  done  in  this  area. 


Reference  17  notes  that,  "The  trends  show  that  compressibility  effects 
reduce  dynamic-stall  delay,  and  at  about  M = 0.6  no  dynamic-stall 
delay  is  evident."  For  this  reason  an  upper  Mach  number  limit  of  0.6 
was  implemented  in  the  dynamic  stall  calculations  for  REXOR.  The  test 
data  obtained  by  Boeing  Vertol  and  given  in  the  references  cited  was  for 
for  the  Mach  number  range  0.2  to  0.6.  As  implemented  in  REXOR  if 
M < 0.25,  the  value  M = 0.25  is  used  in  the  analytic  expression  for  de- 
veloping the  stall  hysteresis  loop. 


Reference  18  notes  that  it  was  found  that,  "airfoils  used  currently 
by  the  helicopter  industry  had  stalling  dominated  by  leading  edge 
stall. 


For  this  type  of  stalling  process,  the  dynamic  extension 
was  proportional  to  the  time  rate  of  change  of  the  angle  of  attack." 
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In  that  reference,  so  as  to  use  static  airfoil  data  as  much  as 
possible,  static  stall  and  dynamic  stall  are  empirically  related  by 
developing  a reference  angle  of  attack  given  by 


“ref  “ ~ 


sign(c 


(7-36) 


in  which. 


Y = log. 


0.601 

M 


(7-37) 


and  is  physically  related  to  dynamic  stall  delay,  a is  identified  as 


M = a 


VUN 


2 2 
+ UC 


(7-38) 


As  noted  in  Reference  19  in  regard  to  dynamic  stall..."  as  a blade 
element  reaches  and  exceeds  the  static  angle  of  attack,  stall  does  not 
occ  jt  as  long  as  a sufficient,  positive  time  rate  of  change  of  the 
airfoil  angle  of  attack,  a,  is  present."  The  experimentally  derived 
equation  for  dynamic  stall  delay  is  given  in  the  reference  as 


dynamic  stall  delay  = y 


(7-39) 


where 


(7-i+O) 


the  blade  element  reduced  frequency. 

Referring  to  the  gamma  expression.,  we  note  that  Y + 0 as  M + 0.601, 
which  is  the  upper  limit  for  Mach  number  values  for  'dynamic  stall  cal- 
culations. Also,  note  that  Y ->-1  as  M + 0.2211,  which  is  approxi- 
mated by  the  value  of  M = 0.25,  the  lower  limit  in  REXOR  for  dynamic 
stall  simulation. 


The  termapgp  given  above  is  also  called  the  dynamic  angle  of  attack 
(Reference  19)  and  given  by  the  notation  odYN* 


7.2.3. 4.1  Lift  Accounting  for  Dynamic  Stall 


Using  the  reference  or  dynamic  angle  of  attack  computed  fromct  , the 
REXOR  program  implements  the  "Fast  Aerodynamic  Table",  Section^1. 2. 4, 


subroutine  and  determines  the  lift  coefficient,  CL  , corresponding  to 
«r,r-n  and  the  freestream  Mach  number  for  the  specified  blade  element 


Oflgjp  and  the  freestream  Mach  number  for  the  specified  blade  element 
and  blade  azimuth  position.  Also  computed  at  the  given  Mach  number 


are  the  (L  for  zero  angle  of  attack  and  the  CL  for  a small  incre- 
ment A a with  respect  to  zero.  Yawed  or  radial  flow  is  accounted  for 


by  computing  the  yaw  angle  of  the  flow  given  by: 


A = Tan" 


© 
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where  Ug  and  represent  blade  spanwise  and  chordwise  components  of 
flow  respectively. 


The  slope  of  the  lift  curve  is  then  found  from: 


fe) 


CL(aREF,M)  - CL(°’M) 


a REF  cos  A 
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It  can  be  argued  from  physical  reasonings  that  the  dynamic  lift- 
curve  slope  cannot  exceed  the  static  life-curve  slope.  As  a check, 
REXOR  also  calculates: 


© 


CL(Aa,  M)  - Cl(0,M) 


(7-43) 


Only  in  the  event  (3Cl/3oi)dyN  is  greater  than  (3CL/3a)0>M  is 
latter  value  used  to  calculate  C^.  Otherwise  Cl  is  calculated  by 


■© 


a + Cl(0,M) 


(7-44) 
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The  ability  of  this  approximation  to  describe  mathematically  the  lift 
hysteresis  characterized  by  dynamic  stall  is  shown  in  Figure  7-5, 
which  compares  analytical  results  with  experimental  two-dimensional 
airfoil  data.  (From  Reference  19). 


The  component  of  the  lift  force  per  unit  span  acting  normal  to  the 
blade  chord  axis  and  including  dynamic  stall  ef  bets  is  then  cal- 
culated from 


AFNC  = CL  c — c°sa 


(7->+5) 


The  total  normal  force  is  determined  by  adding  to  this  term  the  drag 
component,  Cp  c PV2/2  sin  a,  and  the  unsteady  aerodynamic  terms  dis- 
cussed in  the  previous  section.  To  account  for  dynamic  stall  effects 
on  drag,  two-dimensional  drag  coefficient  data  are  used,  but  as  de- 


termined at  a 


not  a.  This  is  consistent  with  Reference  18. 


The  component  of  the  lift  force  per  unit  span  parallel  to  the  blade 
chord  axis  is  found  correspondingly  from: 


CL  C PV 

AFC  = sin  a 


(7-1+6) 


The  total  chordwise  force  is  then  obtained  by  adding  the  correspond- 
ing drag  coefficient  term  multiplied  by  cos  a . 

7.2.3.1+.2  Pitching  Moment  Accounting  for  Dyanmic  Stall 

For  determining  pitching  moments  due  to  dynamic  stall  (see  Reference 
17),  the  reference  or  dynamic  angle  of  attack  given  by  otR-n,„  must  be 
modified.  In  REXOR,  this  is  accomplished  by  multiplying  une  second 
term  by  an  empirical  constant,  K.  Hence, 


= a _ K 


—■  sign(a)^ 
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K is  selected  baaed  upon  the  dynamic  stall  characteristics  of  the  air- 
foil. In  general  it  has  been  found  for  conventional  rotor  blade  air- 
foils that  K should  be  selected  so  that 


(7-48) 


ct  • = a + A a 

REF  REF 


where  Aa  is  of  the  order  of  2.5  degrees.  With  “p^F  calculated  from 
the  above  equation,  the  moment  coefficient  is  determined  from  tables 
such  that, 


CM  " CM^°REF’ 


(7-49) 


A comparison  of  test  and  theoretical  dynamic  from  Reference  19  is 
shown  in  Figure  7-6. 

The  total  pitching  moment  acting  per  unit  span  on  a blade  element  is 
then  given  by: 


T(i)  = - 


CM' 


2 pV 


- PN0  SY(1> 


’quasi-steady' 
aero 
terms 

(Section  1.2.3. 3) 
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where  S,,(x)  represents  the  distance  from  the  aerodyanmic  center  to 
the  blade  elastic  axis,  and  the  quasi-steady  aerodynamic  terms  are 
included  as  described  in  Section  7- 2. 3. 3. 

7*2.4  Coefficient  Table  lookup  - Overview 

In  cataloging  blade  section  aerodyanmic  data,  CL,  CQ  and  CM,  there  are  two 
procedures  available. 

• Curve  fit  the  aerodynamic  data  to  the  specific  airfoil  goemetry 
being  investigated  for  the  range  of  Mach  number  and  angle  of 
attack  to  be  considered. 

• Tabulate  the  data  as  a function  of  performance  and  geometric  param- 
eters, and  interpolate  to  the  exact  conditions  at  hand. 

REXOR  uses  the  second  procedure.  The  data  consist  of  a series  of  synthe- 
sised airfoil  data  for  five  digit  series  airfoils.  The  synthesising  pro- 
cess takes  measured  rotor  data,  and  by  means  of  numerical  integration  per- 
formance formulae,  determines  the  equivalent  two-dimensional  airfoil 
section  properties  required. 

Two-dimensional  airfoil  data  from  Lockheed  tests  and  technical  literature 
were  used  at  a starting  point  to  synthesize  the  data.  Data  for  thickness 
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Figure  7-6.  Dynamic  Stull  - Moment  Coefficient  vs 
Angle-of-Attack  Hysteresis  Loop 


ratios  and  camber  combinations,  not  previously  available,  were  developed 
by  using  transonic  similar  rules  as  indicated  in  Reference  20.  These 
data  were  cross  checked  for  consistency  with  the  data  of  the  many  sources 
noted  in  the  bibliography  of  that  reference.  The  rotor  airfoil  data  syn- 
thesizing process  from  two-dimensional  and  data  included  the  tip 
Mach  relief  and  C and  CL  adjustments  for  rotor  applications  in  a manner 
consistent  with  the  synthesized  rotor  airfoil  data  of  Reference  21. 

The  CL,  Cj)  data  covers  uncambered  thicknesses  ratios  from  6 to  12  percent 
and  design  lift  coefficients  of  0.09  to  0.69  for  cambered  sections  of  the 
same  range  of  thickness  ratio.  For  angles  of  attack  within  j\30  degrees 
the  data  is  also  presented  as  a function  of  Mach  number  to  account  for 
compressibility.  For  large  angles  (above  +30  degrees)  compressibility 
effects  are  ignored.  Blade  element  pitching  moment  data  is  also  arranged 
along  these  lines,  and  is  available  for  an  uncambered  12  percent  thickness 
ratio  section  and  a 2~percent  camber,  8-percent-thick  section. 

7. 2.4.1  Inputs  and  Outputs 

Two  alternate  methods  are  available  in  REXOR  to  use  the  de  a bank 
described  above.  First  a lengthy,  but  generalized  procedure,  uses 
the  data  bank  directly.  From  the  spanwise  design  lift  coefficient 
variation,  the  existing  span  station  angle  of  attack  and  Mach  number 
the  tables  are  multiple  interpolated  to  yield  a and  CD.  The  CM 
coefficient  is  selected  from  one  of  the  two  sections  available. 

The  second  and  the  most  common  method  of  table  usage,  is  to  pre- 
interpolate the  base  tables  described  into  a set  keyed  to  uhe  re- 
quired geometry  and  arranged  for  rapid  access.  Currently  the  data 
take  advantage  of  blade  geometry  characterized  by  either  being  con- 
stant thickness,  camber  or  possessing  a linear  relationship  between 
span  - thickness  and  span  - camber  to  specify  only  thickness  ratio 
and  automatically  specify  the  other  ingredients.  The  data  have  been 
curve  fitted  for  arbitrary  variations  in  camber  again  specified  as  a 
function  of  thickness  ratio.  The  data  sets  referred  to  in  this  sec- 
ond manner  are  known  as  the  fast  aero  tables. 

Note  that  by  proper  substitution,  other  tables  and  lookup  procedures 
may  be  employed  in  REXOR. 

7. 2. 4. 2 Organization  of  Data  in  Tables 

The  base  data  bank  contains  a total  of  31  tables.  These  are  primarily 
referenced  by  design  lift  coefficient  (and  coincidently  camber)  and 
thickness  ratio.  Secondary  sorting  is  according  to  angle  of  attack 
range  and  if  applicable,  Mach  number. 

The  fast  aero  tables  are  a set  of  six  tables,  two  each  for  lift,  drag 
and  pitching  moment  coefficients.  The  double  tables  are  for  angles 
of  attack  below  +30  degrees  (also  a function  of  Mach  number)  and  above 
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±30  degrees.  These  tables  are  generated  for  a specific  blade  con- 
figuration and  referenced  by  thickness  ratio. 

Detailed  contents  of  the  tables  are  covered  in  Section  U of  Volume  II. 
7. 2. It. 3 Accession  Method 


As  mentioned,  the  base  data  set  is  referenced  to  the  desired  Cp,  Cp  and 
Cm  by  a multiple  interpolation  scheme.  The  fast  aero  tables  require  a 
single  tri variant  interpolation  (with  Mach  effects)  to  arrive  at  the 
required  coefficient.  By  making  the  table  entries  at  equally  spaced 
intervals , an  entry  index  is  simply  thn  entry  value  divided  by  the 
spacing  value.  The  trivariant  interpolation  scheme  is  shown  in  Figure  7-7. 


7-2.5  Blade  Element  and  Rotor  Aerodynamic  Loads  Summary 

The  required  loads  for  use  in  the  equations  of  motion  are  in  B£n  axis. 
Development  to  this  form  from  BLE  axis  about  this  quarter  chord  point  is 
covered  in  Section  6.6.U.  The  BLE  axis  form  is: 
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where , 


[’’“-bleJ 


0 COS(a3A  c)  -sin(“3/l.  c) 

0 sin(“3A  o)  cos  (“3/4  c) 


(7-52) 


7.3  INTERFERENCE  TERMS 


7*3.1  Nature  of  the  Phenomenon 


In  the  process  of  producing  lift,  the  various  parts  of  the  rotorcraft  im- 
part a net  momentum  change  to  the  air  mass  opposite  to  the  direction  of  the 
force  produced.  This  induced  air  velocity  from  the  momentum  change  impinges 
upon  other  elements  of  the  rotorcraft  changing  their  aerodynamic  behavoir. 


The  sources  of  interest  are  the  main  rotor  and  wing  (or  lifting  body  charac- 
teristics of  the  fuselage).  The  surfaces  being  affected  are  the  wing  plus 
fuselage  and  horizontal  tail.  The  impinging  velocity  is  expressed  in  Zp 
(fuselage  verticle)  axis  as  a percentage  of  the  source  flow  and  a function 
of  the  wake  angle  of  this  flow. 


A second  interference  velocity  source  is  to  consider  the  circulation  part 
of  the  Theodors en  function.  Here  the  wing  or  wing  equivalent  of  the  fuse- 
lage is  producing  lift  at  the  quarter-chord  point  according  to  the  air 
velocity  at  the  3A-chord  location.  Accordingly,  the  vertical  component 
of  air  velocity  at  the  wing  includes  a component , 


-i.  c q 

2 wing  yF 


(7-53) 


Here  the  wing  quarter  chord  is  assumed  to  lie  on  the  Yp  axis.  This  com- 
ponent is  also  effective  at  the  horizontal  tail  via  the  wing  to  horizontal 
tail  downwash  factor. 


7.3.2  Rotor  to  Wing/Fuselage 


The  downwash  function  (percentage  of  source  flow)  used  in  REXOR  is  a lookup 
table  of  downwash  factor,  FY  , and  idealized  main  rotor  wake  angle  X 

MR-2 


where , 


XMH  “ ta”' 


-1  ^ VH  ~ WiMR  ^ 


(7r5M 


The  table  data  is  linearly  interpolated  to  the  required  wake  angle  value. 
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Figure  7-7*  Trivariant  Linear  Interpolation 
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The  fuselage  reference  downwash  velocity  at  the  wing  (or  equivalent)  then 
is 


w =v_v  f + — C a 
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(7-55) 
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and  taking  time  derivatives, 
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The  total  air  velocity  to  the  wing/fuselage  is 


WING 


/ 2 2 2\ 
= lu  + v + w 1 
y F F WING  ) 


1/2 


(7-57) 


and  the  angle  of  attack  is 


WING 


(7-58) 


The  total  velocity  in  the  fuselage  XZ  plane  is  used  in  the  horizontal  tail 
computations : 


V. 


X„ 


= (u  2 + v 2\ 
^ F WING  / 


1/2 
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7 • 3 - 3 Rotor  to  Horizontal  Tail 

A downwash  factor  F-^^  ^rp  between  the  main  rotor  and  horizontal  tail  is 

computed  in  the  same  manner  as  F^^  ^ frori  the  main  rotor  wake  angle  X^. 

This  data  in  conjunction  with  the  wing  to  horizontal  tail  downwash  factor 
is  used  to  compute  incremental  air  velocities  at  the  horizontal  tail.  The 
fuselage  lift,  drag  and  moment  data,  Section  7.*+,  is  assumed  to  be  complete, 
less  the  downwash  factors.  Hence,  only  incremental  correction  velocities, 
which  in  turn  will  be  used  in  computing  load  increments,  are  needed. 
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Evaluating  the  main  rotor  increment, 


"imi  = “1MR  ( Vht  ' Vw) 


(7-6o) 


7.3.*+  Wing  to  Horizontal  Tail 


The  wing  to  horizontal  tail  downwash  factor  appears  explicitly  as  a quasi- 
unsteady aerodynamic  term.  The  steady  term  is  accounted  in  an  overall 
fuselage  aerodynamic  characteristic  set.  An  airflow  time  delay  from  the 
wing  to  horizontal  tail  is  computed  as 


At  = 


(7-61) 


Using  the  downwash  factor  3 e/a  a the  vertical  airflow  component  at  the 
horizontal  tail  is 
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The  total  incremental  velocity  components  at  the  horizontal  tail  are 


Aw  = - w 
HT  iMR, 


q + JJLw  /y 

^F  Set  WING7  Xn 


(7-63) 


AVHT  " £VT  rF  + hVT  PF 


No  fuselage  - wing  induced  side  wash  is  modeled  in  REXOR. 


7-3.5  Data  Sources 


The  theoretical  downwash  factor  ranges  from  0 at  X = 0,  180  degrees  to  2 
at  X = 90  in  the  fully  contracted  rotor  wake.  Several  sources  of  measured 
data  are  available  to  construct  a distribution  for  a given  configuration. 
Reference  11  gives  isolated  rotor  data  for  field  distances  and  wake  angle 


Reference  22  gives  a good  data  set 


ranges  suitable  for  FXmr_w  and 
for  typical  wing  locations. 


7-4  BODY  LOADS 


7.4.1  Nonrotating  Airframe  Airloads 


The  required  loads  are  computed  in  REXOR  as  the  sum  of  steady- state  forces 
and  moments  plus  loads  arising  from  stability  derivative  type  terms.  The 


steady- state  data  are  formed  in  terms  of  overall  C^,  C^,  and  C for  the 
fuselage,  wing,  and  empennage  assembly.  Typically  these  data  are  from 


wind  tunnel  tests. 


The  static  body  loads  are: 
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where : 


1 2 
0 = — P S V 

1A.  2 WING  X„ 

Li 


(7-65) 


The  wing  area,  S^^_,  and  chord,  are  actual  or  the  equivalent  of  the 

lifting  fuselage.  Alternately,  they  may  be  the  reference  length  and  area 


used  for  the  available  wind  tunnel  data.  Cp  , Clt  and  are  linearly  in- 


terpolated  from  input  data  tables  of  C^,  C^,  C^,  versus  angle  of  attack. 

The  data  are  interpolated  on  a from  Section  7.3-2.  The  loads  developed 
are  in  wind  axis. 


The  stability  derivative  load  contributions  are  computed  as  a 6 by  6 de- 
rivative matrix  postmultiplied  by  a velocity  component  vector. 
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where, 


Fm(1,  2)  = ~ PSWING  CD  V « Drag  due  to  side  slip  squared. 


(7-67) 


fmn(2’  3)  = 2 PWTI  = side  force  due  to  sidfc  slip‘  (7"68) 
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- Vertical  force  due  to  vertical  velocity 
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Roll  moment  due  to  wing  incidence 
differential. 
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~ Wing  roll  damping.  (7-73) 
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W1 * * * 5’  2)  =T  pSWING  CWING  CM  v = Pitchin®  moment  due  to 

0 side  slip. 
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FMN(6’  3)  =_2  PSWING  bWING  “dJ  = YawinS  moment  due  t0  Side  slip- 

(7-77) 

1 dCN 

W6-  61  =2  PSWI»G  = ^,ing  moment  due  to 

lateral  velocity.  (7-78) 


These  terms  also  produce  forces  and  moments  in  wind  axes. 

The  static  and  derivative  terms  are  added  to  form  the  total  body  loads  and 

transformed  into  fuselage  axes. 
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and 
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The  air  velocities  wWING,  V V , are  defined  in  Section  7.3.2. 

Y WING 


7-^.2  Component  Additional  Airloads 


A total  array,  {QLOADs},  of  non  main  rotor  air  loads  is  computed  in  fuse- 
lage axes. 
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The  first  component  is  described  above.  The  tail  rotor  load  vector  and 
propeller  load  vector  are  developed  in  the  following  sections. 


7.5  TAIL  ROTOR 

A number  of  different  levels  of  aerodynamic  presentation  accuracy  and  axis 
of  representation  may  be  used  for  tail  rotor  computations.  In  line  with 
the  stated  objectives  of  REXOR,  a linear  aerodynamic  approach  is  used.  A 
shaft  axis  reference  is  used  for  the  analysis.  In  this  system,  the  air 
velocity  quantities  involved  are  easy  to  visualize.  Also,  the  flapping  and 
feathering  motions  are  the  true,  measureable  quantities. 

7.5.1  Formulations 

First,  consider  the  airflow  quantities  available  in  fuselage  axis.  Fig- 
ure 7-8.  Note  the  tail  rotor  axes  align  with  the  fuselage  axis  system. 

The  detailed  velocity  components  are  given  in  Figure  7-9.  Constructing  the 
blade  element  tangential  (U^)  and  perpendicular  (Up)  components, 


UT  = sin*TE 


(7-87) 


UP  “ VF  + Sr  ^F  WiTR  r8  " uf8  cos  ^TR 


(7-88) 


where, 


w is  the  tail  rotor  induced  velocity 
iTR 


V = _ v + £ j)_w 
TR  F TR  VF  iTR 
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Expressing  the  blade  element  angle  of  attack  as  a small  angle  of 
approximation, 


aTR  9+Ut 
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where, 


6=  STH  - A1  C0S  +TR  - B1  Sin*TR 


(7-91) 


A { coupling  is  used  to  minimize  tail  rotor  flapping.  Defining  + 6 _ 

"tR  ^TR 

as  a reduction  in  feathering  for  positive  flapping  gives 
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The  tail  rotor  analysis  assumes  no  coning. 
The  blade  flapping,  3,  is  then 


0 = - a cos  iji  - b sin  ♦ 
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The  tail  rotor  expressions  of  interest  are  the  prime  forces  added  to  the 
fuselage  system.  First  looking  at  the  tail  rotor  thrust.  For  a blade  ele- 
ment we  have 
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where  a is  the  lift  curve  slope,  b is  the  number  of  tail  rotor  blades,  and 
c is  the  blade  chord  (assumed  constant). 


Substituting, 


d TTR  = [i  Pa  b c (0UT2  + Up  UT)dr] 


(7-95) 
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Integrating  for  the  entire  rotor, 
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where  B is  the  finite  airfoil  lift  factor  expressed  as  a so-called  tip 
loss  factor. 

Noting  only  even  functions  contribute  to  the  integrand. 


TTR  2 pa^b  C RVr  (0TR  3 ^R^TR  + 0TR  2 UF 


+ 2 (fiB)TR  VTR  + (fiR)^  up  b -g-  63tr^ 
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Note  the  thrust  is  independent  of  the  longitudinal  flapping,  but  is  a func- 
tion of  lateral  cyclic  shown  as  lateral  flapping  times  delta  3. 

The  required  lateral  flapping  angle  is  obtained  by  equating  the  lateral 
flapping  moment  equal  to  zero. 
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gives 
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To  obtain  the  longitudinal  flapping  angle,  the  longitudinal  rotor  moment  is 
formed  and  set  equal  to  zero. 
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In  formulating  the  tail  rotor  drive  torque,  the  blade  profile  drag  is  ex- 
pressed as 


CD  = CDq  + k V 


(7-102) 


where  is  the  average  lift  coefficient.  Reviewing  the  thrust  equation 
with  a constant  (average)  lift  coefficient  gives 
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The  drive  torque  is  expressed  as  the  reaction  to  turning  the  tail  rotor 
shaft.  The  sense  of  rotation  is  clockwise,  facing  a left-hand  mounted  tail 
rotor. 


■H 


2 — ]_  2 

p c b r UT  CD  + p a b c a <j>  r 


(7-10U) 
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integrating 
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The  remaining  load  term  in  X^p.  Using  the  same  formulation  methodology , 
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Making  small  angle  approximations, 


d XTR  “ 2 P^b  C^TR 


integrating 
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The  induced  velocity  is  calculated  from  simplified  momentum  balance. 
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Normally,  the  thrust,  flapping,  and  induced  velocity  equations  are  solved 
as  an  iterative  set.  In  REXOR,  these  equations  are  solved  for  every  pass 
(azimuth  step)  of  the  main  rotor,  and  the  tail  rotor  se"C  convergence  is 
assumed  a priori. 


Note  that  the  pitch-flap  coupling  does  not  appear  in  the  expressions  de- 
veloped. This  is  due  to  the  equivalence  of  flapping  and  feathering, 
coupled  with  the  absence  of  lateral  flapping. 


7*5-2  Airloads  - Control  Settings 


The  force  and  moment  terms  are  assembled  for  use  in  the  overall  body  loads. 
Section  7*^.  The  pilot  control  is  the  rudder  pedals  toe 
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7.6  PROPELLER 


7*6.1  Formulations 


A pusher  propeller  is  modeled  in  REXOR,  and  is  referenced  to  the  fuselage 
axis  as  shown  in  Figure  7-10.  The  forces  and  moments  generated  by  the 
propeller  are  added  to  the  |F^|  vector  of  total  fuselage  forces  and 
moments.  These  components  are  transformed  to  the  principal  reference  set 
for  use  in  the  equations  of  motion. 
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The  propeller  aerodynamics  are  not  calculated  from  blade  theory,  but  in- 
stead are  interpolated  from  performance  date,  iu  the  form  of  nondimensional 
thrust  and  power  curves.  See  Figure  7-11*  These  curves  are  in  terms  of  a 
thrust  and  power  coefficient  scheme  unique  to  propeller  work.  Defining 
these  coefficients: 


J = Advance  Ravio  = 


NPDP 
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= Thrust  coefficient  = 
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Cp  = Power  coefficient  = 


Qp 


pnp2  dp5 


(7-114) 


where : 


Np  = rotational  speed,  revolutions/second 


Dp  = propeller  diameter,  feet 


Vp  = axial  inflow  velocity,  feet/second 


P = air  density  slug/feet' 


In  REXOR,  the  fuselage  inertial  X velocity  v^  is  used  for  V.  The  rotational 
speed  is  times  the  main  rotor  to  tail  rotor  gear  ratio  less  the  gear 
(fuselage)  inertial  roll  rate. 


case 


7.6.2  Airloads  - Control  Settings 


The  coefficient  curves  are  entered  with  the  calculated  advance  ratio  and 
command  prop  blade  angle.  The  curves  are  bi variant  interpolated  to  the 
required  and  Cp.  These  quantities  are  then  dimensionalized  as  shown 


aoove . 
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LIST  OF  SYMBOLS 


The  symbols  used  in  the  REXOR  equations  are  quite  numerous.  In  order  to 
keep  the  catalog  of  symbols  to  manageable  proportions  the  following  list 
is  divided  according  to  the  discussion  in  Section  U.  Namely,  a list  of 
basic  symbols  is  given,  followed  by  subscripts,  superscripts,  and  post- 
scripts. Nonconforming  cases  of  usage  together  with  complicated  or 
obscure  subscripting  are  fully  annotated  in  the  basic  list. 


SYMBOLS 


0 


“l 

[A] 


A 


1,2,3 


In 


3n 


IS 


IS 


[C] 


"M 


arbitrary  vector 
speed  of  sound 

2 

acceleration  vector,  ft/sec 

longitudinal  component  of  blade  first  harmonic  flapping,  rad 
generalized  mass  element  matrix 
modal  variables 

generalized  displacement  of  nth  blade , first  mode 
generalized  displacement  of  nth  blade,  second  mode 
generalized  displacement  of  nth  blade,  third  mode 
cosine  component  of  blade  first  hai*monic  cyclic,  rad 
number  of  main  rotor  blades;  arbitrary  vector 
dissipation  function 

sine  component  of  blcde  first  harmonic  cyclic,  rad 

blade  segment  chord,  ft 

damping  matrix 

aerodynamic  drag  coefficient 

aerodynamic  lift  coefficient 

aerodynamic  pitching  moment  coefficient 

power  coefficient 

thrust  coefficient 
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"X, Y , Z 


>,0,* 


1,2,3 

C(k) 


d 

dr 

dt 

d/dt 

(d/e) 


0 


(d/e). 


El 


fiMR 


X,Y,Z 


4>  ,0  ,ip 


0PH 


%Y,Z 


G 

M 


G 

GJ 

I,. 


linear  damping,  ]b/ft/sec 

rotary  damping,  ft-lb/rad/sec 

blade  bending  to  feathering  couplings 

lift  deficiency  function 

infinitesimal  increment 

increment  in  rotor,  radius,  ft 

increment  in  time,  sec 

derivative  with  respect  to  time 

swashplate  to  feather  gear  ratio,  zero  collective 

swashplate  to  feather  gear  ratio  slope  with  collective 

pitch  horn  effective  crank  arm,  ft 

blade  bending  stiffness  distribution,  lb-ft2 

ground  effect  factor  for  main  rotor 

factor;  force,  lb 

force  components  along  X,Y,Z  directions,  lb 
generalized  force  about  ip , 0,  axis 


feathering  mode  generalized  force 
2 


gravity,  ft/sec‘ 

gravity  components  along  X,Y,Z  directions 
gear  ratio 

generalized  force  vector 


fey ro  angular  acceleration  partial  product 

2 


blade  torsional  stiffness , lb-ft ‘ 


Enu  X^2,  slug-ft2 


~ Zmi  ^i^»  slug- ft 
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r 


z 

= 

Zm^ 

2 2 
Z^  , slug-ft 

XX 

= 

Znu 

<v- 

p 

Z^'),  slug-ft 

YY 

Zm. 

l 

<xi2  + 

Z^2),  slug-ft 

ZZ 

= 

Zm. 

l 

(X^  + 

Yi2),  slug-ft 

XY 

= 

Zm. 

l 

X.  Y.  , 

l l 

slug-ft2 

XZ 

= 

Zm. 

l 

X.  z. , 
1 1 

slug-ft2 

YZ 

= 

Zm. 

l 

Y.  Z.  , 

l l 

slug-ft2 

J 

J 

k 


[K] 

K 


mj 

h,Y,Z 
K, 


<|>  ,0  ,i|> 


IB 


OB 


TTI 


ra 


"V 

"h 


m. 

1 


unit  vector 
unit  vector 
advance  ratio 


number  of  blade  radial  stations;  reduced  frequency, 
rad/sec ; unit  vector 


spring  matrix 

blade  spring  matrix  element 

spring  constants  along  X,Y,Z  direction,  lb/ft 
spring  rates  about  <J>,  8,  <J>  axis,  ft-lb/rad 
location  inboard  feather  bearing,  ft 
location  outboard  feather  bearing,  ft 


radial  location  of  intersection  of  precone  and  feather 
axis , ft 


tension  torsion  pack  length,  ft 

rolling  moment,  ft-lb 

mass  of  element,  slugs 

summed  fuselage  coordinate  mass,  slugs 

summed  hub  axis  mass , slugs 

mass  of  ith  particle  or  blade  segment,  slugs 
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£ 


a 


■i 


! 


mSP 

swashplate  summed  mass , 

slugs 

M 

pitching  moment,  ft-lb; 

= Em.  , 

slpgs; 

[M] 

generalized  mass  matrix 

s 

MrR 

generalized  mass  matrix 

element 

Mx 

= Em.  X. , slug-ft 

1 

“y 

= Em.  Y.,  slug-ft 
11 

tsif 

= Em.  Z.,  slug-ft 
11 

Hx.y.z 

moments  about  X,Y,Z  axis 

, ft-lb 

blade  torsional  moment, 

ft-lb/ ft 

QLOADS 


Z«J>  ,Z0 


number  of  system  particles 

angular  velocity  about  X axis,  rad/sec;  particle 

main  rotor  pitch  moment  inflow,  ft/sec 

generalized  coordinate;  angular  velocity  about  Y axis, 
rad/sec 

main  rotor  roll  moment  inflow,  ft /sec 
generalized  forcing  function 

aerodynamic  press  ire  times  reference  wing  area,  lb 
total  nonmain  rotor  aerodynamic  loads  matrix 
tail  rotor  torque,  ft-lb 

general  vector;  radius  of  curvature,  ft;  angular  velocity 
about  Z axis,  rad/sec;  notation  for  (X,Y,Z) 

static  blade  shape 

vector  displacement  of  particle  p in  X,Y,Z  axis  system 
vector  displacement  of  x,y,z  origin  in  X,Y,Z  system 
gyro  damper  coupling  ratios 


. .1*  - > - •-).  > ■■■  iWfr'i-i  rir'i  y-f  ■r.'irf  r\  r fcr/n'i'Aiiyit)  ill  l 


C,P ,S,T 


1,2,3 


a ■ W J /t-,7  ■ 


Laplace  variable,  path  of  motion  of  particle  p 


blade  spline  length  along  neutral  axis  locii,  ft 


kinetic  energy,  ft-lb 


transformation  of  coordinates  matrix 


tension  in  tension  - torsion  pack,  lb 


velocity  in  X direction,  ft/sec 


potential  energy  function,  ft-lb;  strain  energy,  ft-lb 


air  velocity  on  blade  element,  ft/sec 


velocity  in  Y direction,  ft/sec 


trajectory  velocity 


velocity  in  Z direction,  ft/sec 


main  rotor  collective  inflow,  ft/sec 


tail  rotor  collective  inflow,  ft/sec 


motion  in  X direction,  ft;  blade  span  location 


coordinate  direction;  axis;  def2ection,  ft;  location,  ft; 
cross  product 


blade  radial  station  of  sweep  and  jog,  ft 


trajectory  path,  ft 


tail  rotor  longitudinal  force,  lb 


motion  in  Y direction,  ft 


coordinate  direction;  axis;  deflection,  ft;  location,  ft 


tension  torsion  pack  outboard  end  modal  coefficients 


difference  between  Y direction  locations  of  eg  and  neutral 
axis  points  of  blade  element,  ft 


SP 


TTO 


1,2,3 


JOBL 


OF 


JOSP 


JOTTI 


motion  in  Z direction 

coordinate  direction;  axis;  deflection,  ft;  location,  ft 

relative  swashplate  vertical  displacement  with  respect  to 
the  hub , ft 

tension -torsion  pack  outboard  end  modal  coefficients 

teetering  rotor  undersling,  ft 

hub  set  distance  above  fuselage  set,  ft 

hub  set  distance  above  swashplate  set,  ft 

blade  vertical  offset  at  outboard  end  of  tension  - torsion 
pack,  ft 


angle  of  attack,  rad 


a 


2 


6 

B 

B 

B 


FA 

PHn 

0 


Y 

YT 

6 

63TR 

9e/'9a 


Q 

9 

e 

A 

p 


angle  of  attack  with  hub  set,  rad 

sideslip  angle,  rad 

blade  feathering  angle,  rad 

feathering/pitch-horn  bending  or  dynamic  torsion 
generalized  coordinate  displacement 

blade  droop  relative  to  precone  angle,  rad 

blade  sweep  angle,  rad;  dynamic  stall  delay,  sec 

trajectory  path  angle  with  E set,  rad 

limit  deflection,  rad;  freeplay,  rad;  small  increment 

tail  rotor  pitch  - flap  coupling 

downwash  factor  of  wing  on  horizontal  tail 

vector  notation  of  <J>,  0,  ip 

rotation  about  Y axis,  rad 

collective  blade  angle,  rad 

sideslip  at  blade  element,  rad 

3 

air  density,  slugs/ft 
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SUBSCRIPTS 


time  constant,  sec;  natural  period,  sec 


feathering  axis  precone,  rad 


rotation  about  X axis,  rad 


feathering  angle , rad 


feathering  angle  of  blade  element  of  nth  blade,  rad 


blade  root  reference  feather  angle,  rad 


blade  torsion,  rad 


sum  of  blade  twist  and  torsion,  rad 


wake  angle  of  main  rotor,  deg 


rotation  about  Z axis,  rad;  sideslip  angle  with  hub  set,  rad 


control  input  axis  rotation  from  swashplate , rad 


pitch  lead  angle,  deg 


trajectory  path  yaw  with  E set,  rad 


main  rotor  apparent  airflow  angle,  rad 


rotational  speed,  rad/secj  angular  velocity,  rad/sec; 
natural  frequency,  rad/sec 


partial  derivative,  derivation 


arbitrary  coordinate  set  a 


due  to  aerodynamics 


arbitrary  coordinate  set  b 


associated  with  blade  elastic  bending 


blade  element  coordinate  system 


blade  reference  axis  system  for  the  nth  blade 
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}■ 

i 

I 

I 


I 


c 

CG 

CORR 

DW 

DYN 

E 

ENG 

EST 

F 

FA 

FB 

Fn 

FR 

G 

GEN 

GFB 

GSP 

GUB 

H 

HT 

i 

IB 

J 

jog 


associated  with  pilot  control  input,  chordwise 

associated  with  center  of  gravity  location 

corrective,  correction 

referring  to  downwash 

referring  to  dynamic  component 

earth  axis 

associated  with  powerplant  - engine 
estimated 

fuselage  axis;  associated  with  blade  feathering 
referring  to  blade  feather  axis 
associated  with  feedback 

associated  with  feathering  of  the  nth  blade 
due  to  friction 

referring  to  gyro  or  gyro  coordinate  system 
associated  with  gas  generator  section  of  powerplant 
associated  with  gyro  control  feedback 
gyro  to  swashplate  connection 
relating  to  gyro  gimbal  unbalance 

referring  to  hub  or  principal  reference  axis  system 

associated  with  horizontal  tail 

referring  to  inflow,  particle 

referring  to  inboard  feather  bearing  location 

spring  matrix  index 

associated  with  blade  attachment  joggle 
associated  with  gyro  end  of  feedback  rod  linkage 
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J 


associated  with  feedback  rod  coming  from  the  nth  blade 

generalized  mass  index 

associated  with  lead-lag  damper 

signifying  limiting  value 

blade  mode  index,  spring  matrix  index 


associated  with  main  rotor 


blade  number  index 


referring  to  blade  segment  neutral  axis 

newly  determined  value 

normal  (to  airflow)  component 

pertaining  to  nonrotating  value 

referring  to  outboard  feather  bearing  location 

value  from  previous  time  step 

associated  with  propeller;  perpendicular  blade  component 

referring  to  pitch  horn 

generalized  mass  index 

referring  to  rotor  axis  system 

associated  with  blade  feather  reference  value 

referring  to  control  gyro  feedback  lever  moment 

iferring  to  blade  spanwise  velocity;  general  mode;  static; 
structural;  shaft 

referring  to  blade  segment  shear  center 

referring  to  swashplate 

command  to  swashplate 

referring  to  swashplate  limit  stop 
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STEADY  steady  component 

SW  referring  to  blade  sweep  angle  location 

T associated  with  trajectory  path  relating  to  E axis; 

tangential  blade  component;  blade  torsion;  blade  twist 

TP  associated  with  the  tail  rotor 

TRIM  initial  or  trim  value 

TT  associated  with  tension  torsion  pack 

TTI  referring  to  inboard  end  of  tension  torsion  pack 

TTO  referring  to  outboard  end  of  tension  torsion  pack 

TW  associated  with  blade  twist  (built  in) 

UB  relating  to  control  gyro  unbalance 

UNSTEADY  associated  with  unsteady  component 

VT  associated  with  vertical  tail 

WING  associated  with  the  wing 

X relating  to  component  in  X direction 

Y relating  to  component  in  Y direction 

YA  relating  to  aerodynamic  component  in  Y direction 

Z relating  to  component  in  Z dijjection 

ZA  relating  to  aerodynamkf^omponent  ir  Y direction 

0 (nought)  associated  with  collective  value,  coordinate  axis 

value,  with  respect  to  principal  reference  axis,  blade 
root  summation 

1,2,3  with  respect  to  blade  modes  1,  2,  or  3 

IS  first  harmonic  component  shaft  axis  feathering 

lA  c with  respect  to  blade  1/U  chord 

3A  c with  respect  to  blade  3A  chord 


’ V i.’jfrwi-n—  - 


associated  with  the  feathering  mode  of  the  nth  blade 
relating  to  component  in  the  $ direction 
relating  to  component  in  the  0 direction 
relating  to  component  in  the  direction 


SUPERSCRIPTS 


referring  to  inertial  reference 

matrix  transpose 

(bar)  average  quantity 

(prime)  slope  with  respect  to  blade  span 

(dot)  time  derivative  of  basic  quantity 


(double  dot)  second  time  derivative 


matrix  inverse 


POSTSCRIPTS 


vector  quantity 


blade  radial  station  index 


blade  number  index 


